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Abstract 

We give a fairly complete analysis of outer billiards on the Penrose 
kite. Our analysis reveals that this 2 dimensional dynamical system 
has a 3-dimensional compactification, a certain polyhedron exchange 
map defined on the 3-torus, and that this 3-dimensional system ad- 
mits a renormalization scheme. The two features, the compactification 
and the renormalization scheme, allow us to make sharp statements 
concerning the distribution, large- and fine-scale geometry, and hid- 
den algebraic symmetry, of the orbits. One concrete result is that the 
union of the unbounded orbits has Hausdorff dimension 1. We estab- 
lish many of the results with computer-aided proofs that involve only 
integer arithmetic. 

1 Introduction 
1.1 Background 

Outer billiards is a dynamical system defined relative to a convex shape in 
the plane. B.H. Neumann [N] introduced outer billiards in the late 1950s, 
and J. Moser [Ml] popularized the system as a toy model for celestial me- 
chanics. See [Tl], [T2], and [DTI] for expositions of outer billiards and 
many references. 
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Figure 1.1: outer billiards relative to K. 

To define an outer billiards system, one starts with a bounded convex set 
K C and considers a point xq E B? — K. One defines xi to be the point 
such that the segment xqXi is tangent to K at its midpoint and K lies to the 
right of the ray xqXi. The iteration xq — )■ xi — )■ X2... is called the forwards 
outer billiards orbit of xq. It is defined for almost every point of E? - K. 
The backwards orbit is defined similarly. 

Outer billiards is an affinely natural system, in the sense that an affine 
map T : P Q carrying the convex shape P to the convex shape Q also 
carries the outer billiards orbits relative to P to the outer billiards orbit 
relative to Q. The reason is that affine maps carry line segments to line 
segments and respect the property of bisection. 

One of the central questions about outer billiards is the Moser- Neumann 
question, which asks if an outer billiards system can have unbounded orbits. 
Here is an abbreviated list of work on this problem. 

• J. Moser [M2] sketches a proof, inspired by KAM theory, that outer 
billiards on K has all bounded orbits provided that dK is at least 
smooth and positively curved. R. Douady [D] gives a complete proof. 

• In Vivaldi-Shaidenko [VS], Kolodziej [Ko], and Gutkin-Simanyi [GS], 
it is proved (each with different methods) that outer billiards on a 
quasirational polygon has all orbits bounded. This class of polygons 
includes polygons with rational vertices and regular polygons. In the 
rational case, all orbits are periodic. 

• In [T2] , Tabachniikov shows the existence of aperiodic orbits in the reg- 
ular pentagon case, and works out a renormalization scheme to explain 
their structure. 

• D. Genin [G] shows that all orbits are bounded for the outer billiards 
systems associated to trapezoids. See 
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In [SI], we settled Moser-Neumann question by showing that outer 
bilhards has some unbounded orbits when defined relative to K{(f)^'^). 
Here (j) is the golden ratio and K{A) denotes the kite with vertices 

(-1,0); (0,1) (0,-1) (AO); (1) 

Figure 1.1 shows an example. 

In [S2], we showed that outer billiards has unbounded orbits relative 
to K{A), when A is any irrational number in (0, 1). 



• Dolgopyat and Fayad [DF] showed that outer bilhards has unbounded 
orbits relative to the half-disk and other "caps" made from shcing a 
disk nearly in half. 

The shapes in [DF] and [S2] are the only ones known to produce unbounded 
orbits, though certainly it now seems that unboundedness is a common phe- 
nomenon. Our monograph [S2] gives more details about the history of the 
problem. 

The set R x Z^dd is invariant under the dynamics on K(A). We call the 
orbits that lie in this set special. In [S2] we gave quite a lot of information 
about special orbits on kites. For instance, we gave a formula for the Haus- 
dorff dimension of the set Ui{A) of unbounded special orbits, in terms of 
something akin to the continued fraction expansion of A. As a special case. 

We also showed that every unbounded special orbit is self-accumulating. This 
is to say that every point of an unbounded special orbit O is an accumulation 
point of O. 

In the 300 page [S2] we only considered the special orbits, for the sake 
of "brevity" . There is quite a bit more to say about the general orbits, and 
our purpose here is to say some of it, at least for the Penrose kite. The 
first phenomenon is that outer billiards on the Penrose kite, an unbounded 
2 dimensional system, has a 3 dimensional compactification. We saw similar 
things in [SI] and [S2]. 

The second phenomenon is that this higher dimensional compactification 
has a renormalization scheme. The renormalization scheme is the main new 
feature of this paper. Its existence allows us to get some precise results about 
the dynamics. We think that a similar scheme exists in great generality, but 
we currently don't have any techniques for investigating it in general. 
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1.2 The Distribution of Unbounded Orbits 



In all that follows, it goes without saying that our results concern only outer 
billiards on the Penrose kite. 

Theorem 1.1 Every orbit is either periodic or unbounded in both directions, 
and the union of unbounded orbits has Hausdorff dimension 1. 

It is convenient to state our remaining results in terms of the square of 
the outer billiards map, which we call ip. The map ip leaves invariant the set 



Rl = [j Rx{y + 2n). 



(3) 



Each Ry is a discrete countable family of horizontal lines. We always take y 
as a point of the circle R/2Z. In [SI] and [S2] we studied the orbits in R^. 

Let Z[(f)\ denote the ring of elements m + n0 where m,n E Z. We will 
sometimes use the notation 



m 
n 



m + ncj). 



(4) 



We define an equivalence relation on points of R/2Z. We say that a ~ 6 if 

b = ±(f>^''a + 2m + 2n(j); k,m,neZ. (5) 

We will explain this equivalence relation in a more natural way in §1.41 

An even length increasing sequence oi < ... < a2„ canonically defines 
a Cantor set C, as follows. Let be the similarity carrying [ai,a2n] to 
[a2k-i, a2k]- Then C is the limit set of the semigroup generated by Ti, ...,Tn. 
For instance, the sequence 0<l/3<2/3<l defines the usual middle-third 
Cantor set. We let denote the set obtained from C by removing the 
endpoints of all the complementary regions. 

Theorem 1.2 The set Ry contains unbounded orbits if and only if y ^ c 
and c G C*, where C is the Cantor set defined by the sequence 
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The set of such y is a dense set of Hausdorff dimension log(3) / log( 
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1.3 The Distribution of Periodic Orbits 



We define the winding number of a periodic orbit p to be lialf tlie number of 
times {ip^{p)} intersects tlie strip 

S = i2x[-2,2] (6) 

Tfiis definition makes sense geometrically. It turns out tliat, at least far from 
the origin, the ■j/'-orbits generally wind around the origin, nearly following a 
large octagon, and return to S after each half-revolution. See §4.21 and §4.21 
The following result says that orbits of high winding number are extremely 
pervasive. In particular, it says that has periodic orbits of arbitrarily 
high winding number provided that y ^ 2Z[(f)\. 

Theorem 1.3 For any integer N, there is a finite subset Wn C 2Z'[0] such 
that Ry has periodic orbits of winding number greater than N if y ^ Wn- 

As with any polygonal outer billiards system, every periodic point is 
contained in a maximal convex polygon consisting of points which all have 
the same period and combinatorial behavior. We call these maximal polygons 
periodic tiles. For convenience, we include K itself as a periodic tile. We call 
the union of the periodic tiles the dynamical tiling, even though it is only a 
tiling of a subset of the plane, and we denote it V. Theorem 1 1 . 1 1 say s that T> 
fills up everything but 1-dimensional set. 

A clean, finitary description of P is a beyond our reach. V is quite com- 
plicated. In particular, it contains infinitely many different shapes. However, 
we will get a near-complete understanding of the portion of V contained in 
the triangle T with vertices 

(0,1); (0,0-3); (0-\0). (7) 

T is bounded by lines extending 3 of the sides of K, as in Figure 1.2. We 
call T the fundamental triangle. 

Figure 1.2 below shows a tiling T of a full measure subset of T by an 
infinite union of kites and octagons, the kites are all similar to each other 
and the octagons are all similar to each other. The similarity factors all have 
the form (f)^'^, where (f) is the golden ratio and k is an integer. We will describe 
T precisely in §3.11 The fractal set which is the complement of the polygonal 
tiles has Hausdorff dimension 1. The set of points in this fractal set, having 
well defined orbits, also has Hausdorff dimension 1. 
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Figure 1.2: The triangle T and its tiling T. 
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Theorem 1.4 Suppose that p ^ T has a well defined orbit. Then p has a 
periodic orbit if and only if p lies in the interior of a tile ofT. Each tile of 
T is partitioned into finitely many periodic tiles. For any N, only finitely 
many of the tiles contain periodic orbits having winding number less than N. 

Remarks: 

(i) Theorem [13] is not quite as sharp as we would like. Experimental evidence 
suggests that each tile of T is itself a periodic tile. Indeed, we plotted Figure 
1.2 using our outer billiards program, which finds the periodic tiles. 

(ii) The set y such that the horizontal line of height y intersects T — T in an 
infinite set is precisely the Cantor set from Theorem 11.21 This is where the 
Cantor set comes from. 

(iii) The fractal set of aperiodic points for outer billiards on the regular 
pentagon has a self-similar structure akin to the one in Figure 1.2. In that 
case, the aperiodic orbits are all bounded. See Tabachnikov [T2] for details. 

1.4 Renormalization 

Now we reconsider the equivalence relation defined in connection with The- 
orem [L2l To say that a ~ 6 means that b = T{a), where T is an affine map 
defined over whose multiplying coefficient is a unit in Z[(j)\ and whose 
action is congruent to the identity mod 2Z[(j)]. Such maps form a group, 
and act on R/2Z with dense orbits. The equivalence classes we defined in 
connection with Theorem 11.21 are precisely the orbits of this group action. 

We think of G2 as related to the level 2 congruence subgroup of 
PSL2{Z), the modular group. In [S2] we discovered that, with respect to 
the special orbits on arbitrary kites, there is a kind of hidden r2-symmetry. 
For instance, the dimension of the set of unbounded orbits on Rf, as a func- 
tion of the kite parameter, is a r2-invariant function. We were not able to 
see the kind of renormalization structure that we establish here, but the way 
we think of things is that the union of all the dynamical systems defined by 
outer billiards on kites is a kind of plane-bundle over the parameter interval. 
We think there is a group r2 acting (in way that meaningfully relates to the 
dynamics) on this bundle in such a way that r2 gives the action on the base 
space and our group G2 here is the restriction of r2 to a particular fiber. 

Having indulged in some speculation, we now return to concrete results. 

Our renormalization works best for orbits which we call generic. Say that 
a point p = (x, y) E B? is generic if it does not satisfy any equation of the 
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form 

ax + by + c = 0; a,b,c & Z[(f)]; aj^O. (8) 

That is, (x, y) is generic if it does not lie on a non-horizontal line that is 
defined over Z[(j)\. The outer billiards map is entirely defined on the set of 
generic points, and preserves this set. So, it makes sense to speak of a generic 
orbit. Theorem 11.91 below explains the sense in which we do not miss much 
by ignorning the non-generic orbits. In §5.8l we discuss the issues surrounding 
the renormalization of non-generic orbits. 

We define two kinds of equivalence relations between orbits. 

• Let (O) denote the graph of an orbit O = {pn}, namely the subset 
{{n,pn)} C -R^. We call two orbits Oi and O2 coarsely equivalent if 
there is a bi-lipschitz map h: ^ such that h{{Oi)) and (O2) 
lie in /T-tubular neighborhoods of each other. We call K the coarse 
equivalence constant. 

• We call the orbits Oi and O2 are locally equivalent if, for each point 
Pi G Oi, there is a point p2 G O2, open disks Ai and A2, and a 
similarity S' : Ai — )■ A2 such that S{pi) = S{p2) and S conjugates the 
first return map ip\Ai to one of the two return maps ip^^\A2, at least 
on generic points. In particular, ^(Oin Ai) = 02 1"! A2. Also, S carries 
P n Ai to "D n A2 modulo the operation of subdividing each tile into 
finitely many smaller polygons. 

Theorem 1.5 Suppose that yi ~ y2 are two parameters in the same G2 or- 
bit. Then there is a bisection between a certain subset of the orbits in Ry_^ 
and a certain subset of the orbits in Ry^ . These subsets contain all generic 
unbounded orbits and also all generic periodic orbits having sufficiently high 
winding number. Corresponding orbits are both locally and coarsely equiva- 
lent, and the coarse equivalence constant only depends on {yi,y2) o-nd not on 
the individual orbit. 

The last statement in Theorem 11.51 is important for the periodic orbits. 
Every two periodic orbits are coarsely equivalent, so we need some kind of 
uniformity in order to make a meaningful statement. 

Here are two applications of Theorem 11.51 

Theorem 1.6 Every generic unbounded orbit is self- accumulating in at least 
one direction. 
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Theorem 11.61 is probably true for the non-generic unbounded orbits as well, 
but our techniques fall a bit short of this result. 

Theorem 1.7 Any generic unbounded orbit is locally and coarsely equivalent 
to a generic unbounded orbit that intersects the fundamental triangle T. In 
particular, in a small neighborhood of any generic point p with an unbounded 
orbit, the dynamical tiling in a neighborhood p is isometric to a neighborhood 
of T, modulo the addition or removal of a countable set of lines. 

The ambiguity concerning the countable set of lines comes from the set 
containing the non-generic orbit. The main thing that is missing in Theorem 
II. 71 is a description of the dynamical tiling in the neighborhood of points that 
do not have well-defined orbits. Figure 1.3 shows what the dynamical tiling 
looks like in a certain region whose lowest vertex is (3, 0), a point which turns 
out to be a fixed point of renormalization in a sense that we will make precise 
in our Fixed Point Theorem from 55. 




Figure 1.3: The dynamical tiling near (3,0). 
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1.5 Fine Points 



In this section we present a few more result that seem a bit more specialized 
than the ones above. 

For all r > 0, let Yr denote those y such that Ry has an unbounded 
orbit that intersects the disk of radius r about the kite vertex 0). 

Theorem 1.8 For all r > the set Yr is a nowhere dense set having Haus- 
dorff dimension log(3)/ log(0^). 

Comparing Theorem ll.Sl and Theorem ll.2l we can say qualitatively that most 
unbounded orbits stay far from the origin. In particular, for any compact 
subset, we can find an unbounded orbit that avoids this set. This is much 
different than what happens for the special orbits. For instance, our result 
[S2, Erratic Orbit Theorem] says in particular that every unbounded orbit 
in R1 returns infinitely often to every neighborhood of the kite vertex (0, 1). 

Our next result quantifies the sense in which there are fewer non-generic 
unbounded orbits than there are generic unbounded orbits. 

Theorem 1.9 The union of non- generic unbounded orbits has H aus dor ff di- 
mension at most log(3)/ log(0'^), and Ry contains non-generic unbounded 
orbits only if it contains generic unbounded orbits. 

Here is one more result about periodic orbits. 

Theorem 1.10 Suppose that y = m + ncf) with m and n odd integers. Then 
there is some open neighborhood V of the line R x {y} that contains no 
unbounded orbits. 

Theorem 11.101 says that these special lines are contained in the interior of 
the closure of the union of the periodic tiles. One might say that these lines 
are paved over with periodic tiles. 

Here is a nice special case of Theorem 11.21 

Theorem 1.11 Suppose that y = m + ncf) with m and n integers. Then Ry 
has unbounded orbits if and only if m is odd and n is even. 

We wonder which y G Q{4>\ are such that i?^ has unbounded orbits. 
Theorem 11.21 reduces this to an arithmetical question about the Cantor set 
C. 
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I. 6 Discussion 

Outer billiards is a non-compact 2- dimensional system. We exhibit a 3 dimen- 
sional compactification of (a certain first return map of) the outer billiards 
system. The compactification turns out to be a polyhedron exchange map 
defined on a 3-torus S. This result is very similar to the Master Picture 
Theorem in [S2]. We call our result the Compactification Theorem. 

Once we find the compactification, we will observe that it has a renormal- 
izable structure that is responsible for most of the results presented above. 
Specifically, we will find sets A, i? C S together with a 3 to 1 map R : A B 
such that R conjugates the first return map on A to the first return map on 
B. Our main result here is Theorem 15. 3[ the Renormalization Theorem. 

It is worth comparing the renormalization scheme here to the one [T2] 
for regular pentagons. In that case, one picks a certain bounded subset B 
of the plane and observes that the first return map to B is renormalizable. 
What this means, in part, is that there is a subset A G B and a similarity 
R : A ^ B which conjugates to fs- Here is the first return map to A 
and Jb is the first return map to B. It also means that the subsets A and B 
are large enough to capture all the dynamical properties of the whole map. 

A direct renormalization like this would be very difficult to establish in 
the presence of unbounded orbits, because the first return times to any suf- 
ficiently large compact set would be unbounded. Also, thanks to Theorem 

II. 81 there really is no compact set that "captures" all the relevant dynamics. 
What we do for the penrose kite is show the existence of the renormalization 
scheme on a higher dimensional compactification. Once we make the com- 
pactification, the renormalization is at least vaguely similar to what happens 
for the regular pentagon. 

One difference between the renormalization scheme here and the one for 
the regular pentagon is that the one here involves a 3-to-l covering map 
rather than a similarity. Another difference is that the regular pentagon case 
can be analyzed by hand, just following the orbits of several convex polygons. 
Here we need to keep track of about a million polyhedra just to see that the 
scheme works. This is what seems to make a computer-aided proof necessary. 

The moral of the story is that if one wants to find remormalization 
schemes for polygonal outer billiards, one should first compactify. In some 
sense, this is a lesson I learned from John Smillie. When hearing about my 
earlier work on outer billiards, Smillie guessed that probably there was a 
renormalization scheme behind it. 
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1.7 Overview of the Paper 

In §2 we study a certain 3-to-l covering map R : R/2Z — )■ R/2Z which 
governs the structure of our renormahzation map R. This map is closely 
related to the Cantor set in Theorem 11.21 The work is in §2 is a microcosm 
for the rest of the paper: We give some theoretical arguments to reduce the 
main result to a finite calculation which is too big to do by hand, and then 
we explain a rigorous computer calculation that finishes the proof. 

In §3 we study the tiling T shown in Figure 1.2 and prove a number of 
results about it. 

In §4 we present some background information about outer billiards and 
polyhedron exchange maps. Most of this information also appears in [SI] 
and [S2]. 

In §5 we state the Compactification Theorem, the Renormahzation The- 
orem, as well as several useful auxilliary results. One of these auxilliary 
results, the Fundamental Orbit Theorem, explains the structure of Figure 
1.2. Another result, the Fixed Point Theorem, explains the dilation sym- 
metry in Figure 1.3. The remaining results, the Near Reduction Theorem 
and the Far Reduction Theorem, explain the sense in which renormahzation 
"brings orbits closer to the origin" . 

In §6 we put together the material from §2-5 to deduce all the theorems 
mentioned in the introduction. We warn the reader that the order we prove 
these results is rather different from the way we have stated them. For 
instance, it takes almost every other result in order to prove Theorem 11.11 

In §7 we explain the main computational ideas we use in our proofs. These 
algorithms perform fairly standard tasks - e.g. detecting whether a point is 
contained in the interior of a polyhedron. Later chapters will refer back 
to the methods explained in §6. Given that the computational algorithms 
perform fairly standard tasks, the reader won't lose much understanding of 
the overall proof if they just skim the material in §7. On the other hand, we 
think that §7 might be very useful for someone who would like to reproduce 
our results of to prove similar results in a related stting. 

In §8 we prove a technical result, known as the Pinwheel Lemma, which is 
helpful in proving the Compactification Theorem. Versions of the Pinwheel 
Lemma also appear in [SI], [S2], and [S3], and related results appear in 
other works on the subject by other authors. 

In §9 we prove the Compactification Theorem. Our proof is very similar 
to the proof of [SI, Arithmetic Graph Lemma]. 
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In §10 we prove the Remormalization Theorem by reducing it to an ex- 
phcit calculation involving the polyhedron exchange map. 

In §11 we prove the Fundamental Orbit Theorem and the Fixed Point 
Theorem. 

In §12 we prove the Near Reduction Theorem. 
In §13 we prove the Far Reduction Theorem. 

In §14 we include coordinates for all the polyhedra involved in the poly- 
hedron exchange map, as well as coordinates the sets A and B. 

1.8 Computational Issues 

The general strategy of our paper is to reduce all the results to statements 
about finite partitions of various sets by polygons and polyhedra. These 
partitions sometimes involve a huge number of individual pieces, on the order 
of a million, and so it is necessary for us to use the computer to deal with 
them effectively. 

One source of potential error in a computer aided proof is floating-point 
(or roundoff) error. To avoid any possibility of floating-point errors, we 
perform our calculations using exact arithmetic in the number ring The 
special nature of our constructions allows us to do this. With exact arithmetic 
calculations, the one potential hazard is overflow error. We avoid overflow 
error by checking the sizes of the integers involved after every arithmetic 
operation. We describe the main features of these calculations in §7.11 

We illustrate our structural results in detail in §4, but there is probably no 
way for the reader to appreciate the details of the objects involved in without 
seeing explicit (and interactive) computer plots. We made an extremely 
detailed and extensive Java applet that lets the reader see everything in the 
paper. This applet is available on my website. For the reader who would 
like to do his/her own experiments, we include enough information in the 
appendix so that in principle one could reproduce the calculations. 

We would like to comment on some of the figures in the paper. To illus- 
trate certain definitions which make sense for any kite K{A), we will draw 
^^(1/4) in place of K{(f)~^), because ^^(1/4) is much easier to draw by hand. 
We draw these pictures mainly to give the reader a picture of what is going 
on, and for these purposes a picture of K{l/4) tells the whole story. Note 
that 1/4 — (f)^^ = .0139..., so the pictures we draw for K{l/A) are geometri- 
cally quite close to the ones for K{(j)~^). On the other hand, the computer 
pictures we draw will show K{(j)~^). 
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1.9 Further Results 



There are some other things I've noticed about outer bilhards on the Penrose 
kite. The other things have a different character from the results here. They 
have to do with the patterns one sees in the so-called arithmetic graph asso- 
ciated to the dynamics. In §4.51 1 briefly discuss one of the results, something 
I call the freezing phenomenon. 1 would like to have presented some of these 
results, this paper is already long enough and also I have not worked out 
proofs for these other results. 
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2 The Circle Renormalization Map 

In this chapter, we define and then study the renormahzation map R that 
we mentioned in the introduction. We also define and study the tihng T we 
mentioned in the introduction and plotted in Figure 1.2. The map R and 
the tihng TZ are closely related. 

2.1 Basic Definition 

Let 

T = R/2Z (9) 

Here we give a precise definition of the renormalization map R : T ^ T 
discussed in the introduction. 

We decompose [0, 2] into 5 intervals. 



• 


h 


= [o,r'] = [0,2-0]. 


• 


h 


= [0-2,20-2] = [2-0,1-0-3]. 


• 


h 


= [20-2, 2 - 20-2] = [1-0-3, 1 + 0-3; 


• 


h 


= [2 -20-2, 2 -0-2] = [1 + 0-3,0]. 


• 


h 


= [2-0-2,2] = [0,2]. 



We define R as follows. 



• liy eh then R{y) = 03y. Note that R{h) = /i U ... U h. 

• liy e h then R{y) = y + (jy - 0-2. Note that Rih) = h. 

• liy e h then R{y) = (j)^y -03 + 1. Note that Rih) = /i U ... U I5. 

• liy e h then R{y) ^y-(f) + 0-2. Note that Rih) = h- 

• liteh then R{y) = (p^y - 203 + 2. Note that Rih) = U ... U I5. 

The map it! pieces together correctly on the endpoints of these intervals, and 
induces a degree 3 covering map i? : T — )■ T. 

Note that R preserves the ring Z[(p], because 0^^ G Z[0] for all k E Z. 
Note also that the second iterate R^ is strictly expanding. 
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2.2 The First Descent Lemma 



As we mentioned above, R preserves the set 



A=Z[0]n[O,2] 



(10) 



The goal of this section is to prove the following number-theoretic result. 

Lemma 2.1 (Descent I) The action of R on A is the identity mod 2Z[(j)], 
and for any y & A, there is some n such that R^{y) is one of: 



0; 1; 0, 2-0; -! + </>, -3 + 30, 3-0, 5-30. 



Each collection of points, grouped according to the parity of their coefficients, 
forms a periodic cycle for R. 

One can see directly from the formulas that R is the identity mod 2Z^[0]. 
We concentrate on the proof of the second statement. 
Given y = m + n0 e A we define 



Lemma 2.2 Ify eAO {h U 74) then N{R{y)) < N{y) + 2. 

Proof: If y = m + n0 G /2 then R{y) = (m — 2) + (n + 2)0. In this case the 
result is obvious. If y G I4,, the proof is very similar. 4^ 

Lemma 2.3 Ify eAn {h U Is U I5) then N{R{y)) < N{y)/2 + 8. 
Proof: Suppose first that y E h. Let y — m + ncf). We compute 



Now we observe the following two inequalities. 

\m + 2n\ = |(m + 0n) - (0 - 2)n| < 2 + 0~^|n| < 2 + N{y)/2 

|2m + 3n| = I (2m + 20n) - (20 - 3)n| < 4 + 0-=^|n| < 4 + N{y) /2. 

The last inequalities come from the fact that y G [0, 2]. When y lies in 73 or 
I 5 the argument is the same, except that we have must either add 2 or 4 to 
our estimate to account for the translational part oi R. 4^ 



N{y) = max(|m|, |n|). 



(11) 



R{y) = (1 + 20) (m + n0) = (m + 2n) + (2m + 3n)0. 
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Corollary 2.4 For any y G [0,2] there is some n such that N{R"'{y)) < 20. 

Proof: Suppose that y & A has the property that N{y) > 20. Note that 
R{y) and R^{y) cannot both lie in I2 U 14. Hence, we may combine our last 
two results to establish the bound N{R'^{y)) < y/2 + 10 < N{y). So R^^y) 
has smaller integer norm than does y. 4|k 

The above work reduces the proof of the Descent Lemma I to a finite 
calculation of what happens to those y for which N{y) < 20. We make the 
finite calculation and see that the Descent Lemma I holds for these values 
as well. Indeed, in the next section, we will describe a much more extensive 
calculation. 

We mention a variant which we will need in one place. Let S be the 
involution S{y) = 2 — y. 

Lemma 2.5 Let y G Z[(f)] and let y' = Ri o ... o i?jv(y), where each Rk is 
either R or Ro S . Then y' is one of the values listed in the Descent Lemma 
I provided that N is sufficiently large. 

Proof: The map R commutes with S and all the cycles listed in the Descent 
Lemma I are invariant under S. (The first cycle is really [0] = [2].) The 
result follows immediately from these two facts. 4 



2.3 The Second Descent Lemma 

Let G2 be the affine group defined in §1.41 This group has two components, 
so to speak. We let denote the index 2 subgroup consisting of maps of 
the form 

T(x) = (p^'^x + b; keZ; be 2Z[0]. (12) 

This formula differs from the one in Equation [5] only in that we do not allow 
a minus sign in front of (j)^''. 

Lemma 2.6 (Descent II) Two elements yi,y2 G R/2Z lie in the same 
orbit iff there are positive integers ui and n2 such that i?"^(?/i) = i?"^(?/2)- 
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Proof: All the maps defining R lie in G^. So, if there are positive integers 
rii and n2 such that R^^iyi) = R"'^{y2) then yi and t/2 he in the same G2 
orbit. The converse is the interesting direction. 

For the converse, suppose that yi and 2/2 lie in the same orbit. Con- 
sidering the action of R, we can replace 7/2 by some image y'2 = R^{y2) such 
that ?/i — 1/2 ^ 2Z[0]. So, without loss of generlity, we can consider the case 
where we already know that yi — y2 G 2Z[(j)]. The Descent lemma I proves 
this result whenever our points lie in Z[(j)]. So, it suffices to consider the case 
when neither point belongs to 2Z'[0]. 

We find it more convenient to work with a new map that is closely related 
to R. We define p = R on Ii U U and p = R^ on I2 U I4. Unlike R, 
which pieces together continuously across the endpoints of the intervals in 
the partition, the map p is not defined on the endpoints. However, the points 
we are considering, and their orbits, never hit these endpoints. It suffices to 
prove this result for p in place of R. 

To describe the map p, we let /[m, n] denote the map 

X — )• (p^x + 171 + ncj). (13) 

Then 



• 


On 


h, 


we 


have 


P = 


/[o,o]. 


• 


On 


h, 


we 


have 


P = 


/[2,-2]. 


• 


On 


h, 


we 


have 


P = 


/[0,-2]. 


• 


On 


h, 


we 


have 


P = 


/[-2,-2] 


• 


On 


h, 


we 


have 


P = 


/[0,-4]. 



This is a short calculation which (after many tries) we did correctly. We omit 
the details. 

We write (2/1,1/2) = max(|m|, |n|), where yi — y2 = + ncj). The same 
argument as in Lemma 12.31 and Corollary 12.41 shows (with tons of room to 
spare) that 

(2/1,^2) > 100 =^ {y[,y'2) < (2/1,2/2) - 1. (14) 

Here we have set 2/^ = piVk)- Equation [H] reduces this lemma to the case 
when (2/1,2/2) < 100. In the next section, we will explain our computer- 
assisted proof that The Descent Lemma holds for such choices. 4 
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2.4 A Dynamical Computation 

Notice that our statement that the Descent Lemma II holds for pairs 1/2) 
with {yi,y2) < 100 is not obviously a finite calculation, because it involves 
infinitely many values. However, we will explain how to reduce the problem 
to a finite calculation, which we then make. This situation is typical of the 
results in this paper. The challenge is reduce seemingly infinite statements 
to finite computations. 

Let Q — [0,2]^. The set of pairs (^1,^/2) of interest to us lie on a finite 
number of line segments of slope 1 that are contained in Q. By switching 
the order of the two points if necessary, we can assume that yi < y2- We let 
p act on Q by having p act separately on each coordinate. 

The square Q is partitioned into 25 subsquares 

Qij ^liX Ij (15) 

on which p is entirely defined and a similarity. Our point of view is that p acts 
separately on each Qij, and the action on the various boundaries depends on 
which square we include the boundary in. Anyway, we don't care about what 
happens on the boundaries: As we said above, the points we consider, and 
their orbits, never hit the boundaries. 

Say that a diagonal is a segment of slope 1 contained in Q. We call the 
diagonal small if it is contained in one of the 25 subsquares. The endpoints 
of a small diagonal might lie in the boundary of the subsquare, but this is 
fine with us. Given a small diagonal /, we can define p(/) using the action 
of p\Qij. Note that p{I) is another diagonal, but not necessarily a small one. 

Each diagonal has a canonical decomposition into small diagonals: We 
just take the intersections with the 25 sub-squares. This we have a kind of 
dynamical system defined on lists of diagonals: Given a list of diagonals, we 
first subdivide each member of the list into small diagonals. Next, we let p act 
on all the small diagonals. And so on. As one final nicety, we switch the two 
coordinates of the endpoints, if necessary, so that all our diagonals lie about 
the line y2 — yi — 0. (We do this simply for computational convenience.) 

Say that a good seed is a diagonal of the form 

A{m,n) = Qn{y2-yi^2m + 2n(f)}, (16) 

where 

2m + 2n0 e [0,2]; max(|2m|, |2n|) < 100. (17) 
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We run the dynamical system starting with any good seed and we find that, 
after finitely many steps, the only remaining intervals lie in 

A(0,0) U A(-2,2) U A(4,-2). (18) 

Points (?/i,?/2) e A(0,0) obviously satisfy yi = y2. When (?/i,?/2) e A(-2,2), 
it means that y2 — Vi = <P ~ Iii this case, we either have yi G I2 and 
?/2 = R{yi) or y2 G I4 and yi = R{y2)- Compare the definition of the map R. 

It only remains to deal with those points in A(4, —2). For this purpose, 
we just have analyze the dynamics more carefully. Define 

A'(4, -2) = A(4, -2) n g24. (19) 

When we perform the dynamics on A (4, —2) we find that the following occurs. 

• A(4, —2) breaks up into 5 small diagonals, one of which is A' (4, —2). 

• p maps each of the 4 other small diagonals into A(0, 0) U A (—2, 2). 

• p maps A'(4, -2) back into A(4, -2). 

This analysis shows that the Renormalization Lemma can only fail for a pair 
of points (2/1,2/2) such that p"'{yi) G I2 for all n (and also ^"(1/2) G /4 for all 
n.) But the fixed point of Pnlh is an endpoint of I2 and p is an expanding 
map. Since yi is not this endpoint, we see that p"(?/i) eventually escapes I2, 
and we are done. 

Remarks: 

(i) We perform the calculations with exact arithmetic, as explained in §7.11 

(ii) In §7.21 we explain how we eliminate any possibility of overfiow error in 
our calculations. Even without specific guards against overfiow error (which 
we do have) for all our seeds the dynamical system reaches the 3 end-states 
above very quickly and all integers remain pretty small. 



2.5 The Cantor Set 

In this section, we give some information about the Cantor set C from The- 
orem 11.21 First of all, the main property of C is that both C and C* are 
forward i?-invariant. That is, R{C) = C and R{C'^) = C*. Indeed, this is 
how we discovered the map R. 
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Lemma 2.7 Let y e -^[0]- Then y e only if y = 1 mod 2Z[0] anc? 
y & C — C* on/y if y — m -\- n(j) with m even. 

Proof: Of the 8 values listed in the conclusion of the Descent Lemma I, 
we see that only 1 lies in C*. Given any y e Z[(j)] which intersects C'^, we 
simply note that G C* for all n. But there is some n such that 

is one of the 8 values hsted in the Descent Lemma L But this means that 
R"'{y) = 1. Hence y ~ 1 mod G^. But this means that y = 1 mod 2Z[(f)\. 

For the second statement, we just have to rule out the case that y = —l+(f) 
mod 2Z[(p]. But we check easily that —1 + ^ C. But C is forward R- 
invariant and = —1 + for some n. This situation is impossible. 4 

To each point y G R/2Z we assign a renormalization sequence in the 
digits {1, 2, 3, 4, 5}. The sequence is such that the nth term is k if and only 
if R^iy) e Ik- Not every point has a unique renormalization sequence. A 
point y has a non-unique renormalization sequence if and only if R^{y) G dik 
for some n and some k. Examining the endpoints of our intervals and also 
the conclusion of the Descent Lemma, we see that this happens if and only 
if y = m + n0, where m and n are integers and n is even. In particular, all 
points of have unique renormalization sequences. 

Lemma 2.8 A point lies in C* if and only if its renormalization sequence 
is unique and has no 5 '5 in it. 

Proof: A point y has a non-unique renormalization sequence if and only 
if R"'{y) G dIk for some n and some k. Examining the endpoints of our 
intervals and also the conclusion of the Descent Lemma, we see that this 
happens if and only if y = m + n0, where m and n are integers and n is even. 
In particular, all points of have unique renormalization sequences. 

Note that the renormalization sequence of a point in y G cannot start 
with 5, because = [4>, 2] only shares its bottom endpoint with C, and this 
endpoint lies in C — C*. Since is forward invariant, we see that the 
renormalization sequence cannot have any 5's in it at all. 

Now we know that C* only contains points that have unique renormal- 
ization sequences with no 5's in them. Conversely, suppose y has a unique 
renormalization without 5's. Suppose that y ^ C. Since y ^ I^, we can say 
that y lies on one of the bounded components of i? — C. These components 
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all have diameter 0-3fc+i k = 1,2, 3... The largest component has size 
and is precisely the interior of J2. But y ^ I2 because then R{y) G I5. But, 
iteration of R increases the sizes of all gaps except the largest one. Hence 
R^iy) lies in the largest gap for some n. This is a contradiction. 

The endpoints of the gaps do not have unique renormalization sequences. 
So, the same argument rules out the possibihty that y e C — C*. We con- 
clude that y e C*. 4 
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3 The Fundamental Tiling 



3.1 Definition of the Tiling 

Let T be the tiling of the fundamental triangle T shown in Figure 1.2. In 
this section we define T precisely. As we mentioned in the introduction, T 
is bounded by three sides of the Penrose kite K. The reader can most easily 
understand the definitions we make by referring back to Figure 1.2. 

Let / : — ?■ be the map that fixes the top vertex of T and shrinks 
distances by a factor of 0^. We define Kq = f{K), where K is the Penrose 
kite. Kq is the largest kite in the tiling T. The largest octagon Jq in T has 
4-fold dihedral symmetry, and 3 vertices located at 

(13-80,4-20) (5 -30,6-30) (-3 + 20,-2 + 20). (20) 

This is enough information to characterize Jo uniquely. A calculation shows 
that 

T-Ko-Jo = (Tn U T^2) U Tgi U T32 U T41 (21) 

Here is a similar copy of T, with scaling factor 0~^. Our notation is 
such that Tjk C Ij, the interval used in the definition of the map R. Just 
to pin thing down exactly, we say that the centroid of Tn lies to the left of 
the centroid of Tj2. The bottom triangles Tn and T12, mirror images of each 
other, are not disjoint. However, the similarities carrying T to Tu and T12 
both map Jq to the same smaller octagon. Therefore, we can compatibly 
subdivide each of our smaller triangles into a kite, an octagon, and a union 
of 5 smaller triangles. Continuing this process forever, we get the tiling T. 

Let S denote the union of points in T that do not belong to the interiors 
of any of the tiles of T. Let S'* denote those points of S that do not belong 
to the boundarty of and tile of T. Then S = S** is contained in a countable 
union of lines defined over Z [0] . These are the lines extending the sides of 
the boundaries of the tiles in T. We call S"* the fundamental fractal. We 
have the basic relation 

C* = 1T2{S*y, C = 7T2{S). (22) 

Here it2 is projection onto the second coordinate, and C is the Cantor set 
from Theorem 11.21 
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3.2 HausdorfF Dimension 

The purpose of this section is to prove the following result. Here dim(S') 
refers to the Hausdorff dimension of 5*. See [F] for details about Hausdorff 
dimension. 

Lemma 3.1 (Dimension) Let S he the fundamental fractal and let L he 
any countahle collection of non-horizontal lines. Then dim(S' — L) = 1. In 
particular, dim(S'*) = 1. 

Proof: After we remove Jq and Kq from T we are left with the union in 
Equation [211 We have the equation 

Tn U Ti2 - R^^ (i^o) = Tn U , (23) 

where is similar to T, with similarity factor Figure 3.1 shows this 
operation. 




Figure 3.1: (Tn U T12) - R^i{K^) = Tn U T[^. 

So, removing two open kites and an open octagon from T leaves the 5 
(temporarily renamed) disjoint triangles Ti, ...,T5. We have 

5 

diam(Tfc) = 0-^(4 + 0-^)diam(T) = diam(T) (24) 

fc=i 

The set S is the self-similar fractal which is the limit set of the semigroup 
generated by the similarities which carry T to each of the smaller triangles. 
From here it is an exercise to show that dim(S') = 1. 

Note that C, the Cantor set from Theorem 11.21 is the set y G [0, 2] such 
that the line of height y intersects 5*. Each line in L intersects in a set 
of dimension dim(C) < 1. Since dim(S') = 1 and dim(S' fl T) < 1, we have 
dim(^- L) = 1. ^ 
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3.3 The Horizontal Intersections 



Let C* be as in the previous chapter. Every y G C* is such that 



Ly = Rx {y} 



(25) 



intersects 5** nontrivially. The goal of this chapter is to prove the following 
result. 

Lemma 3.2 (Horizontal) For any y G C*, the set LyCiSisa Cantor set. 
Hence Ly fl 5* is obtained from a Cantor set by removing at most countably 
many points. 

We begin by explaining how the set 5* may be constructed recursively, as 
the nested intersection of finite unions of (overlapping) triangles. 

Recall that T — Jq — Kq = [jTij, a union of 5 triangles. Let pij be the 
similarity such that pijiTij) = T. Define Sq = T and (inductively) Sn such 
that 



Our proof involves an analysis of how these triangles sit with respect to the 
horizontal lines. 

Remarks: 

(i) Note the similarity between the maps pij and the map p considered in 
the previous chapter, precisely, pij acts on the horizontals having heights in 
Jj exactly as p acts on Jj. We will pursue this analogy further in the next 
section. 

(ii) We have to be careful with our definition, because Tn and T12 overlap. 
Here is the justification for what we do. Assume by induction that Sn-i is 
well-defined and that Sn-i H (Tn fl T12) has bilateral symmetry and that pn 
and pi2 have the same action on this intersection. Then Sn is well-defined 
and inherits all these same properties by symmetry. 

Let {kn} be the renormalization sequence associated to y. By Lemma 
12.81 this sequence and has no 5's. Since Sn is a finite union of triangles, the 
intersection 



Pij{Sn n Tjj) — 5, 

Then Sn is a finite union of triangles and 



'n-l- 



(26) 




(27) 
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Kn = Lyn Sn (28) 
is a finite union of disjoint intervals. 

Lemma 3.3 Let J he a maximal interval of An- Then JnA„+i is nonempty 
and contains 2 maximal intervals in case kn = S or {kn, kn+i) = (1,4). 

Proof: Let yo = y and ?/„ = R^iy) and L{n) = Ly^. Our sequence 
starts witfi /cq. We fiave yn G Ik„- 

Recall that 5*0 = T. It follows from induction and fact that pijiTij) = T 
that there is a sequence of maps po, ...,p„_i such that 

Pn-io ... o pq{J n An+m) = L{n) n Sm] m = 0, 1,... (29) 

Here pj = Pkj,mj, where nij G {1,2}. The sequence {pi} is not necessarily 
unique, because Tn and T12 overlap. We don't mind this. 

Since /c„ G {1,3,4}, the set on the right hand side of Equation [29] is a 
nontrivial union of intervals. Hence J fl A„+i has at least one interval. When 
kn = 3 the set L{n) fl Si contains 2 intervals. Hence J fl A^+i contains 2 
intervals in this case. When kn = ^ and kn+i = 4 the line L{n) lies above the 
top vertex of Tn nTi2. Hence L{n) (iSi again consists of 2 intervals. Hence, 
so does J n A„+i. ^ 

Corollary 3.4 Let J be a maximal interval of An- Then JClAm contains at 
least 2 disjoint intervals for m sufficiently large. 

Proof: The sequence associated to y cannot terminate in an infinite string 
of I's. Otherwise, there is another sequence associated to y which terminates 
in an infinite string of 5's. So, the associated sequence either has an infinite 
number of 3's or an infinite number of (l,4)'s. The corollary now follows 
immediately from Lemma 13.31 4ft 

It follows from Corollary 13 .41 that the nested intersection A„ is a Cantor 
set. This completes the proof of the Horizontal Lemma. 
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3.4 Notation 



We have already mentioned the similarities pij. These maps have the prop- 
erty that pij{Tij) — T. These maps would seem to suit us perfectly well, but 

it turns out that there is a slightly more elaborate collection of maps that are 
better adapted to the structure of outer billiards. We introduce these maps 
here. These maps have the advantage that they are all homotheties. 
Define the reflection in the vertical line x — 1: 

r){x,y)^{2-x,y). (30) 

We first change our notation a bit. We let T+ — T and — T, etc. That is, 
we attach the (+) superscript to all the objects associated to the fundamental 
triangle. Next, we define X~ = ?7(X+) for any object X. 

The two triangles and T~ are mirror images of each other. We will 
see eventually that the dynamical tiling intersects T~ exactly in the tiling 
T" . It turns out that the dynamics on 7+ U T~ works out more nicely than 
the dynamics on just T^, even though ultimately all our results are phrased 
just in terms of T"*". 

We define the following ten maps. 

• i?+ = r)n : T+ ^ T+. 

. Rt, = vopii-T,\^T-. 

. R+^riop^,:T+^T-. 

• = V32 : T+ ^ r+. 

• R'li is the isometry carrying T^^ to Tj^. 

• Rio ^v°Rtj°v- 

The maps R^j acts similarly to the map i?^. For instance R^i{T^i) — T'^. 
Recall that the renormalization map R equals the map Ri on the interval 
Jj. The maps Rfj are all similarities preserve the horizontal foliation and act 

, , , I 

on the horizontal lines intersecting T^- as Ri acts on /j. In particular, the 
maps i?4i are better adapted to R4 than the map p4i, which is really adapted 
to the map R\o R^^. 

The maps are orientation preserving isometrics. The remaining 8 
maps are orientation reversing similarities, with expansion constaint (j)^. 
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3.5 The Renormalization Set 

Given subsets A,B cT+UT-,we write A ^ B if A C T^^ and B = 
We use this definition in particular for points. Suppose tliat (pi, qi) is a pair 
of points, botfi at tlie same fieight. We write — )■ (P2,5'2) if Pi — ^ P2 
and (ji — )■ g2- We write p ~ g if 

(p, g) ... (p', g'); p' = q'. (31) 

Define 

T(p) = {g|g~p}. (32) 

Note that T(p) consists of points that are all on the same horizontal level as 
p. We call T(p) the renormalization set of p. 

Remark: In §6.61 we will see that p and q lie in the same outer billiards 
orbit provided that p and q both have unbounded orbits and p ^ q. This is 
a step in our proof that generic unbounded orbits are self-accumulating. 

Let be the set studied in the previous chapter. The goal of this section 
is to prove the following result. 

Lemma 3.5 (Density) Let = L fl , where L is a horizontal line 
whose height lies in C*. For any point p G A+ U A^, the set T(p) is dense 
in A+ U A~. 

the Density Lemma. We call and T~ the distinguished triangles of 
depth 0. Recall that is contained in the nested intersection of sets 5"^. 
We say that a distinguished triangle is a maximal triangle of S^. There 
are 2 x 5*^ distintuished triangles. We call n the depth of the distinguished 
triangle. 

We say that a distinguished triangle is related to a renormalization set 
if it intersects the horizontal line containing the renormalization set. Let 
P{n) be the property that every renormalization set intersects each related 
distinguished triangle depth n. It suffices to prove that the statement P(n) 
is true for all n. 

If r is any distinguished triangle of depth n, then r — )■ r' — )■ r" where one 
of r' or t" has depth n — 1. From this, and from the definitions, we see that 
P{n — 1) implies P{n). It just remains to establish P(0). 
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Say that points p,q & T+ U T~ are distantly placed if one of the points 
hes in T"*" and the other hes in T^. Say that y G is good if -P(O) holds 
for all renormalizations sets of height y. 

Lemma 3.6 Suppose that yi G R^^{y2) and j/2 is good. Then yi is good. 

Proof: Let pi be some point having height yi. We have pi — )■ p2 for some 
P2 having height y2. Since 1/2 is good, there is some q2 such that p2 ~ q2 and 
P2,q2 are distantly placed. Without loss of generality, assume that p2 G T+. 
There is a depth 1 distinguished triangle r"*" such that r+ — )■ and pi G r+. 
Let r~ = p(t+). Then — )► and points in r~ are distantly placed from 
points in r+. In particular, we can find qi G t~' such that qi — ?■ q'2- But then 
Pi and gi are distantly placed and pi ~ gi . Since pi was chosen arbitrarily, 
yi is good. 4 

Lemma 3.7 Suppose that y G fl Is. T/ien ?/ is good. 

Proof: Any horizontal line having a height in I3 intersects the 4 disjoint 
triangles for j = 1,2. Let p be the union of the 4 special maps. Then 

p(T3+) = p(T32) = T+; p(T3+) = p(T3l) = T^. 

It follows from this equation that every renormalization set of height y either 
intersects both of (T3"^,T3"^) or both of (T3"^,T3"^). 4 

Lemma 3.8 Suppose that y G C* fl Ji fl R^^^I^). Then y is good. 

Proof: Any horizontal line having a height in I3 intersects the 4 triangles 
for j = 1,2, in disjoint intervals. (This is true even though the trian- 
gles themselves are not disjoint.) The rest of the proof is as in Lemma [3?71 ^ 

Let y G be arbitrary. As in the proof of Corollary 13. 4[ there must be 
some n such that R^{y) satisfies either Lemma [3.71 or Lemma [3.81 But then 
R^iy) is good. But then y is good as well. Hence, -P(O) holds. But then 
P{n) holds for all n. This completes the proof of the Density Lemma. 
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4 Preliminaries 



4.1 Polytope Exchange Maps 

Definition: For us, a polytope exchange map is a quadruple {M, Xi, X2, 
where M is a flat manifold (possibly with boundary), Xi and X2 are locally 
finite partitions of M into convex polytopes, and / : M ^ M is a piecewise 
isometric bijection which carries X^ to X2. We mean that / is a translation 
when restricted to each polytope P of Xi and f{P) is a polytope of X2. 
Technically, / is not defined on the boundaries of the polytopes. When M 
is a compact manifold, we require that the partitions be finite. 

As a special case, suppose that M = R^/A, where A C is a discrete 
group of translations. We say that a polytope in M is golden if any lift 
to has all vertices with coordinates in the ring Z[(f)]. This definition is 
independent of lift. We call an associated polytope exchange map golden if 
all the polytopes are golden, and if all the translations are defined by vectors 
in {Z[(p])"'. In this paper we will consider two golden polytope exchange 
maps. 

• A 2 dimensional non-compact polygon exchange whose domain is 
the infinite strip E = x [—2, 2]. 

• A 3 timensional compact polyhedron exchange ^, whose domain is the 
torus S = {R/2Zf. 

Fibered Polyhedron Exchange Maps: Let "H be the foliation of S by 
horizontal 2-tori - those obtained by holding the third coordinate constant. 
We say that a leaf of H is golden if its height lies in Z[(p]. We call a golden 
polyhedron exchange map on E fibered if it preserves the leaves of T-L and if 
the restriction to each golden leaf is a golden polygon exchange map. The 
second condition docs not follow automatically from the first; there is an 
auxilliary condition we need on the edges of the polyhedron partition, as we 
now explain. 

Let e be a non-horizontal edge of some polyhedron in the partition. Let 
{xi, Hi, Zi) and {x2, y2, -^2) be the endpoints of e. The 4 quantities 

X2 -xi y2 - vi X2Z1 - X1Z2 y2Zi - yiZ2 , . 

Z2 — Zi Z2 — Zi Z2 — Z\ Z2 — Zi 

all belong on Z[(f)] if and only if every nontrivial intersection of e and a golden 
leaf has coordinates in Z[(f)]. We omit the easy proof of this result. 
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Periodic Tiles: Let {M, Xi, X2, /) be a polytope exchange map, not neces- 
sarily a golden one. Given a periodic point p G M, there is a maximal convex 
polytope Dp consisting of points q such that p and q have the same period 
and the same itinerary. By itinerary we mean the sequence {uk}, defined so 
that the (nfc)th polytope P{nk) of Xi contains f^{p). We have 



Being the finite intersection of convex polytopes, Dp is also a convex polytope. 
We call Dp a periodic tile. We call the union of all the periodic tiles the 
dynamical tiling. When all points of M are periodic, the dynamical tiling is 
another partition of M. 

Our terminology is slightly misleading. In general, the dynamical tiling 
only fills a subset of M. This subset need not even be dense in M. A nice ex- 
ample is furnished by outer billiards on a generic trapezoid. See [G], where 
the issue is phrased somewhat differently. In our examples, thanks to the 
results in this paper, the dynamical tilings are dense. 

Homogeneity: As one might expect, the dynamical tiling has some nice 
homogeneity properties. We say that two points Xi.,X2 G M are tile isomet- 
ric if there are open disks Ai and A2 such that 

• Xj G Aj for j = 1, 2. 

• 2? n Ai is isometric to V D A2. 

Note that Xj need not be the center of Aj and the isometry from from "DPI Ai 
to A2 need not carry Xi to X2. We caU a subset S <Z M tile homogeneous 
if every pair of points in M are tile isometric. 

As a related definition, we say that xi,X2 G M are tile similar if the 
above properties hold, with the word similar replacing the word isometric. 

Lemma 4.1 Every orbit is tile homogeneous. 

Proof: Let po and Pn = f^{po) be two points in the orbit. Let g = f"'. The 
map g is an isometry in a neighborhood of po and g conjugates / to itself. It 
follows from this fact that po and p„ are tile isometric. d|k 




(34) 



k 
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4.2 The Return Map 

The square outer billiards map ip is a local translation. For points far from 
the origin, the dynamics of ■0 is quite simple, forward iterates of ip generally 
circulate counterclockwise around the kite K. Near the origin, the dynamics 
of ip is slightly more complicated. Our paper [S3] considers these dynamics 
carefully, for fairly general convex polygons. 




Figure 4.1: The second return map far from the origin 



The set K' shown in Figure 4.1 is a large compact set that contains K 
well in its interior. Given that the V'-orbits generally circulate around the 
kiet K, at least far from K, it makes sense to consider the return map to 
each half of a suitable horizontal strip. Define 

E = i2x[-2,2]. (35) 
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half this strip is shaded in Figure 4.1. What makes this strip canonical is 
that, far from K and near the consecutive iterates of -0 differ by the 

vector (0, ±4). So, E has just the right width. 

Remcirk: Sometimes we will want to leave off the bottom boundary of E 
and sometimes we won't. It turns out that the first return map is the identity 
on the boundary of E, so we can essentially just ignore these points. 

The definition of the return map to "each half" of E presents some prob- 
lems for points near K. There are "bad points" of E that are too close to K. 
These bad points are in the regions labelled B in Figure 4.2. 




Figure 4.2: A decomposition of E into regions. 

The problem with the bad points is that one cannot really say which side 
of K they are on. We mention this problem in order to justify a modified 
version of the return map, which we now define. Let E_ and E+ be the two 
components of E — — 5. Let ^-t : E-t — ?> E± to be the first return map of 
ip to E±. Just to be clear, is the first return to E+ and is the first 
return to E_. We define ^ — ^+ U ^_ to be the "union map", defined on 
E+ U E_. Finally, to get a map on all of E, we define ^ to be the identity 
on K U B. Note that this definition does not correspond with the action of 
outer billiards on B, but we will deal specially with the points in B whenever 
necessary. 

We call ^ : E — >■ E the return map. 
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Lemma 4.2 (Return) The following is true. 

1. Every ip orbit intersects S+ U S_. 

2. \1/ is well defined on all points of S that have a well-defined orbit. 

3. There is some C such that — p\ < C for all p ^ S with a well 
defined orbit. 

Proof: First of all, we proved the same result in [S2, §2.3], in the context 
of special orbits (on R^) for arbitrary kites. The proof there works in this 
setting with only minor changes. 

Here we explain a different proof. The return Lemma is obvious for points 
of the form {x,y) with |x| > 20. These points just circulate around the kite, 
nearly following a giant octagon before coming back to the strip. For points 
in the strip S having |x| < 20, the calculation we make in connection with 
the Pinwheel Lemma in §8.21 in particular establishes the Return Lemma. 
The calculation we make there simply involves covering [—20, 20] x [—2, 2] 
with 572 convex polygons such that the first return map exists and is well 
defined on the interior of each tile. We also prove that the remaining points, 
the ones in the boundaries of our tiles, do not have well defined orbits. 4|k 

The Return Lemma allows us to work with rather than ip. We state 
the following result in terms of the action of on S+. The same result 
holds with (— ) in place of (+). 

Lemma 4.3 There is a canonical bijection between the unbounded ip orbits 
and the unbounded \E'+ orbits. The bijection is such that a given \E' orbit cor- 
responds to the ip orbit that contains it. Two orbits are coarse equivalent 
if and only if the corresponding ip orbits are coarse equivalent. 

Proof: We associate to each orbit the unique ip orbit that contains it. 
This injective association is also surjective, by Statement 1 of the Return 
Lemma. Given our description of \1/ for points far from the origin, it is clear 
that we can reconstruct the coarse equivalence class of a ■j/'-orbit from the 
coarse equivalence class (defined the same way) for the corresponding 
orbit. Up to a uniformly bounded error, we obtain the graph of the ip orbit 
from the graph of the \E'^ orbit by attaching a centrally symmetric octagon 
of a suitable radius to each point of F that is sufficiently far from the origin. ^ 
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4.3 An Unboundedness Criterion 

Here we establish a useful criterion for the unboundedness of the orbits. 
Let Z+ denote the set of positive integers and let Z_ denote the set of 
negative integers. We say that a subset S C Z is uniformly dense in Z+ if 
there is some N such that every point of Z^ is within N units of Z_^^. We 
make the same definition relative to Z . 

Lemma 4.4 Let O be an infinite ^ orbit. Suppose that there is an open 
horizontal line segment S such that 0(1 S is a nonempty and nowhere dense 
subset of S. Then O is unbounded in both directions. 

Proof: We will assume that O is bounded in the forward direction and 
derive a contradiction. Let A = (j)^^ . Given the locations of the vertices of 
the Penrose kite, we have the following formula. 

^(p) = (2niy4 + 2n2,2n3); m + ns + ns = mod 2. (36) 

When p is far from the origin, the orbit of p stays within a uniformly thin 
tubular neighborhood of a centrally symmetric octagon, as we mentioned in 
connection with Statement 3 of the Return Lemma. The sides of this octagon, 
which depends on p, are always integer multiples of the vectors listed in §4.21 
Moreover, opposite sides have the same length. For this reason, the vectors 
entering into the sum that defines '^{p) — p nearly cancel in pairs, and we 
find that there is a uniform bound to max(|?T,i|, |n2|) in Equaton 

Let (xo,?/o) be some point of O. We can find integers (a„,6„) such that 

^"(a;o, Z/o) = {xq, yo) + 2a„A + 26„. (37) 

The differences |a„+i — a„| and |6„+i — bn\ are uniformly bounded. 

The sequence fl = {onA + bn} is both infinite and bounded. Hence, 
^1 = {^n} has infinitely many values. Our uniform bound on |a„+i — a„| now 
implies that Qi is uniformly dense in at least one of or Z . Combining 
this fact with the fact that Q is bounded, we see that there is some such 
that the union 

n*= \J n + (38) 

i,j<N 

either contains every integer combination of the form aA + 6 G (0, 1) with 
a > 0, or every such integer combination with a < 0. In either case, f2* is 
dense in (0, 1). 
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But f2* is a finite union of translates of f2. We have shown that a finite 
union of translates of Q is dense in (0, 1). But the set Xq + 2f2 is a subset 
of O. Hence, a finite union of translates of O is dense in some line segment. 
Since O is contained in the union of two lines, and the finite union of nowhere 
dense linear subsets is again nowhere dense, O intersects some line segment 
in a set that is not nowhere dense. Since the map is a piecewise transla- 
tion. On 5* is not nowhere dense in S. This contradiction finishes our proof. ^ 



4.4 The Arithmetic Graph 

The arithmetic graph gives us a good way to visualize the orbits of the 
first return map : E — )■ S. Unlike our earlier papers [SI] and [S2], the 
arithmetic graph does not play an important role in our proofs. However, we 
find it very useful to illustrate some concepts with the arithmetic graph. 
For ?/ G (0,2) we define 

Sj, = (i? X {y}) U (i? X {y - 2}) C S (39) 

The set is a union of 2 horizontal lines. (The case y = 0,2 leads to a 
trivial picture, and we ignore it.) In light of Equation[36l it makes sense to try 
to understand our orbits in terms of the triples of integers (ni, n2, ns) rather 
than in terms of the values {2niA + 2n2, n^). The parity of ni+n2 determines 
n^. Thus, the pair of integers (ni, really determines the behavior of \Ef on 
p. 

Fixing some y G (0,2), we define a map My : Z"^ ^ Yiy as follows. We 
first choose some offset value ^ E R. This offset value selects which orbits 
we focus on. Given a point (m, n) E Z'^ we define 

My^^{m, n) = {2mA + 2n + ^,y + T{m, n)). (40) 

Here r(m, n) = if m + n is odd and r(m, n) = — 1 if m + n is even. 

It is a consequence of Equation [33 that preserves the image My^^(Z^). 
That is, given (mo, no) G Z^, there is another (mi,ni) G such that 
o M(mo,no) = M(mi,ni). We write (mo, no) — )■ (mi,ni) in this case. 
Here we have set M = My^^. 

Given (y,0? form the arithmetic graph r{y,^) as follows: We include 
(m, n) as a vertex of F if and only if M{m, n) has a well defined orbit. We 
connect the vertex (mo, no) to the vertex {nii^ni) if and only if (mo, no) — >■ 
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(mi,ni). In this way, we produce a 2-valent directed graph in the plane, 
whose vertices he in . For kites, the graph F is always embedded. We 
proved this for the special orbits (i.e. those on Si) relative to any kite 
parameter in [S2]. The proof for the general orbit is similar, though we have 
not written down the details. 

Figure 4.3 shows a portion of the unique unbounded 7 component of 
the graph F(1,0~^). (The straight line segment at the bottom is just for 
reference.) The straight line is M~^(0). In [SI] we studied this component 
and showed that indeed 7 is unbounded. 




Figure 4.3: The unbounded component of F(1,0 ^). 




Figure 4.4: The unbounded component of r(0 ^,0 ^). 

Figure 4.4 shows the unbounded component of r(0"'^, 0~^). This compo- 
nent is quasi-isometric to the one in the previous picture but looks different 
locally. This is an illustration of Theorem 11.51 in action. 
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4.5 The Freezing Phenomenon 

The arithmetic graph illustrates an interesting phenomemon. Figure 4.5 
shows the portion of a graph that corresponds to the parameter y = 
Notice the long range linear order. This phenomenon becomes more and 
more extreme as ?/ — )■ 0: Longer and longer portions of the arithmetic graph 
contain these nearly linear portions. The components of the arithmetic graph 
follow along the lines, switching onto a new line at every intersection so as 
to avoid collisions. We call this the freezing phenomenon because the orbits 
seem to freeze into a characeristic shape. 




Figure 4.5: Part of 1(0"^ -20"^). 



The freezing phenomenon seems almost to be in contradiction with Theo- 
rem [LSI as we now explain. Let R be the renormalization map. The map E? 
is expanding, so the full preimage of any point y G R/2Z is dense. Hence, 
Theorem 11.51 has the following corollary. 
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Corollary 4.5 Let O be any generic unbounded orbit. Then there is a dense 
set 1/ e (0,2) such that Uy contains an orbit that is both coarsely and locally 
similar to O. In particular, this is true for a sequence {Vn} converging to 0. 

What makes Corollary 14.51 seem at odds with the freezing phenomenon is 
that the arithmetic graph is definitely changing shape as y — )■ 0, but somehow 
the global shape is always coarsely equivalent to some fixed arithmetic graph. 
The escape from the contradiction is that the bi-lipschitz constant implicit 
in the definition of coarse equivalence tends to oo as ?/ tends to 0. 

Given that both the freezing phenomenon and Corollary 14.51 are true, the 
renormalization discussed in our theorems must be expressible directly in 
terms of the multigrid system of lines. This is indeed the case. On the level 
of the multigrid, the renormalization is reminiscent both of the renormaliza- 
tion one sees for Sturmian sequences and the renormalization one sees for 
the Penrose tilings, especially when it is expressed in terms of De Bruijn's 
pentagrids. See [DeB]. So, there really is an underlying connection between 
outer billiards on the Penrose kite and the Penrose tiling. 

The freezing phenomenon works for all kites, and it is part of a larger 
phenomenon, though I don't see the renormalization scheme for a general 
kite. In [S3] 1 wrote some informal notes describing the connection betlen 
the arithmetic graph and the multigrids of lines, but I did not include any 
proofs. 1 don't currently know any. 
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5 Structural Results 



5.1 Compactification 

The map \t' : S — t- S turns out to be an infinite golden polygon exchange 
map. This quasi-periodicity is the driving idea behind our next result. Define 
the flat 3-torus 

g = T^; T = R/2Z. (41) 

Theorem 5.1 (Compactification) There is a fihered golden polyhedron ex- 
change map ^ : S — 7- S, and an injective embedding O : S — )■ S, given by the 
equation 

e(..,) = (l.i.o) + (^^.,) (42) 

which is a semi- con jug acy between \E' and \E'. 

Some terms require explanation. First, the coordinates of O are inter- 
preted as living in T. The translational part of G is somewhat arbitrary. 
We found this choice convenient for the purpose of drawing pictures. To say 
that 6 is a semi-conjugacy is to say that \E' o = O o \[' wherever all maps 
are deflned. Since is injective, we will often identify S with the subset 
0(E) C S. In this way, we think of S as a compactiflcation of S. With this 
interpretation, the semi-conjugacy condition just says that extends \E'. 

The set 0(S) is contained in a countable dense union of parallel planes 
that are transverse to the horizontal planes. We deflne 

= (i?x{y})u(i?x{y-2}). (43) 

for y e (0,2). Given that i/j preserves the set R x (y + 2Z), the map \E' 
preserves for each y. The image Q{^y) is densely contained in a single 
horizontal plane, and is a semi-conjugacy between the map : Sy — )■ 
and a exchange map on the corresponding horizontal plane. 

The semi-conjugacy between \Ef : Ei — )• Si and \E' : Ei — )• Ei is equivalent 
to the Arithmetic Graph Lemma of [SI]. Here Si is the horizontal plane 
of height 1. Similarly, Theorem 15.11 is closely related to the Master Picture 
Theorem in [S2]. See §7 for a discussion. 



40 



5.2 Structure of the Compact ificat ion 



In terms of raw data, we describe the compactification precisely in ^|TU Here 
we highlight the important features. We begin by showing pictures. Each 
picture is a different horizontal slice of S. Let R be the parallelogram with 
vertices 

(-AO); (A 2); (-A+2,0); (A+2,0); A = 0-=^ = (44) 

The solid X = R x [0,2] is a fundamental domain for the action of 2Z^ 
on and conveniently we can choose lifts of our polyhedra so that they 
give a partition of X. Figure 5.1 shows the intersection of X with the plane 
z = 0. The grey polygons represent the intersections of the polyhedra in our 
partition with this plane. 



Figure 5.1: The slice at 0. 

Below we show the slices at the heights 

0; 0-^ 0-1 20-2; 1. (45) 

The slices at heights t and 2 — t are isometric to each other, via the isometry 
which rotates 180 degrees about the midpoint of the parallelogram. Indeed, 
the map {x, y,z) — )• (2 — a;, 2 — y, 2 — z) is an automorphism of the partition. 
The slices values we choose are precisely the values in [0, 1] which contain 
some polyhedron vertices. Here are the remaining slices 
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Figure 5.6: The slice at 1. 



The shading in our pictures has the following explanation. The map \& is 
the identity on a polygon iff the polygon is colored grey. In particular, \E' is 
the identity on the shce z = [0]. 

Now we discuss the action of ^ on the partition. Say that a special affine 
involution of S is an order 2 affine map that preserves each horizontal slice 
and is an isometric rotation in each horizontal slice. The map 

T{x,y,z) ^ {-x,-3z - y,z) (46) 

is one such map. For each polyhedron P in the partition, there is a special 
affine involution Ip such that Ip{P) is also in the partition. We have 

^\p = TOip. (47) 
Lemma 5.2 ^! is a fibered golden polyhedron exchange map, as in j| . 

Proof: A direct calculation shows that every edge of every polyhedron in 
the partition satisfies Equation [331 The equations for the vertices of these 
polyhedra are listed in the appendix. 4 
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There is one other feature of the pictures we would hke to mention. The 
hnes in Figure 5.1 are present in all the pictTircs. These "persistent lines" 
divide S into 10 convex prisms - i.e. polyhedra of the form P x [0, 2], where 
P is a convex polygon in the plane. Hence, the 64 polyhedra in the partition 
can divided into 10 groups, each of which partitions one of the prisms. Figure 
5.7 illustrates this for the slice at height 1. Our Java program allows the user 
to highlight the clusters one at a time. 

In particular, the 64 polyhedra together partition the fundamental do- 
main F, and not just the torus. In other words, the outer edges in all the 
figures are actually part of the polyhedra, and not just artifacts of the way 
we have drawn the picture. A moment's reflection reveals that this issue 
arises whenever one wants to depict a partition of a torus by polygons. 



Figure 5.7: The slice at 1, divided into prisms. 

There is a second way to group the 64 polyhedra into 10 (non-vertical) 
prisms. If P is one of the original 10 prisms, then r(P) is another prism which 
is a finite union of some of the pieces. Here r is the affine involution defined 
above. Thus, the partition of S into 64 convex polyhedra is compatible with 
two interlocking partitions of S into 10 prisms each. 

The existence of the two families of prisms suggests that our polyhedron 
exchange map is actually the square of a piecewise affine map that is defined 
in terms of the two prism partitions. This is indeed the case, but we did not 
find this characterization useful. 
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5.3 The Renormalization Theorem 



We call an open set A C S atomic if A has a finite partition into golden 
polyhedra P U ... U P„ such that the first return map : A A is entirely 
defined, and a translation, when restricted to the interior of each Pk- We call 
the polyhedra Pk the atoms of A. Our definition does not uniquely define 
the atoms, but in practice our atoms will be the maximal ones. 

Suppose that A and B are two atomic sets. We say that a map f : A ^ B 
is an atomic bijection if / bijectively maps the atoms of A to the atoms of 
B, and the restriction of / to each atom is a homothety. 

We say that a layer of A is the set of all points in A whose third coordi- 
nates belong to some interval. We call f : A ^ B an atomic cover if A has 
a partition into layers Ai, A^k and B has a paritition into slabs 5i, P^ 
such that / is surjective and d-to-1 and f : Ai ^ Bj is an atomic bijection 
for each i. Here the index j depends on i. Each index j corresponds to d 
indices i. Assuming that the map / is given, we call the abovementioned 
layers the layers of A and B. 

An orbit (or orbit portion) is generic if it does not intersect any non- 
horizontal plane that is defined over Z[(j)]. 

Theorem 5.3 (Renormalization) There is a pair of atomic sets A,BgT, 
and a 3-to-l atomic covering map R : A ^ B with the following properties. 

1. R conjugates to '^\B or to '^^^\B, according as R acts as a trans- 
lation or a dilation. 

2. R acts on the horizontal planes exactly as the map R acts on R/2Z. 

3. R maps e(S) n A to 0(S) n B. 

4. Any generic orbit portion of length 812 intersects A and any generic 
orbit portion of length 109 intersects B. 

5. For any generic p E A, the orbit of R{p) E B returns to B in fewer 
steps than the orbit of p returns to A. 

6. Any generic orbit which intersects A also intersects B. 

Remark: The constant 109 is optimal. The constant 812 = 703 + 109 is an 
artifact of our proof. The smaller constant 703 would be optimal. See §10.51 
for a discussion of these matters. 
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5.4 Structure of the Renormalization 



In terms of raw data, we describe A, B, and R precisely in §14.31 and §14.41 
In this section we give the reader a feel for these objects. 

5.4.1 The Set B 

We describe B precisely in §14.31 
The partition 

0<(p-^ < 20-2 < 1 < 2 - < 2 - 20-2 < 2. (48) 

divides R/2Z into 6 intervals, Ji, J2, J31, J32, -h, -h- This is the same par- 
tition that defines the intervals Ii,...,!^, except that the value 1 has been 
inserted, so as to split I3 into two symmetric halves. The layers of B (relative 
to R) are the ones corresponding to these intervals. 




Figure 5.8: The slice of -B at ^. 



Figure 5.8 shows a slice of B at the parameter z = 0"^. The 4 dark 
polygons comprose the slice. We are showing the top half of S^. This slice is 
part of the first layer Bq. Each layer Bi decomposes into 4 convex polyhedra, 
which we call branches. Each branch further decomposes into between 9 and 
48 atoms on which the first return map is well defined. 

The fact that each layer of B decomposes into 4 branches is practically 
forced by the structure as we now explain. For z near 0, the restriction of 
to the slice is very nearly the identity. It fails to be the identity only 
in a thin neighborhood of 10 line segments. There is a set of 4 slopes such 
that each line segment has one of these 4 slopes. In this way, the 10 "active 
strips" of are nearly partitioned into 4 groups. (We say nearly because 
these strips intersect.) 
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As 2 — )■ there exist arbitrarily long orbits which remain within a single 
group. If we want to choose B so that every sufficiently long orbit intersects 
B, we must include a polyhedron that nontrivially intersects each of the 4 
groups. The choice of these polyhedra is not uniquely determined by this 
requirement, but we have made a reasonable and efficient choice. Once we 
make our choices for z near zero, the picture is determined for the entire 
interval [0, by a kind of continuation principle: The polyhedra naturally 
open up and follow along as the active strips get wider. 




Figure 5.9: The slice of B (dark grey) at — 17 + 110. 

We say that the freezing phenomenon is the tendency of the long orbits 
to align along thin strips. This phenomenon is another manifestation of 
the one discussed in §4.51 What happens in the ffist interval happens in 
each of the 6 intervals. The parameters 20~^ and 2 — 20~^ also exhibit 
the freezing phenomenon. Figure 5.9 shows the slice of B at the parameter 
—17 + 110. This parameter is fairly near 20^^. The light grey polygons are 
some periodic tiles. These periodic tiles nearly fill up the slice, but they leave 
some thin cracks. These thin cracks line up along the 4 basic directions, and 
the long orbits accumulate in the cracks. The set B fits inside the cracks, 
with one polyhedron per direction. We will list the vertices of the polyhedra 
comprising B in the appendix. The reader can see much better pictures of 
B (and A) using our applet. 
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5.4.2 The set A 



The partition associated to A consists of the 18 intervals we get by puUing 
back the J partition under the action of R. Precisely: 



• 


Within Ii, the partition is 


R- 








• 


Within I2, the partition is 


R- 








• 


Within 73, the partition is 


R- 




..,R- 




• 


Within I4, the partition is 


R- 








• 


Within I5, the partition is 


R- 




..,R- 





The 18 layers of A correspond to these 18 intervals. Again, each natural 
piece is decomposed into 4 The map R carries each branch of each each layer 
of A homothetically to the corresponding branch of the corresponding layer 
of B. As with the 5- branches, the 74-branches are further decomposed into 
atoms. 

5.4.3 The Map R 

Here we highlight the general features. Let Rijk denote the homothety that 
expands distances by (f>^ and fixes the point k). 

1. When resricted to v4(/i), the map is i^no or i?i2o- 

2. When restricted to ^(^2), the map is translation by (0, 0, 20~^). 

3. When resricted to ^(/a), the map is one of Riki for k G {0, 1, 2}. 

4. When restricted to ^(74), the map is translation by (0, 0, — 20~^). 

5. When restricted to A{I^), the map is one of i?io2 or -R112. 

Now we will show the Renormalization Theorem in action. Figure 5.10 
shows a closeup of the set A at the slice 5 — 30 together with part of an orbit. 
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Figure 5.10: Closeup of an orbit. 

Figure 5.11 zooms out, to reveal much more of the orbit in the same shce. 
The set A is at the very bottom. 
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Figure 5.11: The sUce of A at 5 — 30, together with an orbit. 
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Now for the magic trick. Figure 5.12 shows the shce of B at the parameter 
—l + (f) = R{5 — 30), together with an orbit. The orbit in Figure 5.12 looks 
different overall from the orbit in Figure 5.11, but the new orbit intersects B 
in the same way that the old orbit intersects A. In other words, were we to 
take a closeup of just the renormalization set, we couldn't tell whether it was 
the A-set inside the first shce or the B-set in the second slice. The reader 
can see many more pictures like this using our applet. 




Figure 5.12: The slice of B ai —1 + (j), together with an orbit. 

The Renormalization Theorem and the freezing phenomenon together 
combine to give a nice qualitative description of the dynamical tiling associ- 
ated to the exchange map. Consider a parameter zq such that all orbits are 
periodic. (Any parameter zq such that K^^zq) — will have this property.) 
At Zq, the slice is tiled by periodic tiles. As the parameter moves away from 
Zq, some cracks open up, and these cracks are filled with smaller periodic 
tiles associated to longer periodic orbits. If the new parameter is chosen just 
right, the shce is again tiled by periodic tiles: the original large ones and the 
much smaller ones that live in the cracks. Moving the parameter a bit, some 
new cracks open up and even smaller tiles fill in these cracks. And so on. 
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5.5 Renormalization of Orbits 



Genericity: We call a point of S generic if it does not lie in any non- 
horizontal plane that is defined over Z[(f)]. The map G carries the set of 
generic points in E into the set of generic points of S. Moreover, the transla- 
tion vectors defining all lie in Z[(f)\^. Hence \Ef preserves the set of generic 
orbits. All the atoms of A and B have their faces in the planes we have 
excluded. Hence, any generic point in A or i? lies in the interior of an atom. 
See §5.81 below for a further discussion. 

Basic Definition: The Renormalization Theorem allows us to define a 
renormalization operation on generic orbits. Let Oi be a generic orbit that 
intersects A. We know that Oi does not intersect the boundaries of any of 
the atoms. We start with pi G Oi fl A. We let p2 = R{pi) and we let O2 be 
the orbit of p2- This definition is independent of choice of pi thanks to Item 
1 of the Renormalization Theorem. For instance, if we choose the first return 
point p[ = "^^Ipi) for a suitable power of k, then R{p'i) = p'2 = '^^^{p2)- We 
write 

Oi ^ O2 (49) 

when Oi and O2 are related as above. We call O2 the renormalization of Oi. 
Given the nature of the map R, the orbit O2 is also generic. 

Distinguished Orbits: There is one situation where our construction above 
is ambiguous. We call an orbit ^-distinguished if it lies in a horizontal plane 
that contains the top or bottom of one of the A layers. These orbits are 
somewhat of a nuisance. When Oi is a distinguished orbit, we define the 
operation Oi -w O2 by including Oi in one layer of A or the other. We will 
not take the trouble to prove the O2 does not depend on the choice of layer 
because in the cases of interest to us, the uniqueness either doesn't matter 
or comes as a byproduct of our proof. 

Renormalization of Cores: We now mention a more precise kind of renor- 
malization. Say that a A- core is an orbit portion pi, ...,pn such that pi and 
Pn+i lie in A but p2, ■■■,Pn do not. We define i?-cores in the same way. If a 
is an y4-core and /3 is a S-core, we write a ^ /3 if i? maps the endpoints of 
a to the endpoints of /3. In this case, we have Oi O2, where Oi is the 
orbit containing the points of a and O2 is the orbit containing the points of /3. 
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Restriction to the Strip: So far we have been talking about the pic- 
ture in the compactification S, but we can transfer everything over to S. By 
the Compactification Theorem, we can consider E as a subset of S. With 
this interpretation, the action of \l/ is just the restriction of By Item 4 of 
the Renormalization Theorem, the operation Oi O2 preserves the set of 
generic infinite \l/-orbits. 

Coarse Equivalence: The arithmetic graph illustrates the nature of the 
renormalization map on cores. We say that a strand is the arithmetic graph 
of a core, translated so that one endpoint is the origin. The renormalization 
operation on cores gives a map from the set of all A-strands to the set of 
all S-strands. Figure 5.13 plots an A-strand in black and the corresponding 
5-strand in grey. The two strands start at the origin, which is the endpoint 
at right, and they have the other endpoint in common as well. This example 
is what we will call type-1. It corresponds to a layer of A on which R is an 
isometry. 




Figure 5.13: Corresponding A and B strands. 



Figure 5.14 shows another example. This example corresponds to a layer 
of A where i? is a dilation. In this example the origin is at bottom right, 
and we have shrunk the A-strand by a factor of Notice the remarkable 
agreement. This remarkable agreement is the subject of the Far Reduction 
Theorem below, and also the basis for Theorem II. 5 [ 
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Figure 5.14: Corresponding A and B strands. 



There is a total of 2034 = 3 x 678 A-strands and a total of 678 5-strands. 
Using our program, the tireless reader can see pictures like the ones above 
for every pair of corresponding strands. 

5.6 The Fundamental Orbit Theorem 

To state the Fundamental Orbit Theorem, we use the notation established in 
§3.41 Let T be the fundamental triangle, the subject of Theorem 11.41 Recall 
that T is the tiling of T discussed in §3. Recall that Jo and Kq are the largest 
octagon and kite in T. Say that these tiles have depth zero. Say that a tile 
of T has depth n if it is similar to Jq or by a factor of 0"'^"'. Let T{n) 
denote the finite union of tiles having depth n or less. 

Recall that ip is the outer billiards map. It is easy to check that V'"^ is 
entirely defined on the interior of U T~ . 
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Theorem 5.4 (Fundamental Orbit) Let pi G T.^ — T(2) be a generic 
point and let 

p, = r\RfM)) 

Let ai he an A-core that contains pi. Then ai 0^2, where 0L2 is the B-core 
containing p2- In particular Oi O2, where Ok is the orbit containing pk. 

The exponents (±) are written loosely in the Fundamental Orbit Theo- 
rem. The choice of (+) or (— ) in every cases is taken so that the relevant 
maps have the proper domain and range. 

The fundamental orbit says that, on the level of orbits, the maps Rfj 
and the renormalization map have the same action. The Fundamental Orbit 
Theorem is the result behind Theorem 11.41 

5.7 The Reduction Theorems 

Here we present the two results which explain the qualitative action of the 
action of R on orbits. 

We call two orbits Oi and O2 associates if one of 4 things holds. 

• Oi = O2. 

• Oi = O2, the complex conjugate orbit. 

• Oi = '4)'{02)- Here ip' is the outer billiards map. 

• Oi = i;'{02). 

Associate orbits are clearly locally and coarsely similar, equivalent. 

We write Oi — > O2 if 0[ O^, where is an associate of Ok- The two 
relations {~^) and (— t-) are practically the same. We use the latter because 
it allows us to prove our results with significantly less computation. 

Define 

S24 = [-24, 24] X [-2, 2] (50) 

Theorem 5.5 (Near Reduction) Let Oi be any generic infinite orbit that 
intersects S24. Then there is some m such that Oi — )■ ... — ?■ Om, and Om 
intersects T , the fundamental triangle. More precisely, let y G (0, 2) be any 
point such that {/?"(?/)} does not contain 0. Then there is some N with 
the following property. If Oi contains a point within 1/N of the segment 
[—24, 24] X {y}, then Oi — ?■ ... — ?■ Om and Om intersects and m < N. 
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Given an A-coie a, we define 

\a\x = min (51) 

That is, I eta; I measures how close the x-coordinates of a come to 0. We make 
the same definition for /3. 

We say that a has type-1 if the renormahzation map i? is a piecewise 
translation on the layer containing a. We say that a has type-2 if i? is a 
piecewise dilation by 0^ on the layer containing a. 

Theorem 5.6 (Far Reduction) Let a and P be A and B strands respec- 
tively. Suppose that a f3. Then the following is true. 

1. If a has type 1 then \ax\ — C < < \a\,j. + 12. 

2. If a has type 2 then (j)^^\ax\ — C < < (j)^^\a\x + 15. 
Here C is a universal constant that we don't care about. 

We think of the Near Reduction Theorem and the Far Reduction Theorem 
as geometric versions of our Descent Lemma II. The Far Reduction Theorem 
is similar to the theoretical argument we gave in order to reduce the Descent 
Lemma II to a computer calculation, and the Near Reduction Theorem is 
similar to this computer calculation. Indeed, the proof of the Far Reduction 
Theorem is almost the same as the proof of the theoretical part of the Descent 
Lemma II. 

5.8 Discussion 

Recall that A and B are partitioned into atoms such that the first return 
map of \l/ is entirely defined on the interiors of these atoms. Experimental 
evidence leads us to the following conjecture. 

Conjecture 5.7 No point on the boundary of an atom has a well defined ^ 
orbit. 

Conjecture 15.71 is subtle, because there are boundary points on which at 
least 2^^ iterates of both forwards and backwards, are well-defined. In 
particular, thousands of iterates of the first return map '^\ A are well-defined 
on such points. 
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Were the conjecture true, we could eliminate the restriction to generic 
orbits. The vertices of the atoms in the Rcnormalization Theorem all lie in 
Given the formula for the map 0, the preimage in S of these faces 
is a countable discrete set of non- horizontal lines defined over Z[(j)]. For this 
reason, © never maps a point of a generic orbit into the boundary of one of 
the atoms. This is why we work with generic orbits. 
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6 Applications 



6.1 Proof of Theorem 11.51 

We will first consider the case of Theorem 11.51 in which 1/2 = R{yi)- Here R 
is the circle renormalization map. After we finish this case, we will deal with 
the general case. Note that the winding number of an orbit coincides with 
its \E'-period. So, we will work with the \E'-period rather than the winding 
number. We will often just say period in place of \E'-period. 

Suppose that Oi and O2 respectively are generic orbits having sufficiently 
high period. We have already explained the renormalization map Oi -w O2. 
This map is 3-to-l on the level of orbits, and the inverse images of a single 
orbit lie in different horizontal planes. So, in our case, we get a bijection 
between a certain collection of orbits on Ry_^ and a certain collection of orbits 
in Ry^. These collections contain all generic infinite orbits and all generic 
periodic orbits of sufficiently high period. 

Lemma 6.1 Oi and O2 are coarsely equivalent, and the coarse equivalence 
constant is completely uniform - not even dependent on the parameters. 

Proof: We use the language of the Far Reduction Theorem. We decompose 
Oi into A-cores {ak}- At the same time we decompose O2 into 5-cores {(3k}- 
We set things up so that -w /J^ for all k E Z. Let A be either 1 of 
According to the Far Reduction Theorem, the map {x, y) — )■ (Ax, y) carries 
Oik to up to a uniformly bounded error, independent of yi and 1/2. Our 
result follows immediately from this. 4l> 

Remark: We have established the existence of a coarse equivalence constant 
that does not depend on the parameters, and this is even stronger than what 
we are claiming in Theorem II. 5[ The difference is that the parameters in 
Theorem 11.51 might be related by a long chain of renormalizations, and the 
coarse equivalence constants will probably depend on the length of the chain. 
In view of the freezing phenomenon discussed in §4.5[ there couldn't possibly 
be a uniform constant that worked for all equivalent parameters. 

Now we turn to the question of local equivalence. 
Lemma 6.2 Corresponding orbits are locally equivalent. 
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Proof: In the periodic case, each orbit is contained inside a periodic tile and 
there is nothing to prove. So, suppose that Oi -w O2 and these orbits are both 
infinite. We think of the strip S as a subset of S. With this interpretation, 
the action of \l/ on S coincides with the action of \1/ on S. So, in our proof 
we can think of Oi and O2 as orbits of Every orbit is homogeneous, so it 
suffices to produce points pj & Oj, a disk Aj containing pj, and a similarity 
from Ai to A2 which maps pi and p2 and (generically) conjugates \l/|Ai to 
$|A2. 

We can just choose pi E A and p2 E B to be points are related by R. We 
choose Ai small enough so that Ai C A, and then we let A2 = -R(Ai). By 
the Renormalization Theorem, R conjugates conjugates \l/|Ai to \I/^^|A2 on 
generic points. This is exactly what we need. 4 

Now we turn to the general case of Theorem 11.51 Let 0{y) denote the set 
of generic orbits in i?^. Say that a full subset of 0{y) is a set that contains 
all infinite orbits and all periodic orbits having sufficiently high \l/-period. 

Reffection in the x-axis conjugates ^ to and sets up an obvious and 
canonical bijection between Oiy) and 0{—y). (Here we think of —y as an 
element of R/2/Z.) So, for Theorem ll.5[ it suffices to consider the case of 
two parameters that equivalent under the group studied in §2.31 

By the Descent Lemma II, two such parameters yi and y2 are such that 

y = R-^{y,) = R--{y2). 

Let y[ = R{yi) and y'{ = R{y'i), etc. The work above gives us a bijective 
correspondence between full subsets of 0{yi) and 0{y[). Similarly, we get 
a bijection between 0{y'i) and 0{y'(). And so on. Composing all these bi- 
jections, we get a bijection between a full subset of 0{yi) and a full subset 
of 0{y). Call this the first main bijection. From the lemmas proved above, 
corresponding orbits are locally and coarsely equivalant, and the coarse equiv- 
alence constant only depends on the parameter. 

We get the same results for y2 in place of yi. Call this the second main 
bijection. Composing the first and second main bijections, we get the bi- 
jection between full subsets of Oiyi) and 0(?/2) which has all the properties 
advertised in Theorem 11.51 This completes the proof. 
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6.2 Bringing Orbits into View 

Now we explore the consequences of the Near and Far Reduction Theorems. 
We use the same equivalence relation (— t-) as in the Near Reduction Theorem. 

Lemma 6.3 Let Oi be any generic infinite \E' orbit. There is a sequence of 
orbits Oi ... On such that On intersects the rectangle [—24, 24] x [—2, 2]. 
If Oi intersects the region [— A^, A^] x [—2,2] then we can take n < AN. 

Proof: We define the type of an orbit to be the type of the A-cores com- 
prising it. Let Oi be an infinite \1/ orbit. Let O2 be such that Oi -w O2. 

We define \S\x exactly as in Equation [SH except that we use "inf" in 
place of "min" because we are dealing with an infinite set. 

Suppose first that Oi has type 2. and |Oi|a; > 24. Let a be any core and 
let /3 be such that a -w /3. Far Reduction Theorem tells us that 

|/3U<r'|a|. + 15. 

When \a\x > 24, a bit of arithmetic tells us that \/3\x < \a\x — 1- So, when 
\Oi\x > 24, we have \O2\x < \Oi\x — 1. 

Suppose now that Oi has type 2. In this case, given the nature of R, 
we have 01-^02^03, where O2 has type 2. The Compression Theorem, 
combined with the same argument as above, tells us that 

\O2\x < \Oi\x + 12; \Os\x = r'lOal. + 15. 

Therefore 

|O3U<0-'|OiU + 120-^ + 15. 

In this bit of arithmetic tells us that \O3\x < \Oi\x — 1/2 provided 

that \Oi\x > 24. In either case, at most 2 renormalizations the closest point 
of Oi closer by 1/2 units. 4 

Combining Lemma 16.31 with the Near Reduction Theorem, we obtain the 
following result. 

Theorem 6.4 Let Oi be any generic infinite orbit. Then there is some m 
such that Oi — i- ... — i- Om, and Om intersects the fundamental triangle T . 
More precisely, if Oi is an unbounded orbit that intersects [—2,2] x [— A^, A^] 
then we can take n < AN + Cy. Here Cy is a constant that depends only on 
the value of y E R/2Z such that Oi C Ry. 
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6.3 Proof of Theorem [TTH 



Theorem 11.41 deals with the fundamental triangle T and the corresponding 
tiling T. Recall from §3.4l that = T and that T~ is the image of T under 
the reflection in the vertical line x = 1. We really only care about the tiling 
T, but we find it convenient to consider both T"*" and T~ at the same time. 

Lemma 6.5 Let p be a generic point contained in the interior of a tile of 
T^. Then p has a periodic orbit. There is an upper bound on the period of 
p that depends only on the depth of the tile. Finally, the period of p tends to 
oo as the depth of the tile tends to oo. 

Proof: We use the language of the Fundamental Orbit Theorem. We check 
directly that the tiles of T(2) are periodic tiles. We handle the remaining 
tiles by induction. 

Suppose that the first statement of the lemma holds for all tiles having 
depth at most n — 1. Let Ji be a tile of depth n and let p G Ji be a generic 
point. There are indices i and j, and some choice of ±, such that Ji C T^ 
and J' = Rfj{J) is a tile of depth n — 1. Let J2 = 'il>~^{J'). Then the lemma 
holds for all generic points in J2. Let s = ip'^ ° ^fj be the similarity carrying 
J I to J2. By the Fundamental Orbit Theorem, the orbit pi G Ji is periodic if 
and only if the orbit P2 = s{j>i) E J2 is periodic. By induction p2 is a periodic 
point. Hence, so is pi. This proves the first statement of the lemma. 

For the second statement, observe that the process of renormalization 
shortens an orbit by at most a factor of 812. So, the period of generic point 
in a tile of depth n is at most 812". (This is a terrible estimate.) 

For the third statement, we recall Statement 2 of the Renormalization 
Theorem: The renormalization operation shortens the orbit length. So, if p 
lies in a tile of depth n, then p has period at least n. (This is another terrible 
estimate.) 4 

Lemma 6.6 Suppose p eT has a well-defined orbit but p does not lie in the 

interior of a tile of T . Then p has an infinite orbit. 

Proof: If p had a periodic orbit then some neighborhood U oi p would be 
such that all points in U were periodic, with the same period. But U neces- 
sarily contains infinitely many tiles of T. This contradicts Lemma 16.51 ^ 
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Lemma 6.7 No point on the boundary of tile ofT has a well-defined orbit. 

Proof: Lemma [6.61 shows that p cannot have a periodic orbit. On the other 
hand, suppose that p e dJ has an infinite orbit, where J is a tile of T. Then 
for any n, there is an open neighborhood Un such that p G [/„ and the first 
n iterates of \l/ are defined on p. These iterates all act by translation, so 
no point in Un has period less than n. But f/„ contains some generic points 
of J, no matter how large n. This contradicts the second statement of the 
previous lemma. ^ 

Corollary 6.8 Every tile of T is a finite union of orbit tiles. 

Proof: Let J be a tile of T. The same argument as in the first half of 
Lemma 16.61 shows that there is a uniform bound on the period of p. Let 
N be this bound. Since \Ef is a polygon exchange map, there is a finite set 
of lines in the plane such that the first N iterates of \l/ are defined in the 
complement of these lines. These lines partition J into finitely many smaller 
convex polygons, and is well defined and periodic on the complement of 
each of these lines. This shows that J is covered by a finite union of periodic 
tiles. Since no point of dJ has a well-defined orbit, none of these periodic 
tiles crosses the boundary of J. 4 

Theorem 11.41 follows from the results above. We prove one more related 
result in this section. 

Lemma 6.9 Let L be any horizontal line segment such that LflT^ contains 
contains no generic points with infinite orbits. Then there is a uniform upper 
bound to the period of any orbit on L H . 

Proof: When the height of L lies in the interval Ji, L fl T"^ is contained in 
the closure of the union of the two largest octagonal tiles. The result is cer- 
tainly true in this case. In general, we either have R"'{y) = [0] or R^{y) G /i 
for some n. In either case, it follows from the symmetry of T that L n 
is contained in the closure of the union of a finite number of orbit tiles. But 
then our bound follows from Lemma 16. 5[ 4|k 
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6.4 Unboundedness 



Theorems 11.51 and 16.41 are the workhorses in our overall proof. They allow 
us to transfer statements about orbits that intersect T, the fundamental 
triangle, to statements about orbits in general. The results in this section 
illustrate the technique. 

We say that a fundamental orbit is an infinite orbit that intersects the 
fundamental triangle T. 

Lemma 6.10 Every fundamental orbit is unbounded in both directions. 

Proof: Every fundamental orbit is contained in the union of two lines, and 
every line intersects S, the fundamental fractal, in a nowhere dense set. Now 
we apply our unboundedness criterion. Lemma 14.41 ^ 

Now we promote Lemma 16.101 to a statement about all orbits. 
Lemma 6.11 Every infinite orbit is unbounded in both directions. 

Suppose first that O is a generic infinite orbit. By Theorem 16. 4[ we can 
write O = Oi 0„, where On is an infinite orbit that intersects T and 

is coarsely equivalent to Oi. But we have already seen that 0„ is unbounded 
in both directions. Hence Oi is also unbounded in both directions. 

Now suppose that O is an infinite (non-generic) orbit that is bounded in, 
say, the forward direction. Then, by compactness, some point of O lies in the 
accumulation set of O. But then every point of O lies in the accumulation set 
of O. This implies that the set of accumulation points of O is uncountable. 
Indeed, this accumulation set contains a Cantor set. 

The set of non-generic points on any horizontal line is countable, and 
O lies on 2 horizontal lines. (We take the \E'-orbit, as usual.) Hence, there 
is some generic accumulation point p of O. Note that p necessarily has an 
infinite orbit. Hence p has an orbit that is unbounded in both directions. 

In particular, for every there is some e such that any point within e of 
p follows the orbit of p for steps in either direction. Hence, p cannot be 
the accumulation point of an orbit that is bounded in one direction or the 
other. This contradiction shows that O is unbounded in both directions. ^ 
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6.5 Ruling out Some Heights 

In this section we take care of a few annoying technical details that will make 
our arguments in the next section go more smoothly. 

Lemma 6.12 Let y = (p. Then Ry has no generic unbounded orbits. 

Proof: By the Descent Lemma I (or a direct calculation), the point is 
2-periodic with respect to R. We have R{(p) = 2 — 0. Let Li denote the 
horizontal line y = (p. Let L2 denote the horizontal line y = 2 — (p. Let T be 
the fundamental triangle. Note that Li fl T is just a single point, and this 
point is non-generic. At the same time, L2 fl T is contained entirely in the 
periodic tile Kq. Hence neither line intersects T in any point that could have 
a generic unbounded orbit. 

Suppose, on the other hand, that Ry has a generic unbounded orbit Oi. 
By Theorem 16 . 41 there is some m such that Oi — )■ ... — )■ Om and Om intersects 
T. The intersection Om fl T is a generic unbounded orbit by Theorem 11.51 
The operation of replacing an orbit by an associate simply switches from Ry 
to R^-y Hence, Om fl T is a point of either Li fl T or L2 fl T. This is a 
contradiction. ^ 

Corollary 6.13 Suppose that y E C — C* . Then Ry has no generic un- 
bounded orbits. 

Proof: We have y E Z[(f)] by the construction of C. By Lemma [2. 7[ we see 
that y = m + n(j), with m even. But then ?/ ~ or y ~ 0. Here ~ denotes 
G2-equivalence. The map \E' is the identity on R^. So, Theorem 11.51 rules 
out the possibility that ?/ ~ 0. On the other hand. Theorem 11.51 and the 
preceding lemma rule out the possibility that y ~ 0. 

Corollary 6.14 Any generic unbounded orbit that intersects the fundamen- 
tal triangle T contains a point of the form {x, y) where y G C*. 

Proof: By Theorem II. 4[ and the structure of the fundamental fractal S, 
any such point (x, y) must have y E C. But our previous result rules out the 
possibility that y E C — C*. 
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6.6 Proof of Theorem 11.61 



Now we turn to the question of self-accumulation. Lemma 14.31 tells us, in 
particular, that there is a canonical bijection between the orbits in S+ 
and the orbits in S_. We call two such orbits partners. Just for this 
section, we introduce the notation Oi =^ O2 to mean that Oi -w O2, where 
O2 is the partner of 02- The following result is just a reformulation of the 
Fundamental Orbit Theorem. 

Corollary 6.15 Letpi E — T(2) he a generic point and letp2 = Rfj{pi). 
Let Oi and O2 respectively be the orbits of pi and p2. Then Oi ^ 02- 

Proof: If O2 is an orbit in S_, and ip{0'2) intersects S+, then the partner 
of O2 is the orbit of any point of i/j{02). The corollary follows from this fact, 
and from the Fundamental Orbit Theorem. ^ 

Let Uy denote the set of unbounded orbits in Ry. We also have the 
following lemma, which is just a reformulation of part of Theorem II. 5 [ 

Lemma 6.16 (Bijection Principle) Let Oi and O'l be two generic orbits 
in Uy n S+ such that Oi =^ O2 and 0[ =^ O2. Then Oi = 0[. 

Proof: We have 0[ ^ O2 and Oi O2 where O2 partner of O2. On 

the level of infinite orbits, the renormalization map is 3-to-l. At the same 
time, the 3 preimages of an orbit lie in different sets Uy. Since Oi and O'l lie 
in the same set Uy, we have Oi = 0[. 4|k 

Lemma 6.17 Let p E T be a generic fundamental orbit. Then the orbit of 
p is self- accumulating. 

Proof: We already know that the orbit of p is unbounded in both directions. 
By Lemma |6.14[ we have p = {x,y), where y G C*. Let 0{p) denote the 
orbit of p. We refer to the notation and terminology in §3.51 Let T{p) denote 
the remormalization set for p. Let q G T(p). Here p and q lie on the same 
horizontal line, and 

p = Po-^Pi... ^Pn, q = qo ^ qi. . . ^ qn, Pn = Qn- (52) 
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By the Corollary above, 



0{pq) ... 0{pn); 



0(g„_i) 0{qn) = 0{pn). 



(53) 



Since Pn-i and lie on the same horizontal line, the Bijection Principle 
says that = 0(p„_i). It now follows from induction on n that 



By the Density Lemma, T(p) is dense in A = SCiL, where L is the horizontal 
line containing p. By the Horizontal Lemma, A is a Cantor set. Hence p is 
an accumulation point of 0{p). But, 0{p) is homogeneous. Hence, every 
point of 0{p) is an accumulation point of 0{p). This proves that 0{p) is 
self-accumulating. ^ 

Now we prove Theorem ll.6[ Let Oi be a generic infinite orbit. By The- 
orem [^31 we have Oi — )■ ... On where 0„ is an infinite orbit that in- 
tersects the fundamental triangle. But we have already shown that On is 
self-accumulating. By Theorem II. 5[ the orbits Oi and 0„ are locally similar. 
Hence Oi is also self-accumulating. 



We've already proved that every orbit is either periodic or unbounded in 
both directions. It only remains to show that the union U of the unbounded 
orbits has Hausdorff dimension 1. By Lemma 13.11 the set of generic points 
in S has Hausdorff dimension 1. All such points have well-defined orbits 
and, by Theorem 11.41 all the orbits are unbounded. Hence U contains a 1 
dimension set. Hence dim(?7) > 1. 

Let G be the set of generic points in R^. note that — G has Hausdorff 
dimension 1 because it is a countable set of lines. Combining Theorems 11.51 
and 16.41 we see that every generic point p E U has a neighborhood A such 
that A n f/ is similar to a subset of S. Hence dim(f/ fl G fl A) < 1. But 
we can cover any compact subset of U by finitely many such neighborhoods. 
Hence dim(t/ fl G) < 1. Since U — G is contained in a countable family of 
lines, we have dim([/ — G) < 1. Hence dim([/) < 1. 




T(p) c 0{p). 



(54) 



6.7 Proof of Theorem 11.11 
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6.8 Proof of Theorem [Q] 



First we establish part of Theorem 11.91 Let U* denote the set of generic 
unbounded orbits in Ry. 



Lemma 6.18 Uy is empty if U* is empty. 



Proof: Let R be the renormahzation map. Suppose R^ly) = for some 
n. All orbits on the line y = have period 1, and renormahzation decreases 
periods by at most a factor of 812. Hence, there is a uniform bound to the 
period of any generic point in Uy. 

On the other hand, suppose that y is not in the inverse image of 0. Let 
Pi be any generic point that is, say, within 1 unit of y. Let Oi be the 
orbit containing pi. Combining Lemma [6.31 and the second statement of the 
Near Reduction Theorem, we see that there is some uniform m such that 
Oi — 7- ... — 7- Om and Om intersects T. But the horizontal line containing 
Om n T has no generic points with infinite orbits. Hence, by Lemma 16. 9[ 
there is a uniform upper bound to the period of Om- Hence, there is a 
uniform upper bound to the period of Oi. 

Finally, we can take a sequence of generic points approximating a sup- 
posed infinite orbit. This is incompatible with the uniform upper bound we 
have on the periods of this approximating sequence. ^ 

Now we prove Theorem 11.21 Let Y C R/2Z denote those y such that 
Ry contains unbounded orbits. Let y E C*. The horizontal line of height y 
intersects 5* in an uncountable set which must have points with well-defined 
orbits. By Theorem II. 4 [ these points are not periodic. By Lemma [6.101 these 
points have unbounded orbits. Hence C Y. Theorem II. 51 now shows that 
Y contains all points y such that y ~ c and c G C*. 

Now for the converse. Suppose that yi G Y. By Lemma I6.18| we can 
assume that Ry^ has a generic unbounded orbit Oi. By Theorem 16.41 we 
have Oi — )■ ... — )■ Om with Om intersecting T. Let yk be such that Ok is an 
orbit of Ry^. Passing to associates preserves the G2-equivalence class, and 
so does R. Hence yi ^ ym- 



By Theorem 11.51 the orbit Om is generic and unbounded. By Theorem 
11.41 we have ym G C. By Lemma I6.14| we have ym ^ C — C*. Hence 
ym £ C"* and yi ~ ym, as desired. Finally, we compute easily that dim(C) = 



log(3)/log((/.3). 
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6.9 Proof of Theorem [TSl 

As we remarked above, the concepts of winding number and \l/-period coin- 
cide. We will work with \E'-period, as usual. 

Lemma 6.19 Every horizontal line in S contains a dense set of periodic 
orbits. 

Proof: Let U denote the union of unbounded orbits. If L fl ?7 is empty, 
there is nothing to prove. Otherwise, by Theorem 16.41 and Theorem II. 5[ the 
set L n U is locally similar to L' fl U, where L' is a line segment in T, the 
fundamental triangle. But f/ fl T is nowhere dense in each horizontal line. 
Hence U (1 L' is nowhere dense in L'. Hence U (1 L is nowhere dense in L. d|b 

Lemma 6.20 Let y G R/2Z be any value other than 

0; 4-20; -2 + 20; 2 

Then Ry contains a periodic orbit that intersects B, the renormalization set. 

Proof: Let y G R/2Z be any value other than the ones listed. Let H be 
the horizontal plane of height y. The plane H intersects B in an open set. 
Let L C S be the horizontal line of height y. Since Q{L) is dense in H, there 
is an open subset V C L such that 9(V^) C B. By Lemma I6.19[ there is 
a dense subset of V consisting of periodic points. Any periodic point in V 
works for us. ^ 

Lemma 6.21 Let y G R/2Z. Suppose it never happens that R"'{y) = [0]. 
Then Ry contains periodic points having arbitrarily high -period. 

Proof: Our proof refers to the terminology used in the Renormalization 
Theorem. Let yn = R^{y). For our analysis, we identify E with a dense 
subset of S. Define = i?^^ fl S. Let be the horizontal plane of height 

Vn- 

Referring to Lemma I6.20[ the heights of the excluded horizontal planes 
all lie in the inverse image of 0. For this reason, S„ contains a periodic tile 
P whose orbit intersects B in an atom that is a translate of P. 
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Since S„ is dense in S„, the periodic tile P must intersect S„ in an open 
interval. In particular E„ contains a generic periodic point whose orbit 
intersects B. By the Renormalization Theorem, we can find a generic orbit 
0'^_i in such that O'^^^ 0„. The period of 0'^_i is longer than the 
period of 0„. 

0,'j_^ necessarily intersects A. By the Renormalization Theorem, 0^_^ 
also intersects B. So, we can set 0„_i = 0^„^ and repeat the above ar- 
gument to produce a generic periodic orbit 0'^_2 in E„_2 that intersects B. 
And so on. The final orbit Oi has period at least n — 1. But n is arbitrary. ^ 

The full inverse image of [0] is precisely 2Z[(f)\. Hence, if y 2Z'[0], then 
Ry contains orbits of arbitrarily high \l/-period. 

Lemma 6.22 Suppose thaty = m + ncj) where bothm andn are even. Then 
there is a uniform hound on the ^-period of any point of Ry with a well 
defined orbit. 

Proof: We consider generic points first. Suppose Ry contains a generic orbit 
with an enormous winding number - either finite or infinite. There is some 
n such that R"'{y) = [0]. If the winding number of our orbit is too large, 
then we can renormalize this orbit more than n times. But this contradicts 
the fact that \1/ is the identity on Sq, and no such orbit on this plane has a 
renormalization. 

So, every generic point of Ry is periodic and we have a uniform bound 
on the periods. But any point in Ry can be approximated by a sequence of 
generic points. The uniform bound on the approximating sequence immedi- 
ately gives the same uniform bound on the limit, provided that the limit has 
a well-defined orbit, ^(t 

Let y G 2Z[0] be some point for which the following property (*) holds: 

R\y)^[l]; k = 0,...,n. (55) 

Then the same argument as in Lemma 16.211 shows that Ry contains generic 
points having \l/-period at least n — 1. Finally, for any fixed n, there are 
only finitely many points of 2Z[(j)] which fail to satisfy Equation 153 This 
establishes Theorem 11.31 
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6.10 Proof of Theorem [TTO] 

This is one of the more subtle results. We will try to break the proof down 
into small steps. 

Lemma 6.23 Let y be one of the 4 values 

-1 + -3 + 30, 3-0, 5-30. 

Let A be the intersection of the line R x {y} with the fundamental triangle 
T. Then there is a uniform upper bound to the period of any point that is 
sufficiently close to A. 

Proof: Let Ai, A4 be the 4 line segments in question. We check by direct 
inspection that there is some eo > with the following property. The eo 
neighborhood of A^ is contained in a finite union of tiles of T. Indeed, the 
two lines A2 and A4 are disjoint from T, the interior of Ai fl T is contained 
in the interior of a single periodic tile, and the interior of A3 fl T is contained 
in the interior of the closure of two periodic tiles. 4 

Let y = m + ncj) where m and n are odd. Let L = R x {y}. Let N = Ny 
be as in Lemma 12.51 

Lemma 6.24 Let k > N be a fixed integer. There is some e > with the 
following property. Suppose that Oi is an orbit that contains a point within 
e of L. Suppose Oi — )■ ... — >• Ok and that Ok intersects the fundamental 
triangle. Then there is a uniform bound on the period of Oi. 

Proof: By Lemma 12. 5^ and the fact that i?*^ expands distances by at most 
(j)^^, the orbit Ok must intersect the fundamental triangle in a point very 
close to one of the 4 line segments listed in Lemma [6.23[ provided that e is 
small enough. The point here is that, when we keep track of which set Ry 
contains our successive orbits, the renormalization operation implements the 
map R and switching from an orbit to an associate implements the map S. 

If Ok is very close to one of the segments in Lemma I6.23[ then there is a 
uniform upper bound to the period of Ok- But the renormalization operation, 
and the operation of switching to an associate, only can decrease the period 
by a uniformly bounded factor. Hence, there is a uniform upper bound on 
the period of Oi. ^ 
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Lemma 6.25 Let p E L he some point. Let M be any positive integer. 
Then there is some e > with the following property. If Oi is a generic 
periodic orbit that comes within e of p, then there is some k > M such that 
Oi — )■ ... — )■ Ok and Ok intersects the fundamental triangle. The same k 
works for all orbits satisfying the hypotheses. 

Proof: Combining Lemma [673] and the second statement of the Near Reduc- 
tion Theorem, we see that there is some k, which works for all choices of Oi, 
such that Oi — )■ ...O^ and Ok intersects the fundamental triangle. The only 
trouble is that we might have k < M. However, we can renormalize again 
and then, if necessary, apply Theorem 16.41 again. This produces a larger value 
of k. We keep going like this until we arrive at some k > M, and then we 
stop. ^ 

Corollary 6.26 Letp E L be any point. Then there is some e > and some 
constant Z with the following property. If Oi is a generic periodic orbit that 
comes within t of p then Oi has period at most Z . 

Proof: Choose M = N, the constant in Lemma I6.24[ If e is small enough 
then we have some fixed k > M such that Oi Ojt and Ok intersects 

the fundamental triangle. This fixed k works for any orbit that comes within 
e of p. But now Lemma [6.231 applies to Oi. 4 

Suppose that the conclusion of Theorem IL 101 is false. Then we can find 
some p E L which is the accumulation point of points having unbounded 
orbits. By Theorem II. 1[ we also have a sequence of generic periodic points 
converging to a point of L whose period tends to oo. This contradicts the 
corollary we have just proved. This establishes This proves Theorem I l.lOi 

We mention the following corollary of Theorem 11.101 

Corollary 6.27 Let y = m + ncf) with m and n odd. Let K be and compact 
subset of . Then there is some e > 0, depending on K,n,m, such that no 
point of K within e of the line Ly has an unbounded orbit. 

Proof: This is just an application of compactness. 4 
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6.11 Proof of Theorem [TTSl 

Let Yr be as in Theorem 1 1.81 As above, let Uy denote the union of unbounded 
orbits in Ry. 

Lemma 6.28 Yr is nowhere dense. 

Proof: Rather than work precisely with the set mentioned in Theorem ll.8[ 
we define to be the set y ^ Y such that Uy contains a generic orbit that 
intersects the rectangle [— r, r] x [—2,2]. We have Y^-i C F/ C Yr+i, so it 
suffices to prove Theorem 11.81 for Y^ in place of Yr. 
The analysis in Lemma 16.31 shows that 

i?"(y;') U C 1^24; Vn > 4r. (56) 

Suppose that Y^ is not nowhere dense for some r. Since is an expand- 
ing map, and Y/. is not nowhere dense, the left hand side of Equation |56] is 
dense in R/2Z for sufficiently large n. Hence 1^24 dense in R/2Z. But 
this contradicts Corollary 16.271 ^ 

We have already computed that dim(C*) = log(3)/ log(0'^). Observe that 
every neighborhood of the Penrose kite vertex 0) contains a similar copy 
of S. Hence dim(yr) = dim(C*). 

6.12 Proof of Theorem [TTOl 

When we proved Theorem II. 2[ we proved the first statement of Theorem 11.91 
Now we prove the second. 

The set of non-generic points is contained in a countable union of lines. In 
particular, the set of non-generic unbounded orbits is contained in a countable 
union of sets of the form 

Ln{RxY), (57) 

where Y is the set of y such that Ry has unbounded orbits. By Theorem 
11.21 we have dim(y) = log(3)/ log((/)^). But then the union of non-generic 
unbounded orbits is contained in a countable union of sets that have the 
same dimension as Y. This completes the proof. 
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6.13 Proof of Theorem [T7fl 



Almost every statement in Theorem II .71 is an immediate consequence of The- 
orems ll.4[ ll.5[ and 16.41 The one statement that is not immediate is that a 
neighborhood of the dynamical tiling about a generic point with unbounded 
orbit is isometric to a small patch of T. 

The results mentioned above only imply that the abovementioned neigh- 
borhood is similar to a patch of T. However, tracing through our argument, 
we see that the similarity factor is 0^^^ for some integer k. But T is a self- 
similar set with expansion constant 0^. So, we can take /c = in the similarity 
factor between the neighborhood of interest to us and a suitable patch of T. 

6.14 Proof of Theorem [TTTI 

Let y = m + ncf). We already know that 1 G C*. So, by Theorem ll.2[ we 
know that has unbounded orbits provided that m is odd and n is even. 

Conversely, suppose that has unbounded orbits. The set C — C* 
contains points of the form m + n0 with m even and n having either parity. 
Combining Theorem 11.21 and Corollary I6.14[ we see that m cannot be even. 
When m is odd and n is odd. Theorem 11.101 shows that has no unbounded 
orbits in this case. The only case left is when m is odd and n is even, as 
claimed. 
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7 Computational Methods 



7.1 Golden Arithmetic 

A large percentage of our proofs involve exact computer calculations, done 
over the number ring We abbreviate these kinds of computer calcula- 

tions as golden arithmetic. Here we describe how the computer does golden 
arithmetic, for the sake of making our calculations completely reproducible. 
There exist computer packages, such as Pari, which perform these kinds of 
calculations in extreme generality. A software package like Pari is vastly 
more complicated than our own much more limited collection of routines. 

We represent the number ao + ai0 as the integer pair (oq, ai). The com- 
puter can only represent finitely many such numbers, but the large finite set 
of representable numbers is suitable for our purposes. 

The Ring Operations: We have the following obvious rules 

(oo, oi) ± {bo, bi) = {ao ± 60, oi ± bi); 

{ao, ai) X (60, bi) = {aobo + aibi, aobi + aibo + aibi), (58) 

which represent the ordinary ring operations in Z[(f)]. As long as all integers 
stay less than, say, 10^ in absolute value, the computer adds and multiplies 
them correctly. The Galois map r(a -|- 60) = a — 6/0 is a ring automorphism 
of Z[(l)\. In terms of our representation, we have 

r: {a,b) ^ {a + b,-b). (59) 



Division: Since Z[(f)\ is not a field, we cannot generally perform division 
in Z[(j)\. However, it occasionally happens that we know in advance that 
•y = a/f3 lies in Z[0] and we want to find 7 given a and (3. In this cases, we 
find 7 by computing 

a X t(B) , , 

The number in the denominator is an integer, and we find the quotient on the 
right hand side by dividing the coefficients of the numerator by this integer. 
Before returning the value, we verify that the resulting element belongs to 
Z\^ and satisfies the equation a = /3 x 7. 
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Positivity: Now we explain how we check that ao + ai(j) is positive. Our 
method, apphed to + returns a value of true only when ao + ai0 > 0. 
We call ao + ai0 strongly positive if 

ao/ioo + cti/ioi > 0; cto/ioi + aifio2 > 0. (61) 

Here /„ is the nth Fibonacci number. Since the successive quotients of Fi- 
bonacci alternately over and under approximate 0, the fact that both linear 
combinations are positive guarantees that ao + ai0 is also positive. In short, 
a strongly positive element of Z[(f)] is positive. Certainly there are positive 
numbers that are not strongly positive. However, we do not encounter these 
numbers in our calculations. 

In Equation [611 we use the Biglnteger class in Java, which does exact 
integer arithmetic for integers up to many thousands of digits long. We 
might have used the Biglnteger class for all our calculations, but this would 
make the calculations much slower and the computer routines much more 
tedious to program. Thus, we use the Biglnteger class only when needed. 

Equation [61] is the basis for all our computations that involve inequalities 
between elements of Z[(j)]. For instance, to verify that (ao,ai) represents a 
number larger than (6o, &i), we apply the test to (co, Ci), where Cj = aj — bj. 

Golden Structures: The GoldenReal, a pair (ao,ai) as above, is our basic 
object. We also define more complicated objects based on the GoldenReal: 

• A GoldenComplex is a pair x+iy, where x and y are both GoldenReals. 

• A GoldenVector is a tuple of GoldenReals. 

• A GoldenPolytope is a finite list of GoldenVectors. 

The usual operations for these objects are done using the ring operations 
described above. 

Remark: It is worth mentioning that our graphical user interface mainly 
operates with floating point arithmetic, for the purposes of speed. We mainly 
use the special arithmetic when we need to do rigorous calculations for the 
purpose of making a proof. The special structures are all isolated in separate 
files, so as not to interfere with the rest of the program. 
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7.2 Overflow Error 



There is one computational issue that we must face when we do exact in- 
teger arithmetic calculations. The computer does not reliably perform the 
arithmetic operations on very large integers. For example, when multiplying 
together ordinary integers, or ints, our computer tells us that 

100000 X 100000 = 1410065408. 

Each int is a length 32 binary string. One of the bits records the sign of 
the integer and the remaining bits give the binary expansion. The problem 
with the solution to 100000 x 100000 is that it requires more than 32 bits to 
express the answer. 

To give us a bit more flexibility, we use longs. A long is a 64-bit repre- 
sentation of an integer. The calculation above comes out right when we use 
longs in place of ints. 

Here are some conservative bounds on what the computer can do with 
longs 

• The computer reliably computes the value ai ± a2 provided that we 
have max(|ai|, |a2|) < 2^°. 

• The computer reliably computes the value Oi x a2 provided that we 
have max(|ai|, |a2|) < 2^°. 

Whenever we perform one of the basic operations above, we check that the 
conditions above hold. The code is set up to interrupt the calculation if one 
of the conditions above is not met. 

The only other operations we perform are Galois conjugation within Z[0] 
and the guided form of division mentioned above. For conjugation, wc make 
sure that the coefficients of the GoldenReal are less than 2^° in absolute 
value. 

The division operation requires a special explanation. When we perform 
the division operation to find the quotient c = a/6, we do not bound the sizes 
of the coefficients. Rather, when we then verify that indeed a = 6 x c, we 
check the bounds of the coefficients of b and c. This means that the equality 
a = b X c has really been verified. 
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7.3 Some Basic Operations and Tests 



Here we explain some basic operations we perform on polygons and polyhe- 
dra. When we have golden polygons and golden polyhedra, the operations 
are all done using the golden arithmetic described above. 

Positive Convexity: We say that a polygon P is positively convex if P 
is convex, and if the orientation of dP given by the ordering on the vertices 
is counterclockwise. P is positively convex if and only if 

Im(^2i2;3i) > 0; Zij = Zi- Zj. (62) 

for every triple Zi,Z2, z^ of consecutive vertices of P. Note that a positively 
convex polygon need not be strictly convex. 

Containment Test for Polygons: Let P be a positively convex poly- 
gon with vertices Zi,...,Zk- A point w lies in the interior of P if and only 
if 

lm(^{zi - w){zi+i - w)^ > 0; (63) 
holds for every index i, with indices taken mod k. 

Strict Convexity for Polyhedra: To check that a golden polyhedron is 
strictly convex, we exhibit a golden vector V — V{P, v) such that v-P > v'-V 
for all other vertices v' of P. In practice, we search for V amongst all golden 
vectors whose coefficients have real and imaginary parts of the form a + b(j) 
with max(|a|, \b\) < 2. In general, one would need to make a more extensive 
search or else have some special information about P. 

Golden Interior Points: The most natural interior point of a polytope 
is its center of mass. However, the center of mass of a golden polytope might 
not be a golden vector. Here we describe a less canonical way of picking an 
interior point of a golden polytope which results in a golden vector. 

Given two golden vectors Vi and V2, the vector a(Vi, V2) = (1>~^V2 + (I)~^V2 
is a golden vector that lies on the line segment V1V2. Note generally, suppose 
we have a golden polytope with vertices Vi, ...,Vn. We define 

W2^a{V,,V2y, Ws^a{W,,Vs); ... Wn^a{W,,Vn). (64) 

then Wn is a golden vector in the interior of the polytope. 
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Disjoint Interiors I: Here is one method we use to verify that two golden 
polyhedra P and P' have disjoint interiors. We exhibit a golden vector W 
such that V ■ W < V ■ W for all pairs (V, V'), where ^ is a vertex of P and 

V is a vertex of P'. In practice, we search for W amongst all vectors E x E', 
where E is an edge of P and E' is an edge of P'. 

Disjoint Interiors II A second way we will know (without directly testing) 
that two distinct polygons or polyhedra Pi and P2 have disjoint interiors is 
that there is some map f : Pk ^ such that / is entirely defined and 
affine on the interior of Pk, but / does not extend to be defined and affine 
in a neighborhood of any point of dPk- If Pi and P2 did not have disjoint 
interiors, then some point of dPi lies in the interior of Pi, and then the map 
/ is then defined and locally affine in a neighborhood of this boundary point. 
This is a contradiction. 

Coplanarity Test: Let Vi,...,Vk be finite list of golden vectors. Consider 
the successive vectors 

Ni = {V,+i-V^)x{V,+2-V^), (65) 

with indices taken mod k. The vectors are coplanar if and only if iVj and 
Ni^i are parallel for all i. 

Raw Face Enumeration: Here is how we enumerate the faces of a poly- 
hedron that is given in terms of its vertices. Let V{P) be the vertex set of 
a polyhedron P. For each subset S C V{P), having at least 3 elements, we 
first check if is a coplanar set of vertices, as explained above. If 5* is not 
coplanar, we eliminate 5* from consideration. 

Assuming that S is coplanar, let Vi,...,\4 be the vertices of S. Let Ni 
be the vector from Equation [651 Let V = Vi + ...Vk- Also, let kS denote the 
result of scaling all vectors of by a factor of K. Likewise define kP. We 
try to check that either V ■ Ni < W ■ Ni for all vertices W of kP or that 

V ■ Ni>W ■ Ni for all vertices W of P. If one of these two things is true, 
then the center of mass of the convex hull of S is contained in dP, and this 
means that 5* is a subset of a face of P. The reason why we scale everything 
by k is that we want to work entirely in Z[(f)\. (To avoid scaling, we could 
have used the golden interior point described above in place of V , but we 
didn't.) 
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Edge Enumeration: Given tlie list of vertices and the raw list of faces, 
we find the edges as follows. A pair of distinct vertices (fi,f2) is an edge of 
P if any only if P has 4 distinct vertices ^1,^2, 1^1,^2 such that (^1,^2,^1) 
and {vi,V2,W2) are two distinct faces of P. For each pair (t>i,t>2), we try to 
find the pair (^1,^2) with the property just mentioned. If we succeed, we 
add (fi,f2) to the hst of edges. 

Polished Face List: There is one more step, which we perform for the 
sake of tidiness, and also for the sake of making some of our other compu- 
tations go more quickly. We would like the vertices of each face to be listed 
so that they go in cyclic order around the face. Our convention is that the 
vectors produced by Equation ESI should be outward normals. To order the 
vertices around a face, we first use the raw face list to enumerate the edges 
of the face. We then use the edges, in a fairly obvious way, as a guide for 
placing the vertices. 

Slicing: As an application of our edge and face enumeration, we explain 
how we compute the intersection of a golden polyhedron P and a horizontal 
plane 11 that contains no vertices of P. Every edge of P will satisfy Equation 
[33] and 11 will have height in Z[(j)], so P H 11 will be a golden polygon. 

Every edge of P either crosses 11 at an interior point or else is disjoint 
from n. We enumerate those edges ei, ...,6^ of P which cross 11. We order 
the edges so that and Cj+i lie in a common face. Finally, we set Vi = nncj. 
Then vi, ...,Vk are the vertices of the polygon P fl 11, and they are cyclically 
ordered around the boundary of P fl 11. 

Containment Tests for Polyhedra: Let F be a point and let P be a 
polyhedron. Let Pi, P^ be the faces of P. Let Wi, be such that Wj 
is a vertex of Fj for each j. Finally, let Nj be a vector that is perpendicular 
to Fj. We compute Nj by applying Equation |65] to the first 3 vertices of 
Fj. The point V lies in the interior of P if and only if V ■ Nj < V ■ Wj for 
j = 1, ...,k. Similarly, V lies in dP if and only if V ■ Nj < V ■ Wj for all j, 
and there is equality for at least one index j. 

Remark: The tests we have described above look fairly intensive, but for 
the most part we will be applying them over and over to the same list of 64 
polyhedra - the ones that define our polyhedron exchange map. For these, 
we pre-compute everything and store the results in a look-up table. 
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7.4 Covering Tests 



Here we describe how we show that a given polygon or polyhedron is covered 
by a finite union of other polygons or polyhedra. We will explain the polygon 
case first. We then reduce the polyhedron case to the polygon case by a trick 
involving slicing. 

Let Q be a golden polygon and let Let Pi, P„ be a finite list of golden 
polygons. Here we explain the computational test we perform in order to 
prove that Q c[jPk- 

Eliminating Spanning Edges: Say that Pk has a spanning edge relative 
to Q if some edge of Pk intersects Q in an interior point, but has no vertices 
in the interior oi Q. In practice, we eliminate all spanning edges by adding 
an extra "vertex" of Pk to each edge. Since we want our polygons to remain 
golden, wc pick a golden interior point for each edge. See above. 

Our method is not guaranteed to eliminate all spanning edges, but in 
practice it does. In all cases but one, we also have P^ d Q for all k. In this 
case, the only spanning edges are those that have both vertices on dQ. For 
such spanning edges, the addition of a single vertex does the job. In the one 
remaining case, namely the one involving the rectangle Q — T,* in the next 
chapter, there are no spanning edges at all. 

The Matching Property: Let v be an interior vertex of some Pj. Let 
Pi, Pjfc be the (re- indexed) union of all the polygons P, such that v e dPi. 
When V is a vertex of Pi, we let be the edge of Pj such that v is the 
leading vertex, and we let e~ be the edge of Pj such that v is the trailing 
vertex. When p lies in the interior of an edge Cj of Pj, we let and be the 
two oriented copies of e^, labelled so that the indices of the endpoints of ef 
are cyhcally increasing. We think of these 2k edges as vectors, all oriented 
away from v. For each i G {1, ...,k} we check computationally that there is 
exactly one j G {1, ...,k} such that the two vectors e~ and are parallel. 
Call this the matching property. 

Lemma 7.1 Suppose Pi,...,Pfc is a collection of positively convex polygons 

having pairwise disjoint interiors. Suppose also that each Pk intersects Q and 
has no spanning edges with respect to Q. Finally, suppose that the matching 
property holds for each vertex of Pk that lies in the interior of Q. Then 
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Proof: Let p he a vertex of Pk that lies in the interior of Q. The match- 
ing property, the positive convexity, and the matching property combine to 
show that the hnk of v is a circle. Hence, |J Pk contains a neighborhood of 
p. Let Q' = [jPk- If Q is not a subset of Q', then we can find some point 
q G int(Q) contained in an edge e of some Pk such that Q' does not contain 
a neighborhood of q. Since e is not a spanning edge, we can follow along 
e in one direction until we encounter a vertex of some Pk that lies in the 
interior of Q. Note that Q' cannot contain a neighborhood of any point on 
qp because e can only intesect other edges at vertices. (This comes from the 
disjoint interiors hypothesis.) Now we can say that Q' does not contain a 
neighborhood of p, which is a contradiction. ^ 

Now we turn to the case of polyhedra. Let Q be a golden polyhedron 
and let Pi, C Q be a finite collection of strictly convex golden polyhe- 
dra having pairwise disjoint interiors. Suppose also that every polyhedron in 
sight satisfies Equation [33l This means that slicing the polyhedra by golden 
horizontal planes results in golden polygons. In this section, we explain the 
test we use to show that Q is partitioned by {Pk}- 

The Slicing Method Let P^ denote the intersection of P with the hor- 
izontal plane z = h. As we just mentioned, Ph is a golden polygon for all 
h G Z[(j)]. Let S be the union of all the third coordinates of all the vertices 
of our polyhedra, including Q. 

We choose a subset S* C Z[(f)] that is interlaced with S. This means 
that there is one point of S* between each pair of consecutive points of 5*. 
For each h E S* , we use the 2 dimensional covering test to show that Qh is 
partitioned hj {{Pk)h}- 

To see why this test suffices, consider a vertex v in the polygon partition 
at height h. We would like to see that the link of v remains a circle as we 
continuously vary the height of the horizontal slice. Suppose some polygons 
separate as we increase the height of the slice. Then these same polygons 
overlap as we decrease the height from h. This cannot happen because of 
the disjoint interior condition. 

As we continuously vary the horizontal plane away from h, every slice we 
encounter satisfies the criterion for our covering test, namely every interior 
vertex has a circle link. So, for each open interval I & R — S , and any h E I, 
the polygons {{Pk)h} partition Qk- But this is clearly sufficient to show that 
{Pk} partitions Q. 
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7.5 Polyhedron Exchange Dynamics 

Here we explain how we compute various dynamical quantities related to the 
polyhedron exchange dynamics. 

Computing the Orbit Let E = {R/2Z) as usual. Let Pi,...,Pf5i be the 
polyhedra in the partition of S corresponding to the dynamics. Let Vi, 
be the corresponding vectors associated to each of these regions, so that 
^p + VkfoT pe Pfc. 

Here we describe one step of the polyhedron exchange dynamics. Given 
a point p G E, we use the interior point detection algorithm to find k such 
that p E Pk- We then replace p hj p + Vk and (if necessary) subtract even 
integers from the coordinates oi p + Vk until the result p' lies in the chosen 
fundamental domain for E. The map p — )■ p' is one step of the polyhedron 
exchange dynamics. 

In practice, we use floating point arithmetic to produce the canditate 
itinerary for a given point or region. By this, we mean the sequence of regions 
entered by the orbit of that point or region. Once we have the candidate 
itinerary, we switch to golden arithmetic to verify that the itinerary really is 
the correct one, as explained below. 

There is one flne point of our computer calculations which seem worth 
mentioning. The map from the Compactification Theorem involves divi- 
sion by 2. Since we want to work entirely in Z[(f\ we scale all vectors in 
by a factor of 2 when we perform the golden arithmetic calculations. We 
mention this for the benefit for a reader who wants to survey our computer 
code. 

Domain Verification: Let Pi, Pq^ and Vi, Vq4 be as above. Suppose 
that W is a vector and / = {ii, ik} is a sequence, with each ij in the set 
{1, 64}. We say that the pair (W, I) is feasible if 

• W lies in the closure of Pj^ . 

• W + Vi-^ lies in the closure of Pjj. 

• W + Vi^+Vi2 lies in the closure of Pjg. 

• And so on. 
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Suppose that Q is a polyhedron with vertices Wi, Wk- If (Wj, /) is feasible 
for all i = 1, ...,k, then the first iterates of are defined on all interior 
points of Q and / is the itinerary for such points. This works for the follow- 
ing reason. Since Q lies in the closure of Pjj , each interior point lies in the 
interior of Pjj . Since Q + Vi^ lies in the closure of Pi^ , every interior point of 
Q + Vi^ lies in the interior of Pjj. And so on. 

Maximality Verification: Suppose we are given a polyhedron Q and some 
integer N. Suppose that we have already verified that the first iterates 
are completely defined on the interior of Q. We would hke to check that the 
first N iterates are nowhere defined on dQ. 

We perform the following test for each face F of Q. Let W be some point 
in the interior of Q, as described above. We consider the orbit 

Wi^W; Wk^Wk-i + Vi,. 

and check that there is some k < N such that € dP^. 

7.6 Tile Creation 

Suppose that p G S is some point and is some integer. Here we explain how 
we produce the maximal convex polyhedron Q on which the first iterates 
of ^ are defined. Prom the viewpoint of making a rigorous verification, it 
doesn't matter how we produce Q. Once we have produced Q somehow, we 
use the domain verification algorithm to prove we arc correct. For what it is 
worth, we explain without proof how matter how we produce Q. 
We create Q in 4 steps. 

1. We produce the first A'" iterates of p using fioating point calculations, 
as above. 

2. We examine how each point of the orbit sits inside polygon of E that 
contains it. In the end, we produce a list of numbers, called the tile 
data, which carries the relevant information. 

3. We construct a polyhedron Q' from the initial point p and from the tile 
data. 

4. We replace the floating point coordinates of Q' by the best golden real 
approximations. The result is Q. 
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We emphasize that there is no guarantee that our construction will pro- 
duce the correct answer. It could fail at any stage. However, in practice, it 
works just fine, as later verified by golden arithmetic calculations. 

The second and third steps of the sketch above need more explanation. 

Pseudo Distances There is a list of 18 planes ni,...,ni8 such that each 

face of each is parallel to one of these planes. We have precomputed a 
certain list A^i, A^is of golden vectors such that is perpendicular to Nk 
for all k. Given a polygonal face F of one of the 64 polygons Pm, and a point 
q e -R^, we define the pseudo distance from g to F to be the quantity 

S{q,F)^Nk-{q-V). (66) 

Here k is such that F is parallel to H^, and F is a vertex of F. The quantity 
S{q,F) is independent of the choice of V. We choose the vectors {A^fc} so 
that p lies in the interior of Pm iff S{q, F) > for all F of P^- This is another 
way of expressing our interior point detection test. 

Tile Data: Let pi,...,pN be the forward iterates of p. There is a hst of 
18 planes Hi, His such that each face of each P^ is parallel to one of these 
planes. Initially, let D be the length 18 list {100, 100}. (The number 
100 here is just some convenient large number.) At the kih step of our pro- 
cess, we compute the pseudo-distance from pk to each of the faces in V^^. . If 
any of the numbers we get is smaller than the corresponding number in the 
list D, we replace this number of D by the new, smaller one. When we have 
done all N steps of the process, we have a final hst D. We call D the tile data. 

Tile Construction: Let D be the tile data associated to the pair {p,N). 
Call an index i active if Di < 100. If i is active, it means that some face 
of some Pj^ is parallel to Hj. For each active index, we let H- be the plane 
parallel to H such that the pseudo-distance from p to H- is Di. Let H- be 
the half-plane which contains p and has H^ as boundary. Then Q' is the 
intersection of all the H'^. 

In practice, we compute the intersection (J Hi as follows. For each triple 
of indices {i,j,k) we find the point H^ fl H^- fl H'^. We call such a point 
extraneous if some H^ separates it from the origin. Once we eliminate all 
extraneous intersections, the remaining ones are the vertices of Q' . 
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7.7 Tile Filling 



Our proof of the Renormalization Theorem involves partitioning the sets A 
and B into atoms. From the point of view of giving a proof, it doesn't matter 
how we produce these partitions. We simply have to verify that they work. 
However, it seems worthwhile describing what we actually did. 

Generating the Individual Atoms: We generate the B atoms in the 
following way. Given a point p E B, we compute the points "^^{p) for 
k = 1,2,3... until the orbit returns to B. Then we use the tile creation 
algorithm discussed in the previous section. 

We do not actually create the A atoms by this method. Once we have 
created all the B atoms, we use the action of R, the renormalization map, to 
define the A atoms. During the proof of the Renormalization Theorem we 
then use the domain verification method discussed in §7.5l to verify rigorously 
that the polyhedra we have defined really are A-atoms. 

Probing Vectors: We found by trial and error a list of 50 vectors having 
the following property. Let P be a special polyhedron and let f be a vertex 
of P. Then at least one of our 50 vectors, when based at v, points into P. 
Given that our polyhedron exchange map is defined by translations, the same 
property will hold with respect to any polyhedron - e.g. a periodic tile or an 
atom - produced by the dynamics. We call these vectors the probing vectors. 

Probing a List of Atoms: Suppose we have some list of fi-atoms con- 
tained in some branch Q (i.e. one of the 4 big convex polyhedra in each 
layer) of B. We choose some vertex v of one of the 5-atoms, or of Q. Let 
Pi, Pfc be all the atoms on our list that have f as a vertex. Say that a new 
probe is a probing vector W that does not point into and of the Pj. 

Given a new probe W, we choose a point v' so that v' — v is a small multi- 
ple of W, and then we find the P-atom that contains v' and add it to our list. 

Generating the Complete List: For each branch, we start with the empty 
list of B atoms. (The initial list of vertices is just the list of vertices of the 
branch.) We then iterate the probing algorithm, cycling through all vertices 
repeatedly, until we can find no new probes. This gives us our list of atoms 
filling the branch. We do this for all branches, and this gives us our complete 
list of P-atoms. 
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7.8 Variations on a Theme 



In §7.51 and §7.61 we explained some constructions that we perform relative 
to our polyhedron exchange map. In the next chapter, we will do similar 
things with respect to the square outer billiards map, and also with respect 
to an auxilliary map which we call the pinwheel map. We will not describe 
the corresponding operations in detail, because they are so similar, but we 
will give one example. 

Recall that \E' is the first return map to the strip S, as described in §4.21 
Suppose we want to verify that a given convex polygon Q is a maximal 
domain on which the first return map \E' is entirely defined, and that the 
return takes n steps. We first use floating point arithmetic to compute a 
(candidate) length n itinerary / for the center of mass of Q. We then check 
that (y, I) is feasible, in the sense of §7.51 for all vertices V of Q. Next, we 
check that ip^iV) C S for each vertex V of S. Finally, assuming that Q has k 
edges, we generate a list of k points f i, f^, such that Vj lies in the interior 
of the jth edge of P and for each j there is some m = rrij e {1,...,?t, — 1} 
such that ip"^{vj) lies in the boundary of the relevant domain of definition 
for ip. 
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8 The Pinwheel Lemma 



8.1 The Pinwheel Strips 

We are interested in the first return map \1/ : Il± Il± discussed in §4.21 
Our Pinwheel Lemma below equates the map with a new map 11 that is 
simpler to manage. As a step in the Compactification Theorem, we prove a 
result which we call the Pinwheel Lemma. The Pinwheel Lemma equates 
with a map LI which is easier to analyze and "compactify" . 

We first explain how to associate 4 pairs {Si, Vi) for i = 0, 1, 2, 3 to the 
kite K. Here Si is an infinite strip and Vi is a vector that points from one 
component of dSi to the other. 




Figure 8.1: The strip associated to Cq. 
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Wc label the edges of K as in Figure 8.1 We call these edges Cq, ei, 62, 63. 
We orient these edges so that they go clockwise around K. The labels are 
such that the lines extending the edges hit the "circle at infinity" (in the 
projective plane that compactifies R^) in counterclockwise cyclic order. 

To the edge ej we associate the strip Sk that has the following properties. 
One boundary component of L^. of Sk extends e^. The other boundary com- 
ponent is L'j^ is such that the vertex Wk of K farthest from lies halfway 
between Lk and L'^. So, K extends exactly up to the centerline of Sk- We 
define 

Vk^±2{vk-Wk). (67) 

Here Vk is the head vertex of e^. We choose the signs so that the vectors are 
as shown in Figure 8.2 below. Figure 8.1 shows 5*0 and Vq. 
Let A — (f)'^. Here are the formulas for these strip pairs. 

• Vq = (2, 2) and 5*0 is bounded by the fines x + y + 1 G {0, 4}. 

• Vi — (0, 4) and Si is bounded by the fines x — y — 1 e {0, 4}. 

• (2, -2) and ^'2 is bounded by the fines x - Ay - A e {0, -4/0}. 

• ^3 = {-2-2A, 0) and 6*3 is bounded by the lines x+Ay-A e {0, -4/0}. 

Figure 8.2 shows a rough picture of the strips and vectors relative to 
K[l/4). The picture is very close to what one sees for K{(f)^'^). The black 
kite in the middle is i^(l/4). Figure 8.2 also shows the strip T, = Rx [—2, 2], 
as well as the regions E+ and E_. These regions are denoted (+) and (— ) 
respectively. 

To make our later definitions cleaner, we define recursively define 

V^+i = -Vi, S4+i = Si. (68) 

This gives us pairs (Si.Vi) for all i = 0,1,2,3,4... The strips repeat with 
period 4 and the vectors repeat with period 8. However, note that {S, V) 
and {S, —V) define the same strip maps. So, the strip maps repeat with 
period 4. Nonetheless, it is useful to distinguish between Vi and — —Vi, 
etc. This way, the vectors Vi,...,Vs can be drawn so that they naturally 
circulate counterclockwise around K, as in Figure 8.2. 
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Figure 8.2: The strips and vectors associated to K. 
8.2 The Main Result 

Let S be an infinite strip in the plane. We say that a strip pair is a pair 
{S, V) where 5" is a strip and y is a vector that points from one component 
of dS to the other. Given the strip pair, we have an associated strip map 
E : ^ S, given by the formula 

E(p) ^p + nV eS. (69) 

Here n — n{p) is chosen so that E{p) e S. The map E is well-defined on 
the complement of an infinite discrete set of lines parallel to the boundary 
components of V. Note that the vectors V and —V define the same map. 
The pinwheel map is essentially the composition of strip maps. 
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Referring to the strip pairs {Si, Vj) defined in the previous section, we 
recursively define 5*4+4 = Si and V4+j = V^. This gives us pairs {Si,Vi) for all 
integers i > 0. Of course, these pairs repeat with period 4. We let Ei denote 
the strip map associated to {Si, Vi). This gives strip maps Ei. We let ( be 
the strip map defined relative to the pair (E, (0, 4)). We define the pinwheel 
map 

n^CoEso ...oEi-.E^E. (70) 

Recall that ^ is the first return map on E+ and likewise the first return map 
on Sj. 

We proved versions of the Pinwheel Lemma in [SI], [S2], and [S3], but 
unfortunately none of these results gives us the exact statement we need. 

Theorem 8.1 (Pinwheel) A point of U S_ has a well-defined U-orbit 
if and only if it has a well-defined "if -orbit, and 11 = ^ on all points which 
have well defined orbits. 

Let E5_ denote those points in E± of the form {x, y) with < 20. We see 
by direct inspection that the Pinwheel Lemma holds for all points in Ej- — E*. 
Both maps just circulate the point around the kite, as shown in Figure 8.2. 
See any of the papers cited above for more details about this. 



Figure 8.3: Part of the covering 

To deal with E^, we make a direct calculation. We produce a list of 
572 golden convex polygons Poi---i-P57i whose union covers E^. Figure 8.3 
shows some of these polygons. The central dark polygon is the Penrose Kite. 
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The two triangles which share the top and bottom vertex with K belong to 
S — E-t, and our dynamical statements do not apply to them. We do not 
consider them to be on our list. 

We use the method described in §7.81 to establish the following result for 
each polygon P on our list. 

• \E' is entirely defined on the interior of P. 

• is nowhere defined on any non-horizontal edge of P. 

• \1/ is the identity on the interior of any horizontal edge of P. 

As discussed in §7.4[ this result implies that our polyhedra have pairwise 
disjoint interiors. We then use the method in §7.4l to show that the rectangle 
[—20, 20] X [—2, 2] is covered by the union of our polygons and the Penrose 
kite and the two special triangles. This implies that our union of polygons 
covers T,*^UT,*_. 

Next, we computationally establish exactly the same results for 11. Fi- 
nally, we check that 11 = \E' for one interior point of each polygon. Since 
n and \l/ are entirely defined on each polygon, we see that 11 = \l/ on the 
interiors of each of our polygons, and also on the interiors of the horizontal 
edges. None of the remaining points have well-defined orbits for either map. 

8.3 Discussion 

For the purposes of making a rigorous proof, it doesn't really matter how we 
generate our polygons. However, it seems worth saying a word about how 
we do it. First of all, we use the method described in §7.6[ with in place 
of to generate the individual tiles. 

To generate the collection of polygons, we let N be the set of points of 
the form 

(lUo 1^ + ')' ^ = -2500, ...,2500; j = l,...,199. (71) 

The small "fudge factor" e = 10~^ is present to guarantee that we choose 
points on which the map 11 is defined. For each p & N we compute P{p), 
and then we weed out redundancies and remove the ones that are disjoint 
from E^. 
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9 The Compactification Theorem 



9.1 Discussion 

Recall that S = i? x (—2, 2] and 11 : E — )■ E is the pinwheel map. In light of 
the Pinwheel Lemma, it suffices to prove the Compactification Theorem for 
n in place of We first compare the Compactification Theorem with the 
corresponding results in [SI] and [S2]. 

Let El C E denote the union of two horizontal lines R x {±1}. Let Ei 
denote the 2-torus slice of the 3-torus E, obtained by interecting E with the 
plane of height 1. See Figure 4.6. The map n carries Ei to itself. 

Our result [SI, Arithmetic Graph Lemma] is equivalent to the restricted 
version of the Compactification Theorem, when Ei replaces E, and Ei re- 
places E. However, the result in [SI] has a somewhat different emphasis. 
There we were concerned with how Ei controls the structure of the corre- 
sponding arithmetic graphs. In [S2, Master Picture Theorem] we proved 
a result that is equivalent to [SI, Arithmetic Graph Lemma] for all kite 
parameters. The result in [S2] is stated in the same general terms as the 
Compactification Theorem. 

The proof in [S2] is different than the proof in [SI]. It is more conceptual 
and also much easier to generalize. In spite of this, we will prove the Com- 
pactification Theorem following the ideas in [SI]. We do this partly because 
we think of this paper as a sequel to [SI], and also because we would like to 
take the opportunity to revisit that proof and give a cleaner exposition. The 
proof we give is specially adapted to the Penrose kite, and would not easily 
generalize. 

9.2 The Covering Property 

The first order of business is to show that the 64 golden polyhedra defining the 
map really do partition the torus E into strictly convex polyhedra. (This 
would be the first step in any approach we took to proving the main result.) 
We use the tests described in §7.31 to show that the polyhedra involved are 
strictly convex and have pairwise disjoint interiors. We also verify that they 
satisfy Equation [33 Finally, we use the slicing method described in §7.41 to 
show that our polyhedra really do partition E. 
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9.3 The Proof Modulo Computations 

Let F be the fundamental domain for S that we defined in §5.21 Recall that 

n = c o Eg o • • • ^1 (72) 

where Ei, ...,Es and ( are the strip maps defined in §8.21 In the next section, 
we define the following objects. 

• An affine map Q : ^ R"^. See Eq. [H2] 

• A locally affine map C : ^ S. See Eq. M 

• An affine embedding It : F ^ R^, for any t G ^[0]. See Eq. IHH 

• A set Lk C R^ of parallel hyperplanes, for = 1, 8. See Eqs. [85H86] 

• Piecewise affine maps Ei, Eg : R"^ — )• R'^. See Eqs. I57H551 

Associated to Ej is an infinite family Lj of parallel lines, such that Ej 
is defined exactly in the complement of the lines. At the same time, Ek is 
defined precisely in the complement of Lk. We write ~ ^ is V—W G {4Z)'^. 

Lemma 9.1 The following is true. 

1. ^o& = &o(. 

2. CoT = C 

3. ( o If. = Identity for all t. 

4. Eio It is independent oft. 

5. T{Lk) = Lk. 

6. e-\Lk) = Lk. 

7. EkoT^Ek. 

8. Qo Ek EkoQ. 

Here T is and element of (AZ)^, considered as a map of R"^, and k = 1, 8. 
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For the next result, we interpret Q{p) as a point of F. When the point 
hes in the interior of F, it is the unique representative of Q{p). Otherwise, 
we make some choice amongst the possibihties. 

Lemma 9.2 There is a A-element set Q C Z[(f)\ with the following property. 
For any p G S there is some t ^ Q such that that It o B(p) ~ 0(p). 

Define the following piecewise (locally) affine map. 

IL = CoEso ...Ei: R'^ ^t. (73) 

Lemma 9.3 Let P be a special polygon. H o 1^ is everywhere defined and 
locally affine on the interior of P, but nowhere defined on the boundary of P. 

We prove the above 3 results in the sections following this one. 

Lemma 9.4 Suppose that p lies in the interior of S and Q{p) lies in the 
interior of a special polyhedron. Then U is well-defined on p. 

Proof: Let t G Z[(f)\ be the value guaranteed by Lemma [9. 2^ and let / = If. 
We define ^ 

go = /oe(p); Po = eip). (74) 

By Lemma [12] we have po ~ go- We inductively define (jk = -Efc(gfe_i) for 
k = 1, ...,8. Applying Lemma to the open special polygon that contains 
Q{p), we see that the points qi, gs are well defined, meaning that qk-i ^ L^. 

Since Li is (4Z^)-invariant and po ~ qo, we have po ^ Li. Hence, we may 
define pi = Ei{po). By Statement 7 of Lemma |9TT| we have pi ~ gi. Shifting 
the indices and repeating the same argument 7 more times, we find that we 
can inductively define points pk = E^ijik^i) for k = 1,...,8. That is, 11 is 
defined on p. Put another way, Pk-i ^ for = 1, 8. 

Let Po = p. We have po ^ Li. But then 6(po) ^ Li because Li is {AZY- 
invariant and Q{po) ~ Pq- By Lemma |9.H we have p^ ^ Li. Hence Ei is 
defined on po- This allows us to define Let pi = Ei{po). Shifting the indices 
and repeating the same argument 7 times, we can inductively define the 
points pk = Ek{pk-i) for = 1, 8. In short, the composition Eg o ... o Ei 
is defined on p. 

Finally, the set of lines R x {2Z) is invariant under all our strip maps. 
Since p does not he on these lines, neither does ps. Hence ( is defined on pg. 
Hence, H is defined on p and H(p) = ({Ps)- ^ 
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Lemma 9.5 no0 = no0 whenever both sides are defined. 

Proof: We write Vi ~ V2 when Vi - V2 e (4Z)^. Let p G S. Repeatedly 
using Statements 7 and 8 of Lemma 19.1^ we have 

Eg o . . . o El o Q{p) ~ o Eg o . . . o Ei{p). (75) 

By Statement 2 of Lemma I9.H we get equahty above when we apply ( to 
both sides of Equation [7S1 That is, 

n o e(p) = CoQ o Eso... Ei{p) =* eo(oEso...o Ei{p) = eo u{p). (76) 

The starred equality is Statement 1 of Lemma 19.11 ^ 

Lemma 9.6 Suppose pi,p2 lie in the interior of H and 0(pi),0(p2) in 
the interior of the same special polyhedron. Then o n(pj) — Q{pj) is the 
same for j = I and j = 2. 

Proof: Let P be the special polygon whose interior contains Q{pi) and 

e{p2). 

no4oe(pi)-e(pi) = no/,^oe(p2)-0(P2) = no4oe(p2)-0(P2). (77) 

Here tj = t{pj), in the sense of Lemma \9.2\ and the subtraction takes place 
in the abelian group {R/2Z)^. The first equality comes from the fact that 
n o /^^ is locally affine and entirely defined on P. The second equality has 
the following explanation. By Statement 4 of Lemma 19. the map Ei o is 
independent of t. Hence H o is independent of t. 

Noting that C(Pj) = Pjy we use the identities in Lemmas 19.11 and 19.51 to 
conclude that 

e o H(p,) - e(p,) = e o u{p,) - e o c{p,) = n o e(p,) - c o e{pj) = 
H o i,^ o eip,) - Co o Qipj) = H o i,^ o e(p,) - e(p,), (78) 

Combining Equations [77] and [7H1 we get the equation we seek. ^ 
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Recall that \1/ is a piecewise translation on S. For each of the 64 polyhedra 
P that partition S, we produce a point p G S such that 

e{p) G interior (P); 6 o n(p) = (79) 

Indeed, this is how we defined \l/ in the first place. If g G S is any other point 
such that Q{q) lies in the interior of P then, by Lemma [9.61 we have 

e o n(g) - e(g) = e o u{p) - e(p) = ^{p) -p = ^{q)-q. 

This shows that G o n = \I' o B for all points p G S such that Q{p) lies in the 
interior of a special polyhedron. 

To finish the deduction of the Compactification Theorem, we just have 
to understand what happens in exceptional cases. 

Lemma 9.7 Suppose that p lies in the interior of T, and <d{p) does not lie 
in the interior of a special polygon. Then U is not defined on p. 

Proof: We will revisit the proof of Lemma 19.41 and use the notation set up 
in that proof. We will suppose that 11 is defined on p and then derive a 
contradiction. Since 11 is defined on p, the compositions Ei, E2 oEi, etc. are 
defined on p. Hence, the points pi, ...,ps are well defined. But then the same 
argument as in the proof of Lemma 19.41 (but done in reverse) shows that the 
points pi, ...,ps are alljiefined. Hence, the maps Ei, ...,Es are all defined on 
Eq. Finally, the map C is everywhere defined on R^. Hence H is defined on 
Po, a point in the boundary of Pq- This contradicts the second statement of 
Lemma [9731 4 

It remains only to deal with the points on (9S. The only issue is the action 
of the map (, and this only comes up at the end of our proof of Lemma [9.41 
The problem is that ( is undefined on the set R x 2Z. These lines are 
invariant under the action of Ei, Hence, if p G dT,, then ps & Rx 2Z 
and ( is not defined on pg. However, recall that S = (—2, 2] x R. The bottom 
boundary is left off. Given this convention, we define 

C(p8) = (x8,2)GaS. (80) 

Here is the first coordinate of pg. Once we make this definition, our proofs 
above go through without a hitch. 
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9.4 Proof of Lemma 19.11 

Recall that 

e{xi,X2)= (l' ^'O) + ~^"Y~^'^2) ^ C(a;i,a;2) = (xi, [x2]4)- 

(81) 

Here [V]2 means that we take the image of V in {R/2ZY, and [^2)4 means 
representative of X2 in R/AZ that lies in (—2,2]. 
We define 

0(^1, 0:2) = hEi + X2,Xi - 0:2, — 7 — , — 7 — . (82) 



This is exactly the same map we used in [SI]. Next, we define 

C{xi,X2, X3, X4) = (l' ^' 0) + ^(^3 + X4, X2, xi - X2) e {R/2Zf. (83) 
For any t G -Z'[0] we define We define 

/t(xi, X2, Xs) = (2?/2 + 2/3, 2?/2 + 2/3, Z/l, 2/l) + (2/3, -2/3, — " i^, -— + ^), 

(2/1, 2/2, 2/3) = (a;i, X2, X3) - (^1, ^, oj. (84) 

Now we define the families of parallel hyperplanes. We introduce the 
functions ■ R^ R as follows. 

5fi = a;2 + l; g2 = X2+X3 + 2+(f); gs = xi + Xi + 2+(j); 5(4 = + (85) 

We define 

Lk+4 = Lu = 9k\^Z). (86) 

Statements 1-4: 

1. A direct calculation shows C°0 = © and = QoC^. Hence C°0 = 0°C- 

2. It is immediate from the formula that C o T = C for any T G {AZ)^. 

3. A direct calculation shows that ( o is the identity. 

4. It is immediate from the equations that is (4Z'^)-invariant. 
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For X G -R — 4Z let r(x) denote the greatest element in 4Z' that is less 
than X. Define 

7fc(X) = X (r o G 40Z. (87) 

The (x) in Equation [571 means multiplication. Here X = (xi, X2, Xs, X4) is a 
point of — Lk- Finally, we define 
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73 
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-74 

V / 



(88) 



Here we have simplified our notation slightly by writing E\ = Ei{X) and 
7i = 7i(-^)) 6tc. Each E^ is the sum of a linear transformation and a piece- 
wise constant vector- valued function. The linear part is defined on all of R^, 
but the vector-valued function changes valued when one passes through a 
hyperplane of L^. We define ii^4+fc = Ek. 



Statement 5: The maps Jq and It agree up to post-composition with the 
translation form V ^ V + {0,0, —t, t). Inspecting the formula for Ei, we see 
that the only expression involving the third and fourth coordinates is X3 + X4. 
Hence Ei o is independent of t. Hence H o is independent of t. 
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Statement 6: We can deduce the formulas for the hues C by looking 
at the formulas for the strip pairs given in §8.11 We have L^^^ = L^, so we 
just need to consider Lk for k = 1,2, 3, 4. Let A = 

• Li consists of those lines satisfying xi — X2 — 1 G AZ. 

• L2 consists of those lines satisfying xi — Ax2 — AG A(f)~^Z. 

• L3 consists of those lines satisfying xi + Ax2 — A ^ A(f)~^Z. 

• L4 consists of those lines satisfying xi + X2 + 1 € AZ. 

The fact that 0~^(Lfc) = is obvious for k = 1,4. For k = 2 we compute 

g2 o 6(xi, X2) = (pxi — (p~'^X2 + 4 — = — Ax2 — A). 
The result is obvious from here. The proof for /c = 3 is similar. 

Statement 7: Let 61,62,63,64 be the standard basis vectors of R'^. let 

u{j, k, X) = E,{X + 46fc) - Ej{X). (89) 

Statement 7 is equivalent to the statement that u{j, k, X) G (AZ)^ for all j, k 
and X. 

We see directly that u{l, k, X) G (4Z)^ for A; = 1, 3, 4, and we compute 
u{l, 2, X) = (0, 0, -4/0, 4/0) + (0, 0, 40, -40) = (0, 0, -4, 4). 

This proves the result for Ei. The proof for E4 is similar. 

We see directly that u{2, k,X) G (4Z)^ for j = 1,4, and we compute 

u{2, 2, X) = (0, 0, 0, 4/0) + (0, 0, 0, -40) = (0, 0, 0, -4), 

u{2, 3, X) = (0, 0, 0, 40) + (0, 0, 0, -40) = (0, 0, 0, 0). 

This proves the result for E2. 
For E^, we compute 

w(3, 1, X) = (40, 4/0, -4/0, -40) + (-40, -40, 40, 40) = (0, -4, 4, 4) 

m(3, 2, X) = (-4, 0, 4, 4); m(3, 3, X) = (-4, -4, 8, 4). 
«(3, 4, X) = (40, 40, -40, -4/0) + (-40, -40, 40, 40) = (0, 0, 4, 4). 
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Statement 8: We want to prove that 

QoEk^EkoQ, 

where ~ denotes equivalence mod (4Z^)^. Recall that the strip map E^. is 
defined in terms of a pair {Sk, Vk), where Sk is a strip and is a vector. The 
two lines of dSk are consecutive lines of L^. Hence, by Lemma \9.1\ the set 
Sk = Q{Sk) lies between two hyperplanes in L^. Hence E^ is an affine map 
on Sk- Indeed, Ek is linear on Sk because 

7(0, 0, 0, 0) = 0; (0, 0, 0, 0) = 9(0, 0) G 6(5^) = Sk 

A direct calculation shows that the linear part of Ek is the identity on the 
2-plane G(i?^). Hence Ek is the identity on S. In summary 

Q{p) = EkO Q{p) Vj9 G interior(E). (90) 

Let Vk = Q{Vk). Given the piecewise affine nature of our maps, the value 
of Ek{X + Vk) — Ek{X) is independent of the choice of X G R^ — Lk- Compare 
the proof of Statement 7. A direct calculation shows that EkiVk) G {4Z)^. 
Therefore 

EkiX + Vk)r^ EkiX); \/XeR^-Lk. 
Iterating this formula, we get 

EkiX + nVk)--EkiXy, yXeW-Lk. (91) 

Here n is an integer. 

Every point p' G — Lk has the form p + nVk for some p in the interior 
of Sk- Note that Ek{p') = p. We have 

QoEkip') = eip)=EkoQ{p)r^ 

Ek{Q{p) + nVk) = Ek o Q{p + nVk) = Ek o e{p'). 
This proves Statement 8. 
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9.5 Proof of Lemma 19.21 

Let p = (xi, X2) G S. When we interpret 6(p) as a point of the fundamental 
domain F, we are choosing integers Ai, A2, such that 

e{p) = (2A, + l,2A,+ ^-,2As) + Q,^^^,x,) e F. (92) 

Note that 

^3 e {0, 1} (93) 

because X2 G (—2,2]. Using Equation [921 and the definition of the map It, 
we compute 

~ / 2yl 2yl \ 

It o e(p) = e{p) + (4A2 + 4^3, 4A2, 2A, + ^ + t, 2A,-^-ty (94) 

The appropriate choice of t from amongst the values 

0; 2; -2/0; 2-2/0 (95) 

leads to 

It o e(p) ~ e(p) 

This proves Lemma [9.21 

9.6 Proof of Lemma 19.31 

The arguments we give here are very similar to what we discussed in §7.5[ 
but the setting is different enough that we will go through the details. 
Given the 64 special polyhedra Pi, Pq^, we define the 256 polyhedra 

P,, = hXPj) C R'; {1,2,3,4} j G {1, 64}. (96) 

Lemma 9.8 11 is defined on the interior of Pij for all i = 1,2,3,4 and all 
j = l,...,64. 

Proof: Let q be the center of mass of Pij. Using fioating point arithmetic, 
we associate to g a certain 8-tuple of integers, as follows. Let qo = q and 
inductively define qk = Ek{qk-i)- We define 

nk = (t^^^lkiVk-i)- (97) 
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This integer locates the component of — that contains Qk-i- The 8- 
tuple of interest to us is (ni, ng). For the purposes of our proof, it doesn't 
matter how we generated the sequence (ni, ...,^8). This sequence just serves 
as a guide for our rigorous calculations. ^ 

Guided by the sequence we have produced, we define the map El to be 
the same map as Ek, except that 40^^ replaces in the formula. We check, 
using golden arithmetic, that the composition E^o . . . E^ is defined on all the 
vertices of Pij. For each such vertex vq we inductively define Vk = Elivk^i). 
We check that 

9k{vk) e [nfc,nfc + 4]. (98) 

Here gk is as in Equation [H5i 

The same argument we gave when we discussed domain verification in 
§7.51 now says that E^ o Ei is defined on the interior of Pij. Indeed, we can 
now say rigorously that the sequence (ni, ng) "works" for all points in the 
interior of Pij. 4 

Lemma 9.9 11 is not defined on any boundary point of Pij for any indices 
i = 1, 2, 3, 4 and any j = 1, 64. 

Proof: We use the notation from the previous proof. For each face F 
of Pij we exhibit an index k such that one of the following two equations 
simultaneously holds relative to any vertex vq of F: 

• Qkiyk-i) = Uk- 

• 9k{vk~i) = rtk+A- 
This shows that 

Fk-i = Ek-io...A{F) 

lies entirely in one hyperplane of Lk. Hence Ek is nowhere defined on Fk-i. 
Hence H is nowhere defined on F. Since F is an arbitrary face, H is nowhere 
defined on the boundary of Pij. 
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10 Proof of the Renormalization Theorem 



10.1 Some Terminology 

We use the notation from the Renormahzation Theorem. Here we estabhsh 
some terminology which will be useful to us in this chapter and in later ones. 

Forweird Atoms: Recall that A is partitioned into atoms. We call these 
atoms forward A-atoms. We define the forward S-atoms in the same way. 

Neutral and Backward Atoms Let P be a forward A atom. There 
is some n > 1, depending on P, such that ^'^(P) is disjoint from A for 
k = l,...,n - 1 but $"(P) C A. We call the polyhcdra $'=(P) neutral A- 
atoms when k — l,...,n — 1. We call \I^"(P) a backward A-atom. A is also 
partitioned into backward A- atoms, and the first return map of "^f"^ is en- 
tirely defined and a translation on each one. 

Chains: Let P be a forward A-atom. We call the sequence P, .... ^"^"^(P) 
a forward A-chain. So, each forward A-chain is a finite union of A-atoms, 
the first of which is forward and the rest of which are neutral. We make the 
same definitions for B. 

Backward Chains: We define backward chains just as we defined the for- 
ward chains, except that we use backward atoms in place of forward atoms, 
and we use in place of ^. We think of the backward atom as being the 
last atom in the backwards chain. For each forward chain there is a back- 
ward chain that shares all its atoms except the first one. Likewise, for every 
backward chain there is a forward chain that shares all of its atoms except 
for its last one. 

Types: We distinguish two types of layers of A. The type-1 layers are 
the ones on which the restriction of R to that layer is a translation. The 
other layers we call type-2. Note that this definition is consistent with the 
notion of type defined in association with the Reduction Theorems. We de- 
fine the type of a A-atom and a A-chain according to the type of the layer 
that contains it. 
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10.2 Statements 1, 2, and 5 

We compute the list of A and B atoms using the method described in §7.71 
Figure 10.1 shows a shce of the set of 5-atoms at a particular height. 




Figure 10.1: The fi-atoms shced at -7 + 40 ^ .5278 

We use the method described in §7.5l to check that each of the polyhedra 
on our candidate list of A-atoms is indeed a A-atom. We do the same thing 
for the i?-atoms. In the case of the 5-atoms, we also verify that each face F 
of a S-atom has one of 3 properties. 

• The first return map is not defined on F. That is, some small 
power of \E' maps F into one of the faces of one of the 64 special poly- 
hedra that define \Ef. 

• F lies in dB. 

• The first return of F to -B lies in dB. 

Since these properties hold for every edge, all the B atoms are maximal. 
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Remark: At some point we also checked that all the A-atoms are max- 
imal, but we did not save the computer code for this. This check is not 
necessary for our arguments. 

Given that the i?-atoms are all maximal, they have pairwise disjoint inte- 
riors. We also check that they are strictly convex. We now apply the slicing 
method of §7.41 to each branch of each layer of B. We conclude that the 
5- atoms partition B. That is, we have obtained the complete list. 

As we mentioned in §7.7[ we obtain the candidate list of A atoms by 
pulling back the B atoms by the action of the renormalization map R. From 
this way of constructing things, we automatically know that the A-atoms 
are strictly convex and have pairwise disjoint interiors. (We don't need the 
maximality for this.) We now apply the slicing method to A just as we did 
for B. We conclude that A is partitioned into the A atoms on our list. 

We have set things up so that, on each type-1 layer, R establishes a 
bijection between the forward A-atoms and the forward i?-atoms. We check 
that this bijection is compatible with action of the first return maps \l/|y4 and 
Since we have already verified that everything in sight is an atom, we 
just need to check the action of the relevant maps on one interior point of 
each atom. For the check, we choose interior golden points as in §7.31 

We have set things up so that, on each type-2 layer, R establishes a 
bijection between the forward A-atoms and the backward 5-atoms. In the 
same way as in the type-1 case, we check that the action of R is compatible 
with action of the first return maps "^{A and 

These calculations prove Statement 1 of the Renormalization Theorem. 

Statement 2 of the Renormalization Theorem is really just a description 
of the action of R. We have constructed R so that it acts on the horizontal 
slices as R acts on R/2Z. 

Statement 5 is purely combinatorial. We simply check, in each case, that 
each i?-chain has a shorter length than the corresponding i?-chain. 

Remark: For the record, we mention that there are 678 B atoms and 
3 X 678 = 2034 A-atoms. We also compute that the longest A-chain has 
length 703 and the longest 5-chain has length 109. In other words, for 
generic points in A, the first return map is defined in at most 704 iterates 
and for generic points in B the first return map is defined in at most 110 
iterates. 
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10.3 Statement 3 



When we lift to we see that it consists precisely in the planes 

spanned by the vectors 

(20-\l,O); (0,-2,1) 

which intersect the line {0} x Rx {0} in points of the form (0,t,0) where 
t G ^ [(/>]. This follows from a routine calculation. Another calculation shows 
that each choice of R mentioned above preserves this set. 

10.4 Statements 4 for B 

Recall that B is partitioned into 6 layers, as discussed in §5.41 It turns out 
that some periodic tiles intersect some layers of B but not others. To make 
our proof go cleanly, we introduce the notation of a -B-periodic tile. A B 
periodic tile is an intersection of the form 

P n (^{R/2Z) X {R/2/Z) X Ja) (99) 

where Jx is the 6 intervals Ji, Jg in the partition defining the layers of B. 

Using essentially the same method as in §8.31 we generate a candidate list 
of B periodic tiles whose orbits avoid B. In practice, we simply enumerate 
all the periodic tiles whose orbits avoid at least one layer of B, and then we 
keep track of which layers each orbit intersects. We check that the period of 
the longest tile on the list is 33. 

Let A G {1, 2, 3, 4, 5, 6} be an integer that describes a layer of B currently 
of interest to us. Let i G {1, 64} be an index that describes one of the 64 
special polyhedra that partition E. 

Let B{i, A) be the union of all the -B-tiles and forward/neutral 5-atoms 
that are contained in the region 

P(z, X) = Pin (^R/2Z) X {R/2/Z) X Ja)) (100) 

By construction, these tiles have pairwise disjoint interiors. Say that (i. A) is 
good if the tiles in B{i, A) partition the region in Equation llOOi 

Lemma 10.1 Statement 4 holds for B provided that every pair (i. A) is good. 
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Proof: Let pi, ...,pio9 be a (non-periodic) generic orbit portion of length 109. 
Consider the point p = piog. There is some pair (i, A) such that p G P{i, X). 
Since {i, A) is good and p is generic, p lies in the interior of either a B periodic 
tile or in the interior of a forward/neutral i?-atom. In the former case, p has 
period at most 33. This is a contradiction. If p lies in the interior of a forward 
5-atom, then p E B,hj definition, and we are done. 

So, suppose p lies in the interior of neutral i?-atom. Since the longest 
5-chain has length 109, there is some m < 109 such that G B. But 

then pio9-m G B. This covers all the cases. 4 

Our computer code is equipped to verify that all indices are good, but 
we can get away with doing less computational work. The indices {i, A) 
are automatically good for i G {1,...,22}. These correspond to the special 
polyhedra on which acts trivially. So, we only have to consider indices 
i>22. 

We can eliminate more indices by considering the dynamics of \I'. For 
instance, \I'(P26) = As- Hence, if (25, A) is good, then so is (26, A). In 
Lemma 114. 1^ we prove that the list of indices 

/ = {23, 25, 32, 40, 41, 44, 46, 53, 61, 62} (101) 

is such that every well-defined and nontrivial orbit intersects Pi for some 
i E I. For this reason, we only have to check that {i, A) is good for i G / and 
AG {1,2,3,4,5,6}. 

There is one more savings we can make. Everything in sight is invariant 
under the involution 

ix,y,z)^ (2,2,2) -{x,y,z). (102) 

This involution permutes the special polyedra, commutes with preserves 
A and B, and permutes the list of atoms. For this reason, it suffices to prove 
that the indices {i, A) are good for i G / and A G {1, 2, 3}. Using the shcing 
method of §7.41 we check that this is the case. This establishes Statement 4 
for B. Incidentally, our calculaton show that our candidate list of 5-periodic 
tiles is correct and complete. 

Remark: The calculation we make is fairly massive, but the reader can 
use our program and see the individual slices plotted as they are tested. The 
reader can also select any index [i, A) and plot the corresponding partition. 
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10.5 Statement 4 for A 

In principle, we could prove Statement 4 for A just as we proved it for B. 
This time there are 18 layers. The problem with this approach is that it 
leads to a massive calculation. For instance, the partition associated to the 
index pair (23, 7) has 2562 polyhcdra. 

Were we to take this approach (and succeed) we would know than any 
generic orbit portion of length 703 intersects A. The bound of 703 is the 
sharp constant. There is an orbit portion of length 702 that avoids A. We 
will take a different approach and get the slightly worse constant 812. The 
benefit we get from the other approach is that the calculation we make is 
vastly shorter. The idea is to use the result we have already proved for B 
and only calculate "the difference" between B and A, so to speak. 

Recall that each layer of B is divided into 4 branches. Moreover, each 
layer of B is partitioned into finitely many smaller pieces, corresponding to 
the layers of A. We call these pieces of B the sub-layers of B. Precisely, 

• Layer 1 of S has 5 sublayers. 

• Layer 2 oi B has 1 sublayers. 

• Layer 3 of S has 3 sublayers. 

• Layer 4 of S has 3 sublayers. 

• Layer 5 of S has 1 sublayers. 

• Layer 6 of B has 5 sublayers. 

The palindromic nature of the list comes from the involutive symmetry men- 
tioned in the previous section. 

Each layer of B is divided into 4 branches. Each branch is a golden 
polyhedron. By slicing each branch with relevant horizontal planes, we can 
say that each sublayer of B is also divided into 4 branches. In summary, B 
has 18 sublayers and each sublayer has 4 branches. We index the sublayers 
by A G {1, 18} and the branches by /3 G {1, 2, 3, 4}. 

Finally we discuss A. We define the A-periodic tiles just as we defined 
the S-periodic tiles. We find a candidate hst of all A-periodic tiles as we did 
in the B case. The longest period in this case is 213. We say that a pair 
(A, /3) is good if the polyhedron corresponding to (A,/?) is partitioned into a 
union of A-periodic tiles and forward/neutral ^-atoms. 
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Lemma 10.2 Statement 4 holds for A provided that every index (A, /3) is 
good. 

Proof: the same argument as in Lemma [10.11 proves that any orbit portion of 
length 703 that starts in B intersects A. Combining this result with Lemma 
IIO-H we see that every generic orbit portion of length 812 = 703 + 109 inter- 
sects A. 4 

To finish the proof of Statement 4 for A we just have to prove that every 
index (A, (3) is good. Using the involutive symmetry discussed in the previous 
section, it suffices to prove this result for indices A < 9. Using the slicing 
method of §7.41 we verify that the indices (A,/3) are good for A G {1,...,9} 
and /3 G {1, 2, 3}. This completes the proof. 

Remark: This time, our calculations do not verify that our list of peri- 
odic A tiles is complete. Rather, all we know is that each polyhedron on the 
list is indeed a A-periodic tile, and that our list of yl-periodic tiles whose 
orbits intersect B is complete. We can still say something useful, even with 
an incomplete list. Statement 4 for A implies that S is partitioned into the 
forward/neutral A-atoms and a finite union of A periodic tiles, some of which 
perhaps are not on our list. If we have a generic orbit portiont that is not 
contained in a periodic orbit, then it must intersect one of the forward or 
neutral A-atoms. But the longest A-chain has length 703. this establishes 
that a generic orbit portion of length 703 intersects A provided it does not 
lie on a periodic orbit. The constant 703 is sharp because of the A-chain of 
length 703. 

10.6 Statement 6 

Any generic orbit that avoids B is contained in the orbit of one of the B- 
periodic tiles on our list. We check that all these orbits avoid A. Hence, any 
orbit that avoids B also avoids A. 
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11 Critical Renormalization Calculations 



11.1 Discussion 

The goal of this chapter is to prove the Fundamental Orbit Theorem and a 
second result which we call the Fixed Point Theorem. Both results involve 
fixed points of the renormalization operator. They are the cleanest applica- 
tions of the Renormalization Theorem. In the next chapter, we will give some 
messier applications, including the proof of the Near Reduction Theorem. 

We first discuss the general nature of our calculations. We use the nota- 
tion from the previous chapter. Let A be the set of A-chains. Likewise define 
B. The Renormalization Theorem gives us a canonical 3-to-l map from A to 
B. This map does not depend on any conventions concerning the boundaries 
of the layers of A. Let x : ^ — > B denote the map just described. 

Since S can be considered as a dynamically invariant subset of E, we can 
transfer the structure to E. For instance, an A- atom of S is a connected 
intersection of E with the A-atoms. Every polygon in sight is a golden 
polygon. We deal with these polygons and polyhedra using the methods 
from §6. 

Given an A-atom X and a S-atom Y, we write X — )■ y if = \Y]. 

Here [X] is the forward A-chain containing X and [Y] is the B-chain con- 
taining Y. We take the 5-chains to be forwards when they are discussed in 
connection with type-1 A-chains, and backwards when discussed in connec- 
tion with type-2 74-chains. 

Important Observation: By construction, a generic point never intersects 
the boundary of an 74-atom or the boundary of a 5-atom. 

11.2 Proof of the Fundamental Orbit Theorem 

Referring to the notation in the Fundamental Orbit Theorem, we find a 
partition of each by a union of yl-atoms and periodic tiles such that 

r^r'(4(r)), (103) 

for every A-atom r. The periodic tiles all lie in T(2). 

For future reference, we refer to the above mentioned partitions as the 
A- partitions. The left side of Figure 11.1 shows the partition for T^^. 



110 



Figure 11.1: The partition of TJJ and an edge-covering. 



The y4-partition of T^^ is a union of 41 = 34 + 7 A-atoms and periodic 
tiles. The periodic tiles are all shaded. The dark polygon on the side is part 
of the Penrose kite. The horizontal lies through the periodic tiles correspond 
to the tops and bottoms of the layers of A. Thus, for instance, even though 
the large octagon is a single periodic tile, we treat it as 4 separate tiles in 
our proof. Each of the 3 triangles T^, and T^i are also partitioned into 
41 = 43 + 7 polygons. These partitions look exactly like the one in Figure 
11.1, up to reflection in a vertical line. The triangle T^^ is partitioned into 
4 = 3 + 1 polygons and T^i is partitioned into 3 = 2 + 1 polygons. It turns 
out that we do not need to consider the 4 remaining triangles, as we explain 
below. The reader can see all the partitions, and in color, using our applet. 
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Lemma 11.1 Suppose that Equation \10S\ holds for all the A-atoms in the 
A-partitions . Then the Fundamental Orbit Theorem is true. 

Proof: Let Oi be a generic orbit that intersects T^ but avoids set T(2). 
Suppose first that Oi is not a distinguished orbit. Since all the periodic tiles 
in the A-partitions lie in T(2), some point pi G Oi intersects an 74-tile r 
of the y4-partition. Since pi is generic, it cannot intersect the boundary of 
an y4-atom. Hence pi intersects the interior of some A-atom r of the A- 
partition. But then, by Equation 11031 we have Oi -w O2, where O2 is the 
\E'-orbit containing p2 = ^lj~^{Tj^{pi). 

When Oi is a distinguished orbit, we attach it to one of the slabs or the 
other. The argument above is the same, except that what we get is that pi 
must lie in the interior of one of the horizontal edges of the A-partition. But 
then Equation 11031 gives us the same result as for the ordinary orbits. Were 
we to attach the special orbit to the other slab, we could still use Equation 
11031 but it would involve a different A-atom. 4 

To prove the Fundamental Orbit Theorem, it only remains to check Equa- 
tion [TD3] for each of the A-atoms in the A-partitions. Each individual check 
is a finite calculation, similar to what we did in the last chapter. We now 
discuss this calculation. 

First of all, we reduce the calculation from 10 triangles to 6 triangles, as 
follows. We have the relations 

vl/(Tfi) = T+; = TA; ^iT,\) = T3+; ^(T.^) = T3+. (104) 

Using these relations, we see right away that the truth of the Fundamental 
Orbit Theorem for each triangle on the right hand side of the relation implies 
the truth of the Fundamental Orbit Theorem for the triangle on the left hand 
side. Thus, we only need to check the 5 triangles T^j and the triangle T^^. 

First of all, we check that each relevant triangle T^ is indeed partitioned 
by the polygons we have listed. The check is similar to what we did in 
§6. As for verifying Equation 1103^ we have to make the same verification 
141 = 4 X 34 + 3 + 2 times. 

Let P and Q respectively be the polygons on the left and right hand sides 
of Equation 11031 To check that P — t- Q, in each of the 141 cases, we do the 
following calculation. 
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• We check that P = 0(-P) is contained in a A- atom. This amounts to 
producing the two integers m and n, with m < < n, such that is 
entirely defined on the interior of P for all k G [m, n ], and ^''{P) C A 
for k = m,n but not for k G {m,n). We verify these things using the 
methods discussed in §7.5[ Define Pq = \&'"(P). 

• We check that Q = &{Q) is contained in a 5- atom. This is the same 
kind of check as in the previous case. Let Qo = "^"^{Q) be the iterate 
that is contained in B. (Depending on the layer, m is either non- 
positive or non- negative.) 

• We check that R{Po) = Qq. 

We simply perform the 141 calculations and observe that they all work. 
This completes the proof of the Fundamental Orbit Theorem. 

Remark: Note that we are not quite verifying that P and Q are atoms, 
but only that they are contained in atoms. Similarly, we are only verifying 
that P and Q are contained in atoms, and that their boundaries are all con- 
tained in lines defined over Z[(j)]. This suffices for our purposes. We mention, 
however, that we produce our A-partitions by slicing the A-atoms and peri- 
odic tiles by the image of S. Thus, were to make the extra check, it would 
work out. 

The reader can use our applet to survey the calculations. Here we give an 
example illustrating the rough size and complexity of the calculation. There 
is a small isosceles A-tile in the partition of T^i that contains the rightmost 
vertex of T^. This triangle is perhaps too small to see in Figure 11.1. The 
vertices of this triangle are 

(997-6160,-236 + 1460) (1984-12260,-3 + 20) (10-60,-3 + 20). 
In this case 

Po = §-^°«(x); go = ^'\Q). 

Things don't get much worse than this, because the maximum chain length 
is 703. 
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11.3 The Fixed Point Theorem 



It turns out that there is a second region that plays a role that is similar to 
the fundamental triangles. Let Qq be the quadrilateral with vertices 

(6-20,0-=^); (3,0); (20,0-'); (3,20-^) (105) 

Qo is a small kite whose line of symmetry is the vertical line x = 3. The bot- 
tom point of Q is (3, 0). See Figure 11.2 below. Let D denote the homothety 
that fixes (3, 0) and expands distances by a factor of 0^. Define 

Qn = D-"iQ). (106) 

The family {Qn} is a shrinking family of kites that limits to (3,0). Let p 
denote the reflection in the vertical line x = 3. 



Theorem 11.2 (Fixed Point) Let Oi be any sufficiently long and generic 
orbit that intersects Qi. Then Oi -w O2, and O2 is such that 

V''(02)ngo = poD(Oingi). 

Here O2 is the reflection of O2 in the x-axis andip' is the outer billiards map. 



In particular Oi — ?• O2, in the sense of Theorem\6.4 



Proof: The proof here is essentially the same as the proof we gave for the 
Fundamental Orbit Theorem. Let A : — t- be the dilation by 0^ such 
that A(3,0) = (—3, —2). We have the equation 

^' o A(p) = p o D{p) (107) 

for all p E Qi. 

We exhibit a tiling of Qi by a union of 48 polygons, 42 of which are A- 
atoms and 6 of which are periodic tiles. For each A-tile r in the partition, 
we verify, by the same means as above, the equation 

r^A(r), (108) 

Our result follows from Equation I107[ Equation I108[ and the same argument 
as in Lemma 111.11 ^ 
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Figure 11.2: Qi and its partition into atoms and periodic tiles. 



It is a consequence of the Fixed Point Theorem that the point (3, 0) is 
a point of dilation symmetry for the dynamical tiling associated to outer 
billiards on the penrose kite. Figure 1.3 shows what is going on. The white 
regions in Figure 1.3 are small copies of portions the tiling T and the shaded 
regions are increasingly small periodic tiles. The reader can probably spot 
Qi in Figure 1.3: It is flanked by the two largest dark tiles. 
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12 Proof of the Near Reduction Theorem 



12.1 Discussion 

The main structural relation between the y4-atoms and the S-atoms is the 
3-to-l map from the set of A-chains to the set of i?-chains. Now we mention a 
second piece of structure. Though we did not formally check this in general, 
in the cases of interest it turns out that each 5-atom is partitioned into a 
finite union periodic tiles and A-atoms. Thus, we have a kind of auxilliary 
dynamical system defined on the level of A-atoms. Starting with a A-atom 
T, we can consider all those S-atoms - there are finitely many, and all in 
the same chain - such that r — )■ r'. Next, we subdivide each such r' back 
into y4-atoms. And so on. All our proofs in this chapter amount to tracing 
through this auxilliary dynamical system. 

Our computer program is designed so that the user can trace through all 
the calculations manually, plotting the relevant tiles. 

The way to view this chapter is that we have one main calculation, con- 
tained in §12. 4^ that takes care of almost all the orbits that come fairly close 
to the origin. The remaining calculations, such as those in §12. 3i deal with 
the several small regions that the big calculation does not cover. 

12.2 Points Very Near the Kite 

In our next result, the notation Oi O2 is as in Theorem 16.41 Let 

= [0,r] X [0,2]. (109) 

Lemma 12.1 Let r = 2 + Let Oi be any generic infinite orbit that 

intersects Sr- either Oi intersects or Oi — ?■ O2 and O2 intersects T"*". 

Proof: Our proof refers to Figure 12.1. Let r = 2 + Let p denote 

complex conjugation - i.e. reflection in the Figure 12.2 shows a 

covering of Sr by 13 tiles. The left black tile is the Penrose kite, and the 
right vertex of the rightmost dark tile lies on the vertical line y = r. The 
two white tiles comprise T, the fundamental triangle. Let Xi be the lower 
of the two black tiles and let L be the smallest light tile - namely, the little 
kite. Note that Xi and L share an edge and Xi U L is a triangle. 
Short calculations reveal the following. 
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• \1/ is the identity on all the light tiles. 

• The top black tile X2 is such that (V^')"^ ° ^^(^2) C T. 

• The leftmost dark tile Di is such that ^(-Di) C T. 

• The bottom dark tile D2 is such that p o ^[D2) = Xi. 

• The rightmost dark tile is such that p o ^(Da) = XiU L. 




Figure 12.1: A covering of S,.. 



The region Xi is the triangle T^i analyzed in the proof of the Fundamental 
Orbit Theorem. By the Fundamental Orbit Theorem, the orbits intersecting 
Xi renormalize to orbits that intersect T. Lemma [12.11 now follows directly 
from case-by-deductions. We analyze of the cases in detail. The remaining 
cases are similar, and in fact easier. (We chose the most elaborate case.) 
Suppose that Oi intersects D3, the rightmost dark tile. Then an associate 
0[ intersects LU Xi. Being an infinite orbit, 0[ does not intersect the inte- 
rior of L. Being a well-defined orbit, 0[ does not intersect dL U dXi. Hence 
0[ intersects the interior of Xi. But then 0[ ~^ O2, where O2 intersects T. 
This shows that Oi — > O2 where O2 intersects T. ^ft 
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12.3 Two Trouble Spots 



Let a be the triangle with vertices 

(-15 + 110,2-0) (-36 + 240,-19 + 120) (53-310,2-0). (110) 

Let a* be the triangle with vertices 

(21-110,2-0) (42-240,-19 + 120) (-47 + 310,2-0). (Ill) 

a and r are two small A-atoms contained in Qq — Qi. Reflection in the 
vertical line x = 3 interchanges them. Let ip' denote the outer billiards map. 

Lemma 12.2 Let Oi be any generic infinite orbit that intersects a. Then 
we have Oi O2 ~^ O3 and iIj'{0^) intersects T~ . 

Proof: Let di = a. Using the notation above, and the same kind of calcu- 
lations, we check that cri -w erg, where (J2 is the triangle (a translate of ui) 
with vertices 

(-11 + 90,-2 + 0) (-32 + 220,-23 + 140) (57-330,-2 + 0). 

We check that cti is a union of one periodic tile and 2 A-atoms, o"2 and a^. 
Next, we check that o'2 — )■ and as — )■ cr5, where aQ = a^U is the triangle 
with vertices 

(5-40,-6 + 40) (18- 120,7-40) (-3 + 0,2-0). 

A direct calculation shows that ip'{(y&) C . An argument just like the one 
given in Lemma [11.11 finishes the proof. 4ft 

Lemma 12.3 Let Oi be any generic infinite orbit that intersects a*. Then 
we have Oi O2 O3 O4 and ip^O^) intersects T^ . 

Proof: The argument here is very similar to what we did in the previous 
lemma. We omit the details. ^ 
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12.4 The Big Calculation 



Now we prove a result similar to the ones in the previous section, except that 
the computational part of the proof is much more intensive. 

Given any p = {x, y), define \p\x = \x\. In our next result, the significance 
of the number 2(j)~^ is that it is precisely the horizontal width of Qq. 

Lemma 12.4 Let r = 2 + 0^^. Let pi G — Sj. — Qi — cr — a* be a point 
with a generic infinite orbit Oi. Then Oi O2 O3, and there is some 
A; G {1, 2, 3} and some qu G Ok such that \qk\x < bfcU ~ 20~^. 

Proof: Say that a symmetric return tile is a maximal polygon in the main 
strip S on which both ^ and are defined. For each point p in a symmetric 
return tile r, the expressions 

are independent of the choice of p. If /(r) = it means that ^ is the identity 
on r. Otherwise, by Equation[3Sl we have min |i^'-|-(r)| < 2(j)^^ . Accordingly, 
we only have to consider points in symmetric return tiles r such that Kj^{t) 
and K_{t) are both positive. Call such tiles positive symmetric return tiles, 
or PSRT's for short. 




Figure 12.2: Four of the PSRT's. 
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It turns out that there are 23 PSRT's that intersect 5*24. The PSRT 
closest to the origin is Qq. Figure 12.2 shows the first 4 PSRT's, shaded, as 
well as a portion of the partition of the top half of S into symmetric return 
tiles. (The righmost shaded piece is a union of two tiles.) The reader can see 
the remaining tiles using our applet. 

We partition the union of 23 PSRT's by 848 polygons, each of which is 
either an A-atom or a periodic tile. Since we are only interested in infinite 
orbits, we ignore the periodic tiles. Let r be an A-aiom. in our partition, and 
let [r] be the chain containing r. Say that r is slack if there is some other 
yl-atom r' C [r] such that 

max \p'\n, < max \p\^ - 20"^. (112) 

p'er' per 

Otherwise, we call r taut. A taut atom is one that is nearly as close as 
possible to the origin within its chain. Note that r' need not be in our 
partition. 

Ignoring the atoms in Qi, which we do not need to consider for this 
lemma, we find that there are 123 taut atoms. With 12 exceptions, we take 
each taut atom r and exhibit a 5-atom r' such that r — )■ r' and the pair 
(r, r') satisfies Equation 11121 Call the 12 exceptional atoms super taut. The 
12 super taut atoms cannot be brought substantially closer to the origin by 
a single renormalization. 

To deal with the 12 super taut atoms, we trace through the second renor- 
malization, as follows. Let r be a super taut atom. 

• We find a i?-atom r' such that t ^ t' . In practice, we take r' to be 
the atom in its chain which is closest to the origin. 

• The atom r' is partitioned into a finite union r^, of A- atoms and 
periodic tiles. The value k depends on r. The largest value is 27. 

• For each periodic tile we exhibit a i?-atom such that t[ — )■ r" and 
(r", r) satisfies Equation 11121 Again, we take r^' so that it minimizes 
the distance to the origin within its chain. 

The analysis above handles the points in 10 of the 12 super-taut atoms. 

The two exceptional super taut atoms are a and a*, which we have de- 
liberately excluded in our hypothesis. 4 
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12.5 Proof of the Near Reduction Theorem 



The Near Reduction Theorem concerns generic infinite orbits that intersect 
the rectangle S24 = [—24,24] x [—2,2], but our results above concern the 
smaller region 5*24 = [0, 24] x [0, 2]. However, replacing an orbit Oi by on of 
its 3 associate orbits, we can assume that the orbit intersects 6*24 . Let Oi be 
such an orbit. 

Since Oi is a generic infinite orbit, we have Oi O2 O3.... These 
orbits intersect the strip E = x [—2, 2] but they might not intersect the 
smaller region S = x [0, 2]. However, replacing our orbits with associates, 
if necessary, we get a new sequence Oi — )■ O2 — > O3... in which each orbit 
intersects S. 

Remark: We allow ourselves the liberty of passing to associates mainly 
for convenience. We might have proved the Near Reduction Theorem for the 
original sequences of renormalizations, but it is more computational work. 

Say that an orbit O is passed Qi if there is some p E O such that 

IpI^ < min \q\^. 
q&Qi 

Lemma 12.5 The Near Reduction Theorem holds for a generic orbit Oi that 
is passed Qi. 

Proof: Iterating Lemma 112. 4[ we find that there is some m such that one 
of three things is true. 

1. Om intersects S'^ for r = 2 + 0^^. 

2. Om intersects cr U a*, the union of two special A-atoms. 

In either case, the lemmas in the preceding section combine to show that 
Om+k intersects T+ for some k G {0, 1, 2, 3}. ^ 

Now let's consider the general case when Oi is not necessarily past Qi. 
We can still iterate Lemma 112. 4[ but now we arrive at the possibility that 
Om intersects Qi for some m. The map \E' is the identity on the line y = 0, 
so Om does not intersect this line. Hence, there is some e > such that Om 
avoids the e neighborhood of the point (3, 0). But then there is some n such 
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that Om intersects Qn but not Qn+i- Applying the Fixed Point Theorem, 
we see that O^+i intersects Qn~i but not Qn, and Om+2 intersects Qn-2 but 
not Qn-i- Continuing in this way, we see that Oi — )■ ... — )■ Om — > ••• — ?■ Oh, 
where Oh intersects Qq but not Qi. 

If Oh intersects 5*^ U a U cr* then we have the same situation as in Lemma 
112.51 Otherwise, we can apply Lemma [12.41 one last time to produce an orbit 
Oh+k which is passed Qi. This works thanks to the following three facts. 

• Oh intersects Qq in an interior point (because the orbit is well-defined). 

• \Oh+k\x < \Oh\x - 2(f)~^. 

• Qq has width 20~^. 

We finish the proof by applying Lemma [12.51 to Oh+k- 

For each orbit Oi, the Near Reduction Theorem produces some integer 
m = m(Oi) such that m such that Oi — )■ ... — )■ Om, and Om intersects T+. 
There is no uniform m which works for all orbits. The difficulty is that orbits 
that come sufficiently near (3, 0) do not intersect T"*", and iterated renormal- 
ization only gradually moves such orbits away from (3,0), as a consequence 
of the Fixed Point Theorem. On the other hand, our analysis above yields 
the second statement of the Near Reduction Theorem: 

Lemma 12.6 Let y G (0, 2) be any point such that {/?"(?/)} does not contain 
0. Then there is some N with the following property. If Oi contains a point 
within 1/N of the segment [—24,24] x {y}, then Oi — i- ... — > Om OLnd Om 
intersects and m < N. 

Proof: Let's consider iteratively applying Lemma [12.41 to Oi. After we ap- 
ply Lemma [12.41 at most 24/0~^ times, the resulting orbit is either passed Qi 
or else intersects Qi. In the latter case, there is some e such that every point 
of the resulting orbit is at least e from (3, 0). The number e works uniformly 
as long as Oi contains a point sufficiently close to our line segment. But 
then, an additional C log(e) renormalizations produce an orbit that intesects 
T"*". Here C is some uniform constant. ^ 
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13 Proof of the Far Reduction Theorem 



13.1 Deviation Estimates 

Let G : S — 7- S be as in the Compactification Theorem. Recall that an A- core 
is a finite portion a = {po, ...,Pn} of a \E'-orbit with the following properties. 

1. &{po) G A, the first of the renormalization sets. 

2. 6(pi), Q{pn) do not lie in A. 

3. e{pn+i) e A 

Up to translation, there are only finitely many A-cores, one per A-chain. To 
the y4-core a we associate the deviation interval 

{a) = min7ri(pfe -po),max7ri(pfe -Po) . (113) 

. k k 

Here tti is projection onto the first coordinate. That is, we translate the 
A-coie so that its initial point is at the origin, and we measure how far the 
sequence deviates to the left and right of the origin. 

The sequence of points {pk — Po} associated to a is combinatorial in 
nature: In view of Equation[3Sl there is a finite sequence {(ofc, &fc)} such that 

7ri(pfc -Po) = 2v4a„ + 26fc; k = l,...,n. (114) 

Indeed, we can form a connected lattice path by connecting the consecutive 
integers (0,0) = (oq, feo), (^n, &n)- We call this path an A-strand. An A- 
strand is nothing more than an arc of the arithmetic graph, translated to 
that the initial point is the origin. The A-strands are not essential for our 
analysis, but we introduce them because they help us visualize what is going 
on. 

For the purposes of drawing pictures, we will continue the strands to 
the point (a„+i,6„+i), so that the endpoints of the strand correspond to 
points in A. In this way, the arithmetic graph component corresponding to 
an infinite orbit is partitioned into uniformly short arcs, with the endpoints 
being special points that map into A. The distribution of these special marker 
points controls the coarse geometry of the arithmetic graph, and thereby the 
coarse geometry of the orbit itself. 
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We define -B-cores and i?-strands in the same way, except that in the 
type-2 case, the last point of the B-core maps into B. The Renormahzation 
Theorem gives us a 3-to-l- map from the set of A-corcs to the set of 5-cores. 
As in the Far Reduction Theorem, the notation a /3 indicates that the 
B-core (3 is the image of the A-core a under this map. 

Let Nk{S) denote the /c-tubular neighborhood of a set S. Let (f)~^I denote 
the interval obtained by scaling I down by a factor of (f)~^ about the origin. 

Lemma 13.1 (Trivial) There is a universal constant U with the following 
property. Suppose that a ^ (3. 

1. If a has type 1 then {a) C Nu{{l3)). 

2. If a has type 2 then (t)~^{a) C Nu{{/3)). 

3. If a has type 1 then {(3) C Nu{{a)). 

4. If a has type 2 then C Nu{4>~^{a)). 

Proof: Up to translation there are only finitely many cores. Hence, all de- 
viation intervals are contained in a single compact subset of R. 4^ 

Remark: There is nothing special about (f)~^ in the Trivial Lemma. Were 
we to pick any other constant, we would get a similar result. However, we 
place in this result for the sake of comparison with the more effective 
result that we establish below. For the more effective result, the choice of 
is the natural choice. 

Items 3 and 4 serve our needs perfectly. They correspond to the constant 
C in the Far Reduction Theorem that we don't care about. However, we do 
need a more effective version of Items 1 and 2. 

Lemma 13.2 (Deviation) Suppose that a ^ /3. 

1. If a has type 1 then (a) C Nr{{p)). Here r = -20 + 24(/) < 12 - 20. 

2. If a has type 2 then (p-^{a) C N5{{l3)). 
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Proof: We simply enumerate all the integer sequences and check that the 
bounds hold. We use golden arithmetic to make sure that everything is com- 
puted correctly. ^ 

Reiricirk: We could get similar improvements on Items 3 and 4 of the Trivial 
Lemma, but we don't care about those cases. 



13.2 An Estimate for the Markers 

The Trivial Lemma and the Deviation Lemma reduce the Far Reduction 
Theorem to a statement that just concerns points that map into A and B, as 
we now explain. Say that an A-marker is a point p G E such that 0(p) G A. 
We define the type of the point just as we defined the type of an A-coie: It 
depends on which layer of A contains the image. 

We write a 6 if a is an A-marker and 6 is a S-marker, and 

RoQ{a)=b. (115) 

By definition a ^ 6 if and only ii a /3, where a is the A-chain containing 
a and (3 is the S-chain containing b. Below we will prove the following result. 

Lemma 13.3 (Marker) Suppose that a b. 

1. If a has type 1 then \b\x < \a\x + 24>. 

2. If a has type 2 then \b\x < (t)~^\a\x + 10. 

The lower bounds in the Far Reduction Theorem are immediate conse- 
quences of the triangle inequality, the Marker Lemma, and Items 3 and 4 of 
the Deviation Lemma. Let's establish the upper bound \P\x < \ol\x + 10 in 
the type-1 case. Let a and b respectively be the markers in a and ^. Choose 
the point a' G o; that realizes \a\x. Set r = 12 — 20. We have 

\a'\x-\a\xe{a)^Nrm). (116) 

Therefore, there is some 6' G /3 such that 

\b'\x - \b\x < W\x - \a\x + (12 - 20). (117) 
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By the Marker Lemma, 

l&'U- |aU < I6U + 20. (118) 

Therefore 

\b'\x — \o.\x < \o!\x — + 12. (119) 
Adding \a\x to both sides, we get 

\P\x<\h'\x<Wl + l2=\a\^ + l2, (120) 

as desired. The proof of the upper bound in the type-2 case works the same 
way. 

The remainder of the chapter is devoted to proving the Marker Lemma. 
The type-1 case turns out to be fairly trivial, and the type-2 case has an 
arithmetical feel to it. 

13.3 The Type-1 Case 

Let us recall the information given in S|5l First of all, we have 

0(^,1/) = (l,^,0) + Q,^,y) (121) 

The image of is a certain fundamental domain of the form F x [0, 2], where 
F is described in §5.21 For our purposes, we can deal with the simpler region 

VL = [1/2, 3/2] X [0, 2] X [0, 2] C F X [0, 2] (122) 

because (by direct inspection) this region contains both B and A. 

In the type 1 case, namely in the slabs corresponding to the partition 
intervals I2 and I 4^, the formula for R is 

R{x,y,z) = {x,y,z±2(j)-'). (123) 

In this case, we can restrict i? to a smaller domain S C ^ such that R{S) C f2 
as well. 

Given the explicit equation for R, we get the following result. If a = {x,y) 
and a ^ b, then 

b= {xT2(p,y±2(p-^ ±2). (124) 

What makes this equation work is that 20 — 2(j)~^ = mod 2Z. Once we see 
that this formula works, it must be the uniquely correct formula, because 
is an inject ive map. 

Equation 11241 immediatly implies the type-1 case of the Marker Lemma. 
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13.4 The Type-2 Case 



The type-2 case involves the 7 maps Rijk, for various choices oi i,j,k. These 
maps are hsted in §5.41 Define 

S,,, = R7^,in). (125) 

Here 5',^^ is a small rectangular solid that either contains or abuts the fixed 
point k) of Rijk- In the Renormalization Theorem, the map R^k is only 
applied to polyhedra that are entirely contained in Sijk- 

Since Rijk expands distances by a factor of 0^, we see that the sidelength 
of Sijk in the x-direction is and the sidelengths in the other directions 
are 20~^. 

For the sake of notational convenience, we will give our proof in the case 
of -Rill. The other cases have essentially the same proof. 

Let p = (pi,P2) be some point. Let Qk{p) denote the kth coordinate of 
G(p). We first observe that Qi{p) = 1 if and only if pi = 2N(j) for some 
integer N . So, we can always write 

pi = 2iVp0 + ep, (126) 

where ep is an error term of size at most 0. 

When p = a, we have a better estimate. Since 0i(a) lies within 0~'^/2 
of 1, the error term is at most 0~^/2. We also observe that 02 either lies 
within (f)~^ of 1 or within of —1. We will treat the former case. The 
latter case has the same treatment. Summarizing the discussion, we have 

a={2Na(P + ea,l + 5a); le^l < ^"'A |5a| < 0"'- (127) 
From the equation 

Gkib) - 1 = (Qk{a) - 1) X (128) 

we get 

6= (2Arfe0 + 03e„l + 035„). (129) 

In the formula for G, we take the listed coordinates and suitably translate 
them by even integers until these coordinates all lie in [0, 2]. Hence, there is 
an even integer such that 

e2(a) = iV,0 + + ej2 - 6J2 e [0, 2]. 
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(130) 
(131) 



(133) 



We write = 1 + M„, so that 

G2(a) = l + Na<p + Ma + e,/2 - 5a 1 2 

Similarly, 

62(6) = l + N,<P + M, + 03(efe/2 - 5,/ 2). 
Given our bounds on ea and (5a, we have 

Mae[-Na<p-2,-Na<p + 2]. (132) 

Combining the k = 2 case of Equation 11281 with Equations 11301 and I13H 

we get 

Expanding out the right hand side of this last equation, we get 

(2 + 3(j))Na + (1 + 20)Mb = 
{3Na + 2M,)<f)+{2Na + N,). 
Equating coefficients, we get 

Nh = 3Na + 2M„. 
Combining this last equation with Equation I132[ we find that 

G [(3 - 2<j))Na - 2, (3 - 2<j))Na + 2] = [<f)-^Na - 2, <j)~^Na + 2] 

In short, 

0''|iVa| < \Nb\+2. 

Multiplying through by 2(f), we get 

<p-^x (20|iV,|) <20|iVfc|+40. 
Finally, we mention that 



(134) 

(135) 



(136) 
(137) 



\aL-2(h\Na 



\bU-2<f)\N,\ =e,<<f)/2; 



By the triangle inequality, we get 

bx < (p'^a^ + 40 + 0/2 + 0~^ X (0"V2) < 10. 
This completes the proof of the type-2 case of the Marker Lemma. 



(138) 



(139) 
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14 Coordinates 



14.1 The Polyhedron Exchange Map 
14.1.1 Notation and Conventions 

The polyhedron exchange map is defined in terms of 64 polyhedra Po, -Pes- 
The polyhedra Pq, P21 are what we call inactive, in the sense that \1/ acts 
as the identity on these polyhedra. The remaining polyhedra are what we 
call active. The 42 active polyhedra are such that 



Pj+21 



22, ...,42. 



(140) 



Here / is the involution 



Iix,y,z) = (2,2,2) -{x,y,z). 



(141) 



The map I commutes with Thus \I' has the same action both the lefthand 
and the righthand polyhedra listed in Equation I140[ Accordingly, to save 
space, we will just list Pi, P42. 
Our notation is such that 



ao ai 
bo hi 

Co Ci 



(ao + ai0, 60 + Co + ci<; 



(142) 



We list each polyhedron by its vertices. The vertices are not given in any 
special order. 

14.1.2 The Inactive Polyhedra 



PO : 



3 -1 
1 
2 -1 





-3 


3 







1 




1 







3 


-1 






2 







2 







1 

















2 


-1 




2 







2 







2 








1 

2 



"-3 


CO 




"-1 


to 




"0 


1" 


to 


-1 




2 







2 





2 







2 














PI : 



-1 1 

2 
2 -1 





5 -3 




2 


-1 




1 







-1 1 






1 












1 







2 






2 -1 




2 







2 







2 





5 -3 
1 
2 



3 -2 

2 



to 


-1 


to 





to 





to 


-1 


to 
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5 


-3 




2 


-1 




2 


-1 




1 

































1 












1 




2 
























-1 1 

2 




"2 


-l' 




"5 


-3' 




CO 


-2" 


2 










1 

































P4 : 



"l 


o' 












1 










"l 


o' 




to 


-l' 




"l 


o" 




















1 






























P5 ■- 



"l 


o' 




to 







1 










"l 


o' 




"o 


1" 




"l 


o" 




2 







2 







1 







2 







2 







2 






P6 : 



"l 


o' 




1 







1 










"l 


o' 




"o 


1" 




"l 


o" 




2 







2 







1 






























P7 - 



"l 


o' 




1 







1 








6 -3 
-1 1 
2 





"l 


o' 




to 


-l' 




"l 


o" 




















1 







2 







2 







2 








"l 


o' 




"l 


o' 






1 







2 









1 







1 








CO 






-4 


3 " 


2 






3 


-1 


-1 






3 


-1 


"l 


o' 




"l o" 


1 







2 


2 











-4 3 
4 -2 
-1 1 



-3 
1 
1 





1 







1 

















1 









1 







1 








-3' 






-4 




3 " 


1 














-1 






3 




-1 


'l 


o' 




"l 


o" 










1 





2 
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" 6 


-3 




1 







6 


-3 




PIO = 




























-1 


1 






















"-4 


3 




1 







-4 


3 




Pll = 


3 


-1 




2 







2 









-1 


1 





















11 


-6' 




" 6 


-3 


-3 


2 




-1 


1 










_ 





-9 


6 " 




"-4 


3 


5 


-2 




3 


-1 

















P12 : 



P13 : 



P14 : 



P15 : 



11 


-6 








-3 


2 


-9 


6 


5 


-2 


-3 


2 


" 5 










-3 


2 





"6 -3' 




















"-4 3' 






2 














"o 


l' 

























2 

2 




11 


-6 


-3 


2 








-9 


6 


5 


-2 








-3 


3 " 


2 


-1 





_ 


5 


-3 





1 











"3 -l" 









_ 




"-1 l" 




2 








-3 2 
2 




P16 : 



P17 : 



"l 


o' 




"6 


-3' 




















2 







2 







"l 


o' 




"-4 


3' 




2 







2 







to 







2 








11 

-3 
2 



6 


-3 




6 


-3 


-1 


1 




-1 


1 


2 







CO 


-1 



P18 : 



P19 : 



P20 : 



"6 -3' 









2 




"-4 3" 




2 




2 





"o 


l' 












2 








11 


-6 


-3 


2 


2 





-9 


6 


5 


-2 


2 





-3 


3 " 


2 


-1 


2 








"3 -l' 













2 






"-1 l' 






2 






2 





11 


-6 


-3 


2 


5 


-2 


-9 


6 


2 





5 


-2 



5 -2 
-3 2 
5 -2 



P21 : 



2 -: 

2 
2 





5 


-3 




-3 2 









1 




2 






2 







2 





-3 2 
2 
5 -2 
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14.1.3 Half the Active Polyhedra 



P22 : 





1 







1 









1 







2 









1 
















P23 : 





3 


-2 




3 


-2 




2 


-1 




2 


-1 






























2 









1 







2 







2 
















5 -3 
1 






3 -1 




P24 = 


1 






2 -1 





-5 4 

5 -2 





3 


-1 







1 
























2 







2 













1 




2 
















-1 


1 




1 







1 







2 


-1 




2 


-1 




P25 = 










1 

































2 


-1 




1 







2 







2 
















-1 1 





3 -1 




P26 = 









2 -1 





11 

-3 
5 





3 


-1 




6 


-3 
























2 







2 








11 -( 

-3 2 






6 


-3 





















8 -5 




3 


-2 




3 


-2 




8 -5 


P27 = 























-1 1 




-1 1 




1 



















P28 : 



5 -2 
-3 2 
5 -2 





5 


-2 







1 
























2 







2 








1 

3 -1 



P29 : 



1 

5 -2 





1 







6 


-3 
























2 







2 








6 


-3 


-1 


1 


CO 


-1 



P30 : 



-9 
7 

-3 





■-1 


1 " 




■-4 


3 ■ 




"-9 


6 ■ 




"-1 


1 






























2 


-1 






2 


-1 




4 


-2 




4 


-2 
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P31 : 



5 


-3 



-3 3 
2 -1 
2 -1 








1 




5 


-2 
























4 


-2 




4 


-2 










l" 



















P32 : 



-7 5 
6 -3 
1 





3 


-1 




1 

















1 












1 




5 


-2 





-5 4 

5 -2 





3 


-l" 












2 






"l 


o' 




1 
















P31 



-3 3 
2 -1 
2 -1 








1 




5 


-2 
























4 


-2 




4 


-2 










l" 



















P32 : 





3 


-1 




1 

















1 












1 




5 


-2 







3 


-l" 












2 






'l 


o" 




1 
















-5 4 
5 -2 




P33 : 



"l 


■ 




"l 


o' 




1 







1 







5 


-2 




2 








P34 : 



-9 6 



-4 3 
4 -2 
-1 1 



-4 3 

4 -2 





"l o" 













1 






"1 0" 






1 













3 

-4 
4 






■5 


-2 







1 







1 




P35 = 




























4 


-2 




4 


-2 




3 


-1 







5 


-2 
















2 








-3 3 
2 -1 
2 -1 



P36 : 





"-1 


l' 




"-1 


1" 




2 







2 







2 










1 
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1 


o' 




P37 = 


2 









-3 


2 





-4 3 
3 -1 
-1 1 





"l 


■ 






2 









4 


-2 





-9 6 
5 -2 
1 



-9 


6 


5 


-2 









-4 3 
3 -1 
3 -1 

"-4 3" 
2 




P38 : 



-9 6 
5 -2 
2 



-9 6 
5 -2 
1 



P39 : 



-4 3 
3 -1 
3 -1 





"l 


o' 




"l 


o" 




1 







2 







2 







1 








"-11 8 ■ 




"-6 


5 ■ 




"-6 5 ■ 




"-3 


3 ■ 




"-11 8 ■ 




"-3 


3 " 


P40 = 


2 




2 







3 -1 




2 


-1 




4 -2 




2 







4 -2 




4 


-2 




3 -1 




2 







2 




2 


-1 




P42 : 



-4 3 
3 -1 
3 -1 



1 

2 
2 



1 





2 





4 


-2 


-4 


3 


2 





4 


-2 



Using our computer program, the reader can survey all these polyhedra. 
The program is set up so that the vertices of the polyhedron are displayed 
whenever the polyhedron is selected. 
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14.2 The Map 



Now we describe the action of the map \& on the polyhedron Pj. To each 
integer j we assign a pair of integers Nj according to the following lookup 
table. 



- 


-21 


(0,0) 


22 


-29 


(1,0) 


30 


-32 


(1,1) 


33 


-40 


(0,1) 


41 


-42 


(-1,1) 


43 


-50 


(-1,0) 


51 


-53 


(-1,-1) 


54 


-61 


(0,-1) 


62 


- 63 


(1,-1) 



(143) 



For instance, if j = 41 then Nj = (—1,1). The numbers Nj = {m,n) are 
such that 



^(p) -p= (2m + 2ra0^^*) (144) 

for any point p such that 0(p) G Pj. As in Equation |3S1 the third coordinate, 
which lies in {—2, 0, 2} depends on the parity of m + n. When m + n = 
the third coordinate in Equation 11441 is 0. Otherwise, it is ±2, depending on 
which value yields a point in the strip E. 

as follows. 



The action of \E' on Pj is expressed in terms of the pair 



m, n 



P 



—2m + lOn 




2m — 'on 
In 




(145) 



Thus, for instance, for iV4i, the pair ( — 1,1) yields the vector 

(12-80,-2 + 20,0). 

Notice that A^2i+j = foi' 3 = 22, ...,41. This is consistent with the 
action of the involution / defined in the previous section. 

Now we discuss the allowable transitions for our map. We write a — )■ 6 is 
there exists a point p in the interior of Pa such that \l/(p) G Pt- For a < 22 
we only have a — )■ a. We also have the general symmetry 



a ^ b 



(a + 21) (b + 21] 



(146) 
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which holds for all a > 22. For this reason, we just hst the transitions for 
a = 22, ...,42. 

22 61 23 56 

23 63 58 29 32 26 24 

24 ^ 25 53 37 

25 ^ 31 35 28 27 23 62 61 

26 ^ 25 

27 ^ 63 

28 25 

29 ^ 23 

30 -> 23 

31 ^ 32 

32 ^ 25 30 47 36 34 37 42 38 39 (147) 

33 ^ 32 

34 ^ 46 26 32 

35 ^ 44 40 41 48 

36 ^ 39 

37 32 

38 ^ 32 33 

39 ^ 32 

40 48 23 

41 ^ 43 32 46 

42 ^ 44 40 41 

Lemma 14.1 Every nontrivial generic orbit intersects Pj for one of the in- 
dices j = 23, 25, 32, 40, 41, 44, 46, 53, 61, 62. 

Proof: Our list of indices is invariant under the involution i — > i + 21, which 
corresponds to the involution of the system, which we have mentioned several 
tiles. For that reason, it suffices to solve the following modified problem. 
We take each number 6 > 42 on the above list and replace it by 6 — 21. 
We when show for every a e [22,42] that a = Oq a^, such that 

Qk E L — {23, 25, 32, 40, 41}. This is equivalent to the original lemma. 
Here is the modified hst. 
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22 ^ 40 23 35 

23 42 37 29 32 26 24 

24 ^ 25 32 37 

25 ^ 31 35 28 27 23 41 40 

26 ^ 25 

27 ^ 42 

28 ^ 25 

29 ^ 23 

30 ^ 23 

31 ^ 32 

32 ^ 25 30 26 36 34 37 42 38 39 (148) 

33 ^ 32 

34 ^ 25 26 32 

35 23 40 41 27 

36 ^ 39 

37 ^ 32 

38 ^ 32 33 

39 32 

40 27 23 

41 ^ 22 32 25 

42 ^ 23 40 41 

Note that a^be L when a is one of 26, 28, 29, 30, 31, 33, 37, 39, 42. But 
27 ^ 42 and 36 ^ 39. So, we can ehminate these as welL Erasing the 
ehminated cases, and also erasing the cases corresponding to elements of L, 
we have 

22 ^ 40 23 35 
24 25 32 37 

34 ^ 25 26 32 (149) 

35 -> 23 40 41 27 
38 ^ 32 33 

Note that 34 — > 26 or 36 — )> 6 G L. Since we have eliminated 26, we 
can eliminate 34. Since 35 ^ 27 or 35 ^ 6 G L, and we have eliminated 
27, we can eliminate 35. Since 22 — > 35 or 22 — > e L, we eliminate 22. 
Since 24 — > 37 or 24 — > 6 e L, and we have ehminated 37, we ehminate 24. 
Since 38 — > 33 or 38 — > 6 G L, and we have ehminated 33, we eliminate 38. d|k 
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14.3 The set B 

The set J3 is divided into 6 layers. These layers correspond to 6 intervals 
J[, Jg, which are just re-indexed versions of the levels corresponding to B. 
That is 

J'l — Ji] J2 = J2] J'3 = J31', J4 = J32', J'5 = Ja] J'q = Jb- (150) 

The intervals J(, Jg are defined by the partition 

< < 20-2 < 1< 2 - < 2 - 20-2 < 2. (151) 

Each layer is decomposed into 4 convex polyhedra, which we call branches. 
The first 3 layers are contained in the bottom half of E, namely the region 

{R/2Z) X {R/2Z) X [0, 1]. (152) 

The remaining 3 layers are the images of the first 3 layers under the involution 
/ discussed above. For this reason, we just list the 12 polyhedra in the first 
3 layers. Our notation is such that Pij is the jth branch of the ith layer. 

14.3.1 Layer 1 





"1 







-7 5 




6 -3 




-7 5 




1 







-7 5 


Bll = 


1 







6 -3 




-2 2 




4 -2 




1 







6 -3 

















5 -3 




2 -1 




2 


-1 




2 -1 



6 -3 

1 

2 -1 



B12 -. 





"1 


■ 






1 









10 


-1 





-4 3 

1 

2 -1 





"1 


o' 






1 


















B13 : 



■9 


-5' 







1 













22 
-3 




-13 
3 




-1 

-3 

-4 3 
3 -1 
2 -1 



9 

-10 
10 



-25 
11 

10 



16 



-25 16 
11 -6 
2 -1 



-12 
3 
2 



22 -13 
-5 4 
2 -1 





"9 -5" 




1 




2 -1 



BU : 



1 

2 
-3 2 





1 







-4 


3 




9 


-5 




-4 


3 






2 







2 










1 







1 






















2 


-1 




2 


-1 





-12 
5 
2 
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14.3.2 Layer 2 



B21 



-7 5 
4 -2 
2 -1 





"1 


" 






1 









to 


-1 





-7 


5 


6 


-3 


2 


-1 



-15 
9 
2 



10 
-5 
-1 



15 10 




6 


-3 


1 -2 




-1 


1 


7 -4_ 




-1 


1 


" 6 - 


3' 


"l 





-1 1 


1 





4 - 


2 


4 


-2 





"l 


■ 




B22 = 


1 









4 


-2 







"l 


■ 






1 









2 


-1 





-4 3 

1 

2 -1 

9 -5' 

-2 2 

2 -1 



B23 : 



6 -3 
-1 1 
4 -2 



-15 
4 
4 



10 
-2 
-2 



-15 10 
4 -2 
7 -4 





"-10 


7 


B24 = 


4 


-2 




4 


-2 



-15 10 
4 -2 
4 -2 



6 


-3 


4 


-2 


-1 


1 



-15 10 
4 -2 
7 -4 



-15 
9 
2 



10 
-5 
-1 



-10 





"6 


-3" 




1 







2 


-1 




15 


10 " 




6 


-3 




2 


-1 



14.3.3 Layer 3 



B31 ■■ 



B32 : 



B33 ■■ 



6 

-1 
9 



CO 




! 2 




-2 




"l o' 




1 




1 





1 
1 
4 -2 



1 
1 
1 



14 
-4 
1 





"6 


-3" 




1 







1 






22 
-12 
9 



-13 



-7 5 
4 -2 
4 -2 



-15 10 
4 -2 
4 -2 





"i ■ 






1 






4 -2 






" 6 - 


3" 




-1 1 






1 





22 
-7 
4 



-13 



-15 10 
4 -2 
1 



-7 


5 


4 


-2 


1 






B34 : 



-15 
4 
4 



10 

-2 
-2 



11 
-2 
4 



-15 10 
4 -2 
1 





"-2 2' 






1 






1 





11 -( 

1 
1 
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14.4 The set A and the Renormalization Map 

The set A is divided into 18 layers, 9 of which are contained in the bottom 
half of E. The involution I maps the first 9 layers to the second 9. For 
this reason, we will just explain the first 9 layers. Again, each layer has 4 
branches, and so we have 36 polyhedra in all. 

Each layer of A corresponds to an interval of the form 

Iij = IinR-\j;) (153) 

Here Ji, Is is the partition relative to which the circle renormalization map 
R is defined. These intervals are defined by the partition 

< < 20-2 < 2 - < 2 - 20-2 < 2. (154) 

The notation Aijk denotes the A;th branch of the {i, j)th layer of A. Rather 
than list the coordinates of the Aijk, we will define these polyhedra in terms 
of the renormalization map R. This approach simultaneously defines A and 
R. The map R : A ^ B is a, piecewise golden affine transformation of R^. 
We encode golden affine transformations by 8-tuples of integers. The 8-tuple 
(ai, as) corresponds to the map T{V) — rV + W, where 

r = 01 + 020; 1^ = (03 + O40, O5 + O60, O7 + O80). (155) 

The maps we actually hst are R"^. 
There are 6 special cases. 



• 


^314 


— R ^{Bii) 


, where R~^ — 


(-3,2,4,-2,8, -4,4,-2). 


• 


^364 


= R~^{Bq4) 


, where R~^ = 


(-3,2,4,-2,0,0,4,-2). 


• 


Auk 


= R-\Bis) 


, where R~^ = 


(-3,2,4,-2,8,-4,0,0) for k 


• 


A^ek 


= R~^(Be2,) 


, where R~^ — 


(-3,2,4,-2,0,0,8,-4) for k 


Otherwise, 


we have 






1. 


Aijk 


= R~^(Bjk) 


, where R~^ — 


(-3,2,4,-2,4,-2,0,0) 


2. 


A2jk 


= R~^{Bjk) 


, where R~^ = 


(1,0,0,0,0,0,2,-2) 


3. 


Asjk 


— R~^{Bjk) 


, where, R~^ = 


(-3,2,4,-2,4,-2,4,-2). 


4. 


A^jk 


— R'~^{Bjk) 


, where, R~^ — 


(1,0,0,0,0,0,-2,2) 


5. 


A^jk 


= R~^{Bjk) 


, where, R~^ — 


(-3,2,4,-2,4,-2,8,-4) 
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